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A3STRACT:  A  computer  program  has  been  developed  at  NOL  that  describes 
the  shock  and  blast  loading  characteristics  of  the  detonation  of  a 
high  explosive  projectile  internal  to  an  aircraft  structure;  both 
shock  wave  and  confined-explosion  gas  pressure  loads  are  considered. 
With  certain  modifications,  the  program  can  be  made  applicable  to  any 
internal  explosion  irrespective  of  the  type  of  confining  configuration, 
e.g.,  a  naval  ship  compartment,  land  vehicle,  or  building  structure. 
Discussions  are  given  on  the  general  use  and  content  of  the  program, 
the  input  options  available  in  the  code,  and  the  technical  aspects  of 
the  calculational  methods  used  to  determine  shock  loading  functions, 
confined-explosion  gas  pressure,  venting  of  the’  confined  gases,  and 
damage  propagation  to  other  areas  of  the  aircraft.  Comparisons  of  code 
results  with  available  experimental  data  are  presented  to  demonstrate 
the  Justifiable  confidence  in  the  use  of  the  code  on  aircraft  problems. 
Complete  documentation  of  the  code  is  given  together  with  results  of 
sample  problems  that  show  the  various  features  of  the  code  and  the 
readily  usable  form  of  the  resultant  loading  information. 
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00011).  The  objective  of  this  task  was  to  develop  a  computer  program 
for  describing  blast  characteristics  associated  with  the  detonation 
of  a  high  explosive  projectile  internal  to  an  aircraft  structure. 
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SUMMARY  AND  CONCLUSIONS 

Assessment  of  damage  to  aircraft  structures  from  the  detonation 
of  explosive  projectiles  internal  to  the  aircraft  requires  a  detailed 
knowledge  of  the  dynamic  pressure  leads  applied  to  various  structural 
elements.  NOL  has  developed  a  computer  program  that  is  capable  of 
generating  characteristic  blast  loading  parameters  associated  with 
confined  explosions  in  a  form  readily  usable  by  aircraft  design 
engineers  and  vulnerability  analysts.  Existing  state-of-the-art 
explosion  theory  and  experimental  data  were  used  as  the  basis  for  the 
shock  wave  calculations  available  in  the  code.  An  improved  method 
of  predicting  the  confined-explosion  gas  pressure  that  exists  after 
shock  dissipation  was  developed  especially  for  this  code.  Any  size 
explosion  can  be  treated  by  the  code  for  any  ambient  altitude  condition 
up  to  and  above  50,000  ft,  and  the  code  includes  the  blast  properties 
of  some  29  different  explosives  including  mono,  composite,  and  alu¬ 
minized  varieties. 

The  computer  program  analytically  divides  the  internal  explosion 
Into  two  damaging  mechanisms — the  shock  wave  and  the  confined- 
explosion  gas  pressure.  For  the  shock  wave  it  generates  the  incident 
and  normally  reflected  pressure-time  histories  and  impulses  for  the 
positive  phase  duracion  at  a  specified  distance  from  the  explosion. 
Existing  data  and  theory  were  used  to  develop  the  shock  calculational 
model.  The  code  reduces  the  shock  calculation  for  all  cases  to  the 
reference  data  from  a  free-field,  bare,  spherical  1-lb  TNT  explosion. 
Variables  that  affect  airblast  which  are  included  in  the  code  for 
establishing  an  equivalent  TNT  spherical  explosion  are  (1)  explosive 
weight,  (2)  type  of  explosive,  (3)  cylindrical  charge  geometry, 

(4)  case  weight  of  the  projectile,  and  (5)  ambient  pressure  and 
temperature  at  the  location  in  the  aircraft  where  the  explosion  occurs. 
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'for  an  explosion  internal  to  a  confining  structure  or  compartment, 
a  long-duration  quasi-static  pressure  exists  after  dissipation  of 
the  shock  wave.  The  maximum  value  of  the  pressure,  defined  as  the 
confined-explosion  gas  pressure,  is  dependent  on  these  parameters; 

(1)  weight  of  explosive,  (2)  type  of  explosive  (chemical  composition), 
(3)  volume  of  compartment,  and  (4)  pressure  and  temperature  of  the 
air  initially  in  the  compartment.  Because  of  the  inadequacies  of 
existing  methods  of  calculating  the  confined-explosion  gas  pressure, 
a  technique  was  developed  especially  for  this  program  that  follows 
the  energy  generation  of  the  chemical  reactions  and  the  changes 
in  gas  properties  as  the  confined-explosion  gas  pressure  is  developed. 
In  a  completely  closed  compartment  or  structure,  heat  loss  to  the 
surrounding  walls  is  the  only  mechanism  for  reducing  the  pressure  in 
time,  but  this  phenomenon  is  neglected  in  this  program  because  of  the 
very  long  durations  involved.  However,  for  the  aircraft  structure, 
there  will  be  openings  or  vent  areas  through  which  the  confined  gases 
can  escape  such  as  the  initial  opening  due  to  entry  of  the  projectile 
into  the  compartment  and  any  fragment  penetration  openings.  Also  the 
pressure  can  change  abruptly  due  to  wall  failure  of  the  compartment 
which  introduces  a  new  compartment  volume.  The  computer  program 
calculates  the  variation  of  the  confined-explosion  gas  pressure  with 
time  for  venting  and  such  volume  changes.  Vent  area  and  volume 
changes  are  controlled  by  input  damage  criteria  for  compartment  wall 
failure. 

In  addition  to  the  technical  description  of  the  calculational 
models  contained  in  the  computer  program,  a  user’s  guide  and  complete 
documentation  of  the  code  are  given  in  the  text  of  the  report  and  the 
attached  appendices.  Also  nine  sample  problems  are  presented  that 
demonstrate  the  many  options  and  features  of  the  program. 

Although  no  single  set  of  experimental  data  was  available  to 
compare  with  the  overall  code  performance,  each  individual  calcula¬ 
tional  model  was  tested  against  pertinent  experimental  data.  Suf¬ 
ficient  shock  and  confined-explosion  gas  pressure  data  were  found 
for  comparison,  and  the  agreement  with  code  predictions  in  these 
cases  was  excellent.  It  is  concluded  that  the  calculational  models 
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for  these  most  important  aspects  of  the  internal  blast  loading  can 
be  used  with  justifiable  confidence.  Whereas  the  pressure-time 
decay  of  the  confined-explosion  gas  pressure  due  to  venting  has 
been  verified  with  limited  data,  the  introduction  of  volume  changes 
has  not  been  tested.  The  method  for  treating  instantaneous  volume 
changes  is  based  on  fundamental  thermodynamic  relations,  thus  there 
is  no  reason  to  believe  that  this  section  of  the  calculation  detracts 
from  the  use  of  the  code  for  general  aircraft  internal  blast  problems. 

The  computer  program  has  wide  range  potential  for  use  in  studies  of 
structural  response  of  any  military  or  civilian  system  to  an  internal 
explosion  be  it  aircraft,  naval  ship,  land  vehicle,  or  building 
structure.  Although  it  is  adequate  for  the  aircraft  problems  for 
which  it  is  presently  designed,  there  are  five  areas  in  the  code  that 
require  additional  study  and  possible  modifications  before  its 
generality  can  be  claimed  for  large  explosions  in  large  structure 
compartments  such  as  ship  compartments  and  building  rooms.  These  are 
(1)  multiple  shock  reflections  from  surrounding  walls,  (2)  heat 
losses  to  surrounding  walls  that  might  reduce  the  confined-explosion 
gas  pressure  to  a  significant  degree  for  large  structures,  (3)  variable 
backpressure  to  the  venting  process,  (4)  gas  pressure-time  history 
where  mixing  of  gases  after  wall  failure  occurs  in  a  finite  time 
interval,  and  (5)  subsequent  chemical  reactions  with  the  air  in  adja¬ 
cent  compartments  after  wall  failure  if  complete  combustion  is  not 
achieved  in  the  initial  compartment.  With  modification  to  the  code 
reflecting  the  above  suggested  studies,  it  is  believed  this  computer 
code  could  evolve  as  a  general  service-wide  tool  for  the  investigations 
of  structural  response  to  internal  explosion  loading. 
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CHAPTER  1 


INTRODUCTION 


One  of  the  current  ongoing  tasks  of  the  Aerial  Target  Vulnera¬ 
bility  (ATV)  Program  of  the  Joint  Technical  Coordinating  Group  for 
Munitions  Effectiveness  ( JTCG/ME)  has  been  the  development  of  a 
component  damage  data  bank.  An  item  defined  in  this  data  bank  is  the 
vulnerability  of  aircraft  to  internal  blast  from  high  explosive 
projectiles.  Under  the  direction  of  the  Survivability-Vulnerability 
Branch,  Prototype  Division,  Air  Force  Flight  Dynamics  Laboratory 
(AFFDL),  the  objectives  of  this  task  are  (1)  to  define  the  internal 
blast  loading  characteristics  from  a  high-explosive  projectile, 

(2)  to  determine  the  damage  to  aircraft  structures,  and  (3)  to  assess 
the  vulnerability  of  these  structures  to  internal  blast  effects. 

The  Naval  Ordnance  Laboratory  (MOL)  was  assigned  the  technical 
solution  of. the  first  task  problem  area,  namely,  to  define  the 
internal  blast  loading  characteristics. 

Specifically,  the  objective  of  the  NOL  program  was  to  develop 
mathematical  and  graphical  techniques  for  describing  blast  character¬ 
istics  associated  with  the  detonation  of  a  high  explosive  projectile 
Internal  to  an  aircraft  structure.  Existing  state-of-the  art 
experimental  data  and  explosion  theory  were  to  be  combined  using 
sound  engineering  judgment  to  provide  a  computer  program  capable 
of  generating  characteristic  shock  wave  and  blast  loading  from  an 
explosion  internal  to  an  aircraft.  Execution  of  the  computer  program 
and  the  resultant  loading  functions  were  to  be  in  a  form  readily 
usable  by  aircraft  design  engineers  and  vulnerability  analysts. 

Although  this  task  was  directed  to  the  solution  of  the  aircraft 
problem,  the  concepts,  content,  and  format  of  the  resultant  computer 
code  can  be  related  to  any  military  or  civilian  system  be  it  aircraft, 
naval  ship,  land  vehicle,  or  building  structure.  The  code  was 
structured  to  accommodate  easy  modification  for  any  new  system. 
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With  usage  on  response  problems  for  structures  other  than  aircraft, 
it  is  hoped  that  a  more  complete  internal  blast  loading  computer 
program  will  evolve  for  general  use. 
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CHAPTER  2 

GENERAL  DESCRIPTION  AND  LIMITATIONS  OF 
COMPUTER  PROGRAM 

Assume  that  a  high  explosive  is  detonated  in  a  closed  structure 
of  some  arbitrary  geometry  with  a  small  vent  opening.  If  a  pressure 
sensor  were  to  be  placed  on  the  wall  of  the  structure,  it  would 
indicate  a  pressure-time  history  of  the  type  shown  in  Figures  2.1(a) 
and  (b).  On  an  expanded  time  scale  (a),  one  would  note  the  initial 
peak  reflected  shock  overpressure,  APr,  followed  by  subsequent 
reflected  shock  pulses  from  the  adjacent  confining  wall  of  the 
structure.  The  oscillations  would  dissipate  leaving  a  quasi-static 
overpressure,  AP_,  created  by  the  heated  gases  contained  in  the 
structure;  this  pressure  is  defined  as  the  confined-explosion  gas 
pressure.  On  a  reduced  time  scale  (b),  the  shock  reflections  would 
appear  as  high  spikes  near  time  zero.  The  confined-explosion  gas 
pressure,  AP_,  would  be  clearly  established  on  this  time  scale. 

Even  for  a  completely  closed  structure,  the  gas  pressure  would 
slowly  decrease  in  time  due  to  heat  losses  to  the  surrounding 
structure  walls.  For  a  structure  with  some  openings  through  which 
venting  could  occur,  the  gas  pressure  would  decay  much  more  rapidly. 

An  accurate  description  of  the  pressure-time  history  during  the 
multiple  reflected  shock  phenomena  in  a  closed  structure  of  arbitrary 
configuration  was  far  beyond  the  scope  of  this  program  effort  and 
economically  beyond  the  scope  of  any  three-dimensional  hydrodynamic 
code.  For  this  reason  the  shock  wave  calculations  in  this  program 
are  limited  to  the  incident  and  normally  reflected  pressure-time 
shock,  depicted  in  Figure  2.1  (c),  arriving  initially  at  a  point 
on  the  structure  wall.  It  is  believed  that  the  normally  reflected 
pressure-tine  shock  history  and  associated  reflected  impulse  is 
sufficient  to  provide  a  meaningful  index  in  determining  the  local 
structural  response  to  shock  wave  loading.  Further  it  is  believed 
that  the  predominant  damaging  mechanism  from  an  internal  explosion 
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Is  the  confined-explosion  gas  pressure.  For  most  applications,  the 
structure  wall  can  be  treated  as  though  it  were  given  an  initial 
velocity  by  the  shock  wave  as  an  initial  boundary  condition.  Sub¬ 
sequent  loading  on  the  structure  is  defined  by  the  confir.ed-explosicn 
gas  pressure  which  is  handled  as  a  separate  loading  phenomenon 
completely  decoupled  from  the  shock  wave. 

The  initial  magnitude  of  the  confined-explosion  gas  pressure 
is  determined  from  a  tecnnique  developed  specifically  for  this 
program,  which  will  be  described  in  detail  in  a  subsequent  chapter. 
Relative  to  the  slower  plastic  response  time  of  a  typical  aircraft 
structure,  it  is  assumed  that  the  confined-explosion  gas  pressure, 

AP  ,  depicted  in  Figure  2.1  (d)  is  developed  instantaneously  in  time. 

o 

In  other  words,  the  chemical  reaction  of  the  explosion  gas  products 
with  the  surrounding  air  in  the  initial  confines  of  the  structure 
and  the  heat  transfer  to  the  resultant  gas  mixture  are  assumed  to 
Occur  instantaneously  to  develop  the  confined-explosion  gas  pressure. 
This  pressure  decreases  in  time  due  to  venting  through  available 
openings  created  by  the  initial  entry  of  the  projectile  into  the 
structure,  subsequent  openings  from  fragment  penetrations,  and  any 
normal  structural  openings  such  as  cable  passageways.  Venting  cal¬ 
culations  assume  a  constant  back  pressure  equal  to  atmospheric 
conditions  outside  the  aircraft  since  most  leakage  would  occur 
through  fragment  penctx*ations  in  the  aircraft  skin.  Heat  losses  to 
the  structure  walls  that  could  reduce  the  gas  pressure  are  neglected 
because  (1)  significant  heat  loss  would  require  times  much  larger 
than  the  plastic  response  times  of  the  structures  associated  with  the 
typical  small  aircraft  compartments  and  (2)  pressure  decreases  much 
more  rapidly  from  venting  through  even  the  small  projectile  penetra¬ 
tion  opening  than  from  heat  losses. 

Provisions  are  made  to  account  for  sudden  changes  in  gas 
pressure  due  tc  structural  failure  resulting  in  rapid  expansion  into 

an  adjacent  compartment.  It  is  assumed  for  these  calculations  that 

/ 

the  change  in  pressure  is  instantaneous  relative  to  the  much  slower 
plastic  response  times  of  the  aircraft  structures.  Also  in  general 
the  small  compartment  sizes  in  an  aircraft  structure,  would  indicate 
a  rapid  stabilization  of  pressures  if  comportment  walls  failed. 
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In  summary  to  this  point,  four  assumptions  have  been  made  that 
may  limit  the  use  of  this  computer  code  to  aircraft  only.  They  are.: 
(1)  multiple  shock  reflections  are  neglected — only  the  initial  shock 
serves  as  a  damage  index,  (2)  no  variation  in  venting  back  pressure 
occurs,  (3)  heat  losses  to  structure  walls  are  neglected,  and  (4) 
instantaneous  change  in  pressure  occurs  with  compartment  wall  failure 
propagation.  Whereas  it  is  believed  that  these  assumptions  do  not 
significantly  restrict  the  study  of  aircraft  response  to  Internal 
blast,  general  application  of  this  computer  program  directly  to  other 
structures,  such  as  ships  and  buildings,  may  be  hampered  by  these 
assumptions.  Therefore,  flexibility  in  code  construction  has  been 
provided  to  allow  for  easy  and  efficient  modification  to  those 
sections  that  would  be  affected  by  alterations  to  these  assumptions. 

The  basis  for  this  entire  study  was  existing  state-of-the-art 
theory,  analytical  methods,  and  experimental  data  for  explosions. 

When  directed  to  the  problem  of  internal  blast,  these  in  themselves 
introduce  limitations  In  terms  of  applicability,  and  some  are  open 
to  interpretation  even  by  explosion  experts.  Since  the  prime  users 
of  this  code  probably  will  not  be  people  with  background  in  explosion 
effects,  all  pertinent  explosion  properties  are  self-contained  in  the 
code.  Only  the  type  and  amount  of  explosive  are  required  as  input 
to  the  code.  In  calculations  where  limitations  arising  from  theory 
and  data  deficiencies  are  encountered,  caution  statements  are  clearly 
indicated  in  the  output  statements  for  the  user's  benefit. 

The  next  four  chapters  present  the  technical  aspects  of  the 
caiculatlonal  methods  contained  in  the  code  including  theoretical  and 
experimental  background,  certain  operational  procedures,  and  explana¬ 
tions  of  the  more  important  features  and  options  in  the  code.  Com¬ 
parisons  of  code  predictions  with  available  experimental  data  are 
also  given  in  these  sections.  Chapter  7  and  the  attached  appendices 
represent  a  user's  guide  or  manual  for  the  detailed  content  and 

1 

operation  of  the  code.  Also  a  number  of  sample  problems  are  included 
in  this  code  documentation  section  to  acquaint  the  user  with  the 
various  options  available  in  the  code. 
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TYPICAL  TRACES  FROM  ACTUAL  INTERNAL  EXPLOSION 


(«)  (b) 


CODE  APPROXIMATIONS  FOR  INTERNAL  EXPLOSION 


FIG.  2.1  TYPICAL  PRESSURE-TIME  CURVES  FOR  AN  INTERNAL  EXPLOSION 
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CHAPTER  3 

INPUT  DATA  REQUIREMENTS  / 

Explosive  Parameters.  For  a  typical  high  explosive  projectile^, only 
four  of  its  properties  are  required  as  input  to  the  computer  program. 
They  are:  (1)  weight  of  explosive,  (2)  type  of  explosive  (3)  length 
to  diameter  ratio  of  charge,  and  (4)  metal  case  weight  to  charge 
weight  ratio.  Incorporated  in  the  code  are  the  pertinent  properties 
of  24  types  of  explosives  (Ref.  (1)).  Table  3»1  gives  the  coded 
properties  of  these  24  explosives  plus  3  mono  explosives  that  are 
rarely  used  alone  as  the  main  charge.  The  properties  of  aluminum 
and  a  common  wax  binder  are  also  added.  If  the  desired  explosive 
•  is  contained  in  this  table,  it  is  necessary  only  to  input  the  index 
number  to  specify  the  explosive  type.  If  one  wishes  to  input  an 
explosive  not  j.a  the  table,  an  index  number  of  0  is  used  and  the 
required  explosive  properties  must  be  specified.  However,  if  for 
some  reason  the  energy  equivalent  weight  is  not  known,  a  zero  for 
this  quantity  will  permit  shock  calculations  to  be  made  with  an 
equivalent  weight  of  one — the  same  as  for  TNT.  A  diagnostic  statement 
will  appear  in  the  output  —  "WF ACT  IS  NOT  KNOWN,  1.0  IS  USED5'— which 
means  the  shock  calculations  are  equal  to  those  for  TNT.  If  the 
desired  explosive  is  a  mixture  of  components  in  Table  3.1,  an  index 
number  of  -1  is  used  and  the  weight  fraction  of  each  component  must 
be  specified.  Again  the  energy  equivalent  weight  may  not  be  known, 
but  it  can  be  handled  in  the  same  manner  as  before  by  letting  the 
equivalent  weight  equal  zero.  The  only  restriction  on  the  type  of 
explosive  used  in  this  program  is  that  the  explosive  must  be  of 
C-H-N-0  form  with  aluminum  as  the  only  possible  metallic  additive. 

It  should  be  noted  that  this  restriction  arise-s  from  the  confined- 
explosion  gas  pressure  calculations;  and  as  it  will  be  discussed  in 
subsequent  chapters,  this  restriction  has  no  obvious  theoretical 
basis  but  must  bs  invoked  because  of  the  lack  of  experimental  data 
on  other  metallic  additives  or  non 
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The  computer  program  is  capable  of  making  corrections  for 
cylindrical  charge  shape  factors  for  length  to  diameter  ratios 
(L/D)  between  2  and  10.  This  is  sufficiently  general  to  accommodate 
most  common  anti-aircraft  projectiles.  If  the  L/D  ratio  is  less 
than  2  as  input,  the  code  will  treat  the  charge  as  spherical. 

To  account  for  effects  of  metal  casing  on  the  degradation  of 
the  shock  wave,  it  is  necessary  to  input  the  case  weight  to  charge 
weight  ratio  (M/C)  for  the  weapon.  The  total  case  weight  (case, 
nose,  fins,  and  fuze)  should  not  be  used  in  determining  the  ratio. 
Rather,  it  is  recommended  that  only  the  case  weight  immediately 
adjacent  to  the  explosive  charge  in  the  radial  direction  be  used. 


Initial  Conditions  in  Structural  Compartment.  In  order  to  calculate 
internal  blast  characteristics,  it  is  necessary  to  specify  the 
initial  geometric  properties  of  the  structural  compartment  in  which 
the  explosion  occurs  and  of  the  air  that  Is  confined  in  the  compartment 
Specifically  for  this  computer  program,  they  are  (1)  initial  ambient 
pressure  and  temperature  of  the  confined  air  (air  is  treated  con¬ 
ventionally  as  79 %  N2  and  2\%  C>2  by  volume),  (2)  initial  gas  volume 
of  the  compartment,  (3)  initial  vent  area  for  confined  gases  (would 
Include  opening  from  projectile  entry  and  additional  openings  frpm 
fragment  penetrations),  and  (4)  ambient  pressure  or  backpressure 
against  which  venting  would  occur,  I.e.,  air  pressure  outside  the 
compartment.  If  the  ambient  conditions  of  items  (1)  and  (4)  are  the 
same  as  those  for  air  at  the  altitude  at  which  the  aircraft  is 
located,  only  the  altitude  of  the  aircraft  need  be  specified.  For 
this  case  the  computer  code  has  the  1959  ARDC  standard  atmosphere 
taken  from  reference  (2)  as  a  subroutine  for  determination  of  the 
initial  ambient  pressure  and  temperature. 


Shock  Calculations.  To  make  any  shock  calculation,  the  distance 
from  the  point  of  detonation  to  a  desired  location  in  the  compartment 
must  be  designated.  Specifically  for  this  problem,  this  distance 
is  measured  radially  from  the  point  of  detonation  to  the  point  on 
the  structure  wall  where  shock  pressure-time  information  is  desired. 
As  input  to  the  code  the  total  number  of  different  distances  of 
interest  must  be  specified,  followed  by  a  list  of  these  desired 
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distances.  In  this  manner  the  variation  of  shock ^loading  as  a 
function  of  location  on  a  particular  structural  wall  can  be 
determined.  There  are  two  options  in  the  code  that  must  be  specified 
as  input  depending  on  the  type  of  shock  and  confined-explosion  gas 
calculations  desired.  For  normal  code  operation  where  both  shock 
and  confined-explosion  gas  pressure  calculations  are  of  interest, 
option  1  is  used  with  the  number  and  list  of  distances.  However, 
if  shock  calculations  are  not  desired,  then  option  1  is  used  and  the 
number  of  distances  is  set  to  0,  and  the  shock  calculation  section 
of  the  program  is  bypassed.  On  the  other  hand,  if  one  wishes  to 
examine  only  shock  calculations,  option  2  is  used  with  the  number  and 
list  of  distances,  and  the  confined-explosion  gas  pressure  section 
of  the  program  is  bypassed. 

Volume  and  Vent  Area  Changes.  It  is  quite  possible  that  damage  to 
an  aircraxt  from  an  internal  projectile  explosion  will  propagate 
beyond  the  confines  of  the  initial  compartment  where  detonation 
occurs.  Excessive  shock  loading  or  confined-explosion  gas  pressure 
may  fail  a  compartment  wall  allowing  the  confined  gases  to  propagate 
to  an  adjacent  compartment.  If  this  pressure  remains  excessive, 
additional  wall  failures  may  occur  with  the  subsequent  spread  of 
the  confined  gases.  As  the  gases  propagate  to  different  compartments 
and  occupy  larger  volumes,  the  pressure  is  reduced.  When  the  confined- 
explosion  gas  pressure  has  decreased  to  the  point  where  wall  failure 
does  not  occur,  the  damage  propagation  stops.  Although  a  structural 
response  code  will  eventually  be  interfaced  with  this  blast  loading 
code  to  assess  the  damage  propagation  to  these  pressure  loads,  a 
skeleton  format  is  provided  in  this  code  to  allow  for  an  initial 
examination  of  the  propagation  phenomena.  As  an  example,  take  the 
box  structure  represented  in  Figure  3.1  where  the  initial  explosion 
occurs  in  the  compartment  with  the  circled  X.  Through  wall,  failure 
(sides  with  cross-lines),  the  damage  may  propagate  to  compartments  A, 

D,  and  C.  It  is  necessary  to  specify  certain  conditions  that  control 
the  damage  propagation  such  as  wall  failure  criteria  and  the  amount 
and  condition  of  the  air  in  the  adjacent  compartments.  There  are 
three  options  available  in  the  code  to  specify  desired  characteristics 
of  compartment  wall  failure. 
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The  table  in  Figure  3.1  is  an  example  of  failure  criteria 
option  3.  Interpretation  of  this  input  table  to  the  computer  program 
is  as  follows.  (1)  If  the  confined-explosion  gas  pressure  in  the 
'initial  compartment  is  above  ^5  psia  0.15  sec  after  detonation,  the 
wall  fails  allowing  the  gases  to  mix  with  4  ft^  of  air  at  1^.7  psia 
and  20°C  in  compartment  A  and  providing  an  arbitrary  additional  vent 

P 

area  of  0.005^5  ft  .  (2)  If  the  gas  pressure  after  mixing  with  air 

in  compartment  A  is  above  20  psia  0.60  sec  after  detonation,  the 
next  wall  fails  exposing  compartment  B  which  has  h  ft^  of  14.7  psia, 

p 

20°C  air  and  an  additional  arbitrary  vent  area  of  0.005^5  ft  . 

(3)  If  the  gas  pressure  after  mixing  in  compartment  B  is  above  19  psia, 
wall  failure  occurs  involving  compartment  C,  etc.  Option  3  offers 
both  pressure  and  response  time  control  on  damage  criteria.  If  at  any 
step  the  pressure  is  below  the  tabulated  value  at  the  specified  time, 
propagation  stops.  For  example,  if  the  pressure  is  below  20  psia 
at  0.60  sec  after  detonation,  damage  does  not  propagate  to  compart¬ 
ments  B  or  C. 

Damage  criteria  options  1  and  2  are  simplified  versions  of 
the  above.  Option  1  specifies  only  the  pressure  failure  levels  with 
wall  failures  occurring  instantaneously  in  time.  Option  2  specifies 
only  the  time  of  failure  irrespective  of  pressure  level.  Both  of 
these  options,  like  option  3,  require  as  input  the  volume  and  ambient 
conditions  of  the  air  in  the  various  compartments  and  any  additional 
vent  area 
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INDEX 

NUMBER 


EXPLOSIVE 

NAME 

EQUIVALENT 

WEIGHT 

f 

e 

TNT 

1.00 

TNETB 

1.13 

EXPLOSIVE  D 

0.85 

PENTOLITE 

1.17 

PICRATOL 

0.90 

CYCLOTOL 

1.14 

COMP  B 

1.10 

RDVWAX 

98/2 

1.19 

COMP  A-3 

1.09 

TNETB/AL 

90/10 

1.23 

TNETB/ AL 

78/22 

1.18 

TNETB/ AL 

72/28 

1.18 

TNETB/AL 

65/35 

1.23 

TRITONAL 

1.07 

RDX/Al/WAX 

88/10/2 

1.30 

RDX/Al/WAX 

78/20/2 

1.32 

RDX/Al/WAX 

74/21/5 

1.30 

RDX/AL/WAX 

74/22/4 

1.30 

RDX/Al/WAX 

62/33/5 

1.19 

TORPEX  II 

1.24 

H*6 

1.27 

H8X-I 

1.21 

H8X-3 

1.16 

TNET8/RDX/AL 

39/26/35 

1.24 

ALUMINUM 

0 

WAX 

0 

RDX 

0 

PETN 

0 

TETRYL 

0 

TABLE  3.1 

LIST  OF  EXPLOSIVE  PROPERTIES 
(BASED  ON  DATA  FROM  REF.  1) 


HEAT  OF 
FORMATION 
(CAL/GM) 


-78.40 

-307.1 

-386.3 

-242.8 

-238.5 

22.79 

4.33 

57.00 

24.93 

-276.4 

-239.5 

-221.1 

-199.6 

-62.72 

50.38 


-3.57 

-12.56 

-22.93 

-21.83 

-102.6 

0 

-392.0 

66.16 

-407.1 

16.26 


WEIGHT  FRACTIONS  OF 
COMPONENTS 


H  N 


0.170  10.452  I  0.220 
0.157 


8  |  0.148 
0.333 


0.379  j  0.4J2 
0.177 
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CHAPTER  H 


SHOCK  WAVE  CALCULATIONS 


Base  Data.  The  principal  thesis  of  the  shock  wave  calculations  ir. 
this  computer  program  is  that  a  cased,  cylindrical  charge  of  a  given 
type  and  amount  of  explosive  detonated  at  any  altitude  from  sea  level 
to  at  least  50,000  ft  can  be  equated  to  a  free-fiold  1-lb  TNT 
spherical  explosion  at  sea  level.  Generally  explosion  data  given  in 
handbooks  for  TNT  do  not  provide  sufficient  information  to  yield  the 
pressure-time  history  of  the  shock  at  a  specified  distance  from  the 
explosion;  usually  only  peak  pressures,  positive  phase  durations,  and 
positive  impulses  are  given.  One  must  turn  to  various  hydrodynamic 
codes  to  obtain  such  information  in  lieu  of  extensive  experimental 
data.  However,  such  codes  are  lengthy  and  expensive  to  run  and  do  not 
lend  themselves  to  the  objective  of  this  program.  It  was  decided 
to  use  results  from  a  current  version  of  the  WUNDY  hydrocode  developed 
at  NOL  and  described  in  reference  (3),  to  normalize  these  results  to 
form  a  family  of  pressure-time  curves  for  a  large  number  of  distances, 
and  to  find  the  best  empirical  fit  to  represent  these  results. 

Figure  *1.1  shows  four  representative  curves  developed  by  WUNDY 
that  demonstrate  the  incident  pressure-time  behavior  of  a  free-field 
shock  wave  as  a  function  of  distance,  R.  From  some  25  curves  of  this 
type  that  exist  outside  the  explosion  gas  contact  surface,  it  was 
found  that  the  family  could  be  represented  quite  well  by  the  equation. 


tt  =  AP/APi  =  (1  -  t)  e' 


<1  + 


(4.1) 


where 


(t  -  ta)/td 


a  =  (228/R)  -  0.95 


A  =  0.5 
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and  AP^  *  peak  incident  shock  overpressure 

AP  =  instantaneous  overpressure 

t  e  time  measured  from  detonation 

t  *  arrival  time  of  shock  measured  from  detonation 
a 

t^  =  positive  phase  duration  of  incident  shock  pulse 
R  -  distance  from  detonation  (cm) 

From  the  above  equations,  it  is  seen  that  peak  incident 
pressure,  arrival  time,  and  positive  phase  duration  for  a  given 
distance  are  the  only  parameters  required  for  development  of  the 
shock  pressure-time  curve.  Values  of  peak  incident  overpressure  ana 
arrival  time  are  readily  available  from  WUNDY  code  results  and  have 
been  tabulated  in  the  computer  program  for  108  distances  ranging 
from  the  charge  surface  to  2.342  x  10^  cm  (based  on  l~lb  bare TNT  sphere ). 
Positive  phase  durations  for  distances  outside  the  contact  surface 
are  also  obtained  from  WUNDY.  However,  inside  about  70  cm,  the 
positive  phase  duration  of  the  shock  wave  is  not  completed  before 
interaction  with  the  contact  surface  occurs.  Although  WUNDY  follows 
the  contact  surface  boundary,  the  available  runs  do  not  yield  usable 
results  inside  the  contact  surface.  Therefore,  experimental  data 
from  references  (4)  and  (5)  were  used  to  derive  approximate  positive 
phase  durations  Inside  the  contact  surface.  The  computer  code  assumes 
that  equation  (4.1)  continues  to  hold  inside  the  contact  surface. 

Arrival  time  of  the  contact  surface  is  programmed  into  the  code.  When 
shock  information  is  desired  at  a  point  inside  the  contact  surface, 
the  user  is  alerted  to  the  fact  that  the  pressure-time  history  is  an 
approximation  by  a  diagnostic  or  warning  statement  that  appears  in 
the  code  output— "CAUTION— CONTACT  SURFACE  HAS  ARRIVED.  DATA  ARE 
CRUDE  BEYOND  T(KSEC)  AFTER  SHOCK  ARRIVAL  =". 

As  indicated  previously,  parameter  values  are  tabulated  in  the 
computer  code  for  108  distances.  An  interpolation  method  was  needed 
to  accommodate  any  given  distance.  Plots  of  peak  pressure,  shock 
arrival  time,  positive  phase  duration,  and  contact  surface  arrival 
time  as  functions  of  distance  on  log-log  scales  demonstrated  a  nearly 
linear  slope  over  relatively  small  intervals.  Therefore,  linear 
interpolation  between  the  log  values  of  the  parameters  and  the  log 
values  of  distance  is  coded  into  the  program  as  an  accurate  interpo¬ 
lation  method.  1A 
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Scaling  Equations.  To  relate  a  spherical  TNT  explosion  at  altitude 
to  a  1-lb  spherical  TNT  explosion  at  sea  level,  conventional  Sachs 
scaling  was  used  in  the  computer  program.  (Sachs  scaling  method  can 
be  found  in  many  references  on  airblast  from  explosions,  such  as 
reference  (6).)  The  scaling  relations,  as  they  are  used  in  this 
computer  code,  are  given  as 

Rs  -  Ra(ws/wa)1/3  (Pa/Ps)1/3  (4*2) 

APa  =  APs(W  (4‘3) 

fca  *  ts(Wa/Ws)1/3  (Ps/Pa)1/3  (Ts/Ta)1/2  (**.4) 

Ift  *  Is(Wa/Ws)1/3  (pa/ps>2/3  (Ts/Ta)1/2  (4.5) 

where 

R  *  distance 
W  «  charge  weight 
AP  *  overpressure 
P  *=  ambient  pressure 
t  =  time 

T  =  ambient  temperature 
I  =  impulse 

s  =  subscript  denoting  1-lb  TNT  sphere  at  sea  level 
a  ■  subscript  denoting  TNT  sphere  at  altitude 

Following  the  order  of  these  equations,  for  a  given  distance, 

Ra,  from  a  spherical  TNT  charge,  W  ,  at  altitude  pressure,  P&, 
and  temperature,  Tg;  a  scaled  distance,  Rg,  from  a  1-lb  TNT  spherical 
charge  (W  «1)  at  sea  level  pressure  and  temperature,  Po  and  T  ,  is 
determined.  The  code  then  calculates  a  pressure-time  curve  for  the 
scaled  distance  and  1-lb  TNT  charge  at  sea  level  from  equation  (4.1) 
and  numerically  integrates  this  curve  to  determine  the  incident  impulse. 
With  values  P  ,  t_,  and  I,  equations  (4.3)— (4.5)  define  the  values 
of  these  parameters  at  the  desired  altitude.  More  details  of  this 
method  will  be  presented  in  subsequent  sections. 
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Equivalent  Weight.  In  the  previous  section  on  scaling,  all  charges 
were  TNT  spheres.  Since  typical  projectile  charges  are  cased, 
cylindrical,  non-TNT  explosives,  methods  for  equating  the  blast 
effects  of  a  real  projectile  to  those  of  an  idealised  TNT  sphere 
were  required.  Virtually  all  studies  directed  at  establishing 
equivalent  weights  have  been  conducted  at  sea  level  conditions.  The 
assumption  was  made  for  this  program  that  the  relative  performance  of 
explosives  is  essentially  the  same  at  sea  level  as  at  altitude  sc  that 
equivalent  weights  do  not  vary  with  altitude.  The  equivalent  weight 
relating  various  explosive  compositions  in  bare  spherical  charge  form 
to  airblast  performance  is  defined  in  this  report  as  the  energy 
equivalent  weight,  fg,  and  is  given  as  an  explosive  property  in 
Table  3.1.  Factors  relating  cylindrical  to  spherical  charges  and 
cased  cylindrical  to  bare  cylindrical  charges  have  been  shown  experi¬ 
mentally  to  depend  on  the  peak  incident  overpressure  level  in  a  gross 
sense.  Therefore,  before  evaluating  these  factors,  an  estimate  of 
the  peak  incident  overpressure  for  a  spherical  explosion  scaled  to  sea 
level  is  made  based  on  the  energy  equivalent  weight  alone. 

First,  methods  were  developed  to  determine  the  cylindrical  charge 
equivalent  weight.  A  compilation  of  data  for  bare  Comp  B  cylindrical 
charges  was  taken  from  reference  (7)  for  L/D  ratios  between  2  and  10. 

It  was  found  that  the  data  formed  the  three  curves  shown  in  Figure  4.2. 
The  90°  curve  gives  peak  incident  overpressures  measured  along  a  line 
perpendicular  to  the  longitudinal  axis  of  the  cylindrical  charges; 
the  45°  curve  gives  overpressures  along  a  line  inclined  45°  from  the 
longitudinal  axis;  and  the  0°  curve  gives  overpressures  along  the 
extension  of  the  longitudinal  axis.  The  orientation  of  the  projectile 
with  respect  to  the  aircraft  compartment  structure  would  vary  con¬ 
siderably  depending  on  the  mode  of  attack.  Since  the  design  of  air¬ 
craft  to  withstand  internal  blast  is  of  utmost  concern  for  this  program, 
a  conservative  assumption  is  to  use  the  curves  yielding  the  highest 
pressure.  In  Figure  4.2  this  is  the  90°  curve  up  to  a  scaled  distance 
of  about  8  and  then  the  45°  curve  for  scaled  distances  greater  than  8, 
or  the  resultant  composite  cylindrical  charge  curve  shown  in  Figure  4.3. 
(If  one  is  interested  in  weapon  selection  for  damaging  aircraft,  he 
would  chose  the  composite  curve  of  45°  and  0°  in  Figure  4.2.) 
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Prom  the  same  DRI  study  in  reference  (7)>  bare  spherical  charges 

of  the  same  explosive  Comp  B  were  detonated  yielding  the  spherical 

charge  curve  in  Figure  4.3.  Chosing  a  particular  pressure  level,  one 

can  determine  the  cylindrical  charge  equivalent  weight  (f_  *  weight 

s 

of  sphere/weight  of  cylinder)  for  equal  distances.  In  this  manner  the 
low-pressure  curve  shown  in  Figure  4.4(a)  was  developed.  Empirical 
relations  that  represent  this  curve 

0  <  <  20  ;  f  «  1.45  (4.6) 

20  <  APi  ;  fs  «  0.613  (AP± )°* 287  (4.7) 

have  been  programmed  into  the  computer  code.  It  is  assumed  in  the 
program  that  all  explosives  follow  the  behavior  of  this  experimental 
data  for  Comp  B.  Lack  of  complete  sets  of  data  for  other  explosive 
compounds  makes  this  assumption  necessary. 

In  Figures  4.2,  4.3,  and  4.4,  the  experimental  data  do  not  extend 
to  overpressures  above  100  psi.  To  limit  shock  calculations  to  this 
experimentally  verified  range  is  too  restrictive  for  a  realistic 
problem  where  the  projectile  will  generally  be  relatively  close  to  a 
structure  wall  where  incident  overpressures  above  100  psi  will  surely 
exist.  It  is  improper  to  assume  that  the  cylindrical  charge  equivalent 
weight  will  continue  to  increase  indefinitely  with  pressure  as  given 
by  equation  (4.7).  As  the  distance  from  the  cylindrical  charge 
decreases,  at  some  point  the  charge  will  begin  to  appear  as  an 
infinitely-long  charge  or  a  line  charge,  and  the  equivalent  weight 
will  begin  to  decrease.  Since  no  experimental  data  are  available  to 
provide  guidance  for  determining  equivalent  weights  for  pressures 
above  100  psi,  the  following  method  is  assumed.  From  theoretical  work 
on  line  charges  developed  by  Kirkwood  and  Brinkley  in  reference  (8), 
the  high-pressure  curve  shown  in  Figure  4.4(b)  was  determined.  Note 
that  the  ordinate  is  not  equivalent  weight  as  in  Figure  4.4(a),  but 
rather  a  term  defined  for  this  report  as  comparative  weight  index. 

From  the  L/D  ratio  of  the  charge,  this  index  can  be  converted  to  the 
equivalent  weight,  f  .  If  one  makes  this  conversion  over  the  entire 
curve,  a  family  of  curves  is  generated  and  shown  as  the  various  L/D 
curves  in  Figure  4.5.  The  transition  from  a  line  charge  to  a  cylindrical 
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charge  of  finite  length  is  made  by  extending  the  low-pressure  curve 
in  Figure  4.4 (a),  or  equation  (4.7),  until  it  intersects  this  family. 
Thus  Figure  4.5  represents  the  composite  example  of  the  method  used 
by  the  computer  to  calculate  the  cylindrical  charge  equivalent  weight. 
Although  the  trend  of  the  assumed  method  agrees  with  the  expected 
physical  behavior  of  cylindrical  charges,  any  shock  calculations  based 
on  this  method  must  be  viewed  as  approximations.  Because  of  the 
uncertainty  associated  with  this  approach  and  the  lack  of  experimental 
verification  of  pressure  data  above  100  psi  for  cylindrical  charges 
in  general,  a  diagnostic  or  warning  statement  in  the  computer  output 
appears— "CHARGE  SHAPE  CORRECTION  IS  CRUDE.  PSI  EXCEEDS  RANGE  OF 
EXPERIMENTAL  DATA". 

Secondly,  a  method  to  determine  the  effects  on  airblast  shock  of  a 
metal. casing  surrounding  a  cylindrical  charge  was  required,  i.e.,  a 
casing  equivalent  weight,  fc,  to  relate  cased  and  bare  cylindrical 
charges.  As  found  in  reference  (9),  a  number  of  methods  have  been 
proposed  and  are  given  in  Figure  4.6.  These  curves  alone  give  no 
insight  to  the- best  approximation  of  the  casing  effects.  Assorted 
case  effects  experimental  data  taken  from  references  (10)—  (13)  have 
been  plotted  in  Figure  4.6.  Whereas  it  appears  that  the  equation 

f  «  0.20  +  0. 80/(1  +  M/C);  M/C  *  case  weight/charge  weight 

v 

best  fits  the  experimental  data,  it  is  noted  that  it  slightly  under¬ 
estimates  the  effects  of  the  case  for  much  of  the  data.  Consistent 
with  previous  assumptions  for  the  cylindrical  charge  equivalent 
weight,  the  most  conservative  method  was  sought,  i.e.,  the  method 
that  yields  the  greatest  pressures.  This  was  accomplished  by  combining 
the  upper  two  curves  into  one  as  plotted  in  Figure  4.7  with  a  replot 
of  the  experimental  data  points.  It  is  noted  that  only  one  data  point 
lies  above  this  curve.  Therefore,  the  casing  equivalent  weight  used 
in  the  computer  program  is  expressed  as 
0  <  M/C  <  0.53 

fc  *  -  1  -  (M/C)2/(l  +  M/C)  (4.8) 

(M’  =  M/C  for  values  of  M/C  less  than  1) 
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0.53  <  M/C 

f  «  0.47  +  0. 53/(1  +  M/C)  (4.9) 

C 

(Again,  if  one  were  interested  primarily  in  weapon  selection  rather 
than  aircraft  design,  the  lower  curves  of  Figure  4.6  would  be  more 
suitable.)  It  should  be  pointed  out  that  the  experimental  data  used 
for  casing  effects  were  based  on  measured  incident  overpressures  less 
than  100  psi.  Therefore,  a  warning  statement  appears  in  the  computer 
output  to  alert  the  user  to  the  approximate  nature  of  the  data  for 
overpressures  above  100  psi— "CASE  WEIGHT  CORRECTION  IS  CRUDE.  PSI 
EXCEEDS  RANGE  OF  EXPERIMENTAL  DATA”. 

From  the  above  discussions,  a  cylindrical,  cased,  non-TNT  explosive 
charge  can  be  related  to  a  bare  spherical  TNT  charge  by  the  expression  ^ 

WTNT  '  «  X  fe  X  fs  X  fc  (4.10) 

where 

WTNT  *  weight  of  equivalent  bare  spherical  TNT  charge 

W  «  weight  of  cased  charge 

f  »  energy  equivalent  weight  from  Table  3.1 

v 

f  *  cylindrical  charge  equivalent  weight  from  Figure  4.4 
or  4.5 

f  ■  casing  equivalent  weight  from  equations  (4.8) and  (4.9) 

v 

Free-Fleld  Incident  Pressure-Time  and  Impulse.  Using  equation  (4.10) 
and  the  scaling  equation  (4.2),  the  computer can  scale  a  high  explosive 
projectile  explosion  at  altitude  to  a  free-field,  bare,  spherical  1-lb 
TNT  explosion  at  sea  level.  For  a  specific  scaled  distance,  the 
computer  selects  appropriate  free-field  explosion  data  required  for 
the  pressure-time  equation  (4.1).  It  then  calculates  incident  free- 
field  overpressures  that  correspond  to  equal  time  steps  during  the 
positive  phase  duration  of  the  shock  wave.  The  number  of  equal  time 
steps,  k,  can  be  varied  from  10  to  40;  however,  the  number  should  be 
as  large  as  conveniently  possible  because  these  time  steps  control 
the  numerical  integration  procedure  that  calculates  the  positive 
impulse  as  seen  by  the  following  equation 
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pjfere 

I  *  incident  impulse 
AP  ■  incident  overpressure 
i  *=  step  index 
At  *  time  interval 
td  *  positive  duration 
k  *  number  of  time  steps 

For  the  convenience  of  the  user,  time  is  measured  both  from  the 
instant  of  detonation  and  from  the  Instant  of  shock  arrival  at  the 
desired  distance  from  the  explosion.  The  tabulated  incident 
pressure-time  information  and  the  incident  impulse  are  then  scaled 
to  the  actual  conditions  at  altitude  using  the  scaling  equations 
(4.3) — (4.5) . 

Normally  Reflected  Pressure-Time  and  Impulse.  As  stated  in  the 
general  program  description,  accurate  analysis  of  shock  reflection  in 
a  structure  of  an  arbitrary  configuration  is  presently  beyond  the 
scope  of  this  code.  Normally  reflected  pressure-time  information  and 
normally  reflected  impulse  have  been  chosen  as  loading  indices  for 
studying  structural  response  to  shock  loading.  For  the  snail  structural 
compartments  in  aircraft,  shock  loading  from  the  relatively  small 
explosive  charges  in  anti-aircraft  projectiles  would  be  completed 
before  any  appreciable  plastic  response  of  the  structure  has  occurred. 
Therefore,  the  shock  basically  can  be  treated  as  an  impulsive  load 
on  an  aircraft  compartment  structure,  and  it  is  believed  that  the 
normally  reflected  explosion  data  developed  be  the  computer  code 
will  be  sufficient  to  index  the  response  of  the  structure  to  shock 
loads.  While  this  assumption  is  sufficient,  in  all  probability, 
when  applied  to  small  aircraft  compartment  structures,  it  might  prove 
restrictive  and  limiting  for  code  application  to  response  problems 
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relating  to  large  structures  such  as  ship  compartments  or  building 
rooms . 

Methods  for  predicting  the  peak  normally  reflected  overpressure 
have  been  developed  and  verified  by  experimental  data  over  a  wide 
lange  of  pressure.  The  method  used  in  this  computer  program  is 
based  on  the  reflection  factor  curve  developed  by  Brode  in  reference  (1*0 
and  shown  in  Figure  4.8  as  the  solid  curve.  With  sufficient  accuracy 
this  curve  is  approximated  by  the  combination  of  dashed  curves  shown 
in  this  figure.  In  the  computer  program  the  following  equations  are 
used  to  calculate  reflection  factors  for  the  peak  normally  reflected 
overpressure  at  sea  level  conditions. 


0  <  AP^^  <  200  psi 


f  =  2 


(7X14.7)  +  4  AP. 


n  (J  7)  (14.7)  +  APj^ 
200  <  APi  <  10,000  psi 


(4.11) 


fR  =  -3.18  +  3.97  10glQ  (APi) 


(4.12) 


10,000  <  APi 


fR  =  13 


where  AP^  =  peak  incident  shock  overpressure, 


(4.13) 


The  major  problem  now  is  how  to  relate  this  normal  reflection 
factor  derived  for  the  peak  reflected  overpressure  to  the  entire 
reflected  pressure-time  history  and  reflected  impulse.  Whereas 
hydrocodes  exist  that  will  follow  the  normally  reflected  shock 
phenomena  in  time,  they  do  not  lend  themselves  economically  for  use 
with  this  program.  For  the  lack  of  a  better  method  at  this  time,  it 
is  assumed  that  an  adequate  approximation  to  the  reflected  pressure¬ 
time  history  is  found  by  multiplying  the  pressure  level  of  the 
previously  calculated  incident  pressure-time  history  by  the  reflection 
factor,  fR,  derived  for  the  peak  reflected  pressure.  Thus  the 
computer  program  multiplies  the  incident  pressures  as  they  are 
determined  for  the  1-lb  TNT  sphere  at  sea  level  by  the  appropriate 
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reflection  factor,  fR,  from  equations  (*1.11) — (*1.13)  using  the  peak 
incident  peak  overpressure.  It  then  scales  these  results  to  altitude 
conditions  in  the  same  manner  as  it  scales  the  incident  pressures. 
Likewise,  the  incident  impulse  is  multiplied  by  the  same  reflection 
factor  to  obtain  the  normally  reflected  impulse. 

If  one  wishes  to  evaluate  a  reflection  condition  other  than  normal, 
it  is  possible  to  modify  the  incident  pressure-time  curve  with  a 
reflection  factor  for  an  angle  of  incidence  other  than  normal  (90°). 
Such  reflection  factors  can  be  found  in  Figure  *1-6  of  reference  (15) 
and  in  Figure  3.71b  of  reference  (16).  With  the  uncertainties  involved 
in  this  approach  to  reflected  pressure-time  histories,  caution  should 
be  exercised  in  using  the  referenced  factors  for  other  angles  of 
incidence  with  this  method  except  to  serve  as  an  index  or  for  scoping 
calculations. 
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Comparisons  with  Experimental  Data.  Since  normally  reflected  pressure¬ 
time  and  impulse  information  are  assumed  to  be  the  important  shock 
characteristics  in  terms  of  aircraft  compartment  structure  damage, 
it  is  most  important  that  the  code  predicts  normally  reflected  shock 
phenomena  accurately.  Unfortunately,  documentation  of  experimental 
programs  studying  reflected  shock  data  relatively  close  to  the 
explosion  has  been  difficult  to  find.  The  best  available  set  of 
data  was  found  in  the  BRL  study  reported  in  reference  (5).  This 
report  presented  experimental  curves  for  peak  reflected  pressure  and 
reflected  impulse  based  on  old  and  new  tests  with  bare,  spherical 
pentolite  charges. 

These  curves  are  shown  in  Figures  *1.9  and  4.10  along  with 
calculated  results  from  the  computer  program  depicted  by  the  circles. 
Peak  reflected  pressure  predictions  agree  remarkably  well  with  the 
experimental  curve  in  Figure  4.9.  Reflected  impulses  agree  well 
with  the  experimental  curve  in  Figure  4.10,  but  the  relative  variation 
Is  not  as  good  as  for  peak  reflected  pressure.  In  all  fairness  to 
the  computer  program,  a  close  study  of  the  spread  of  experimental 
data  from  which  t^e  impulse  curve  was  drawn  (Figure  8  of  reference  (5)) 
reveals  experimental  variation  as  large  as  that  observed  from  the 
code  prediction-experimental  curve  comparison.  It  is  important  to 
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note  here  that  these  experimental  data  extend  to  the  very  high 
pressure  range  and  that  the  reflected  impulse  for  some  of  the 
tests  include  effects  inside  the  contact  surface — all  questions  of 
uncertainty  in  the  formulation  of  the  code.  Therefore,  it  must  be 
concluded  that  the  computer  code  yields  predictions  that  agree 
remarkably  well  with  this  set  of  experimental  data,  which  certainly 
provides  confidence  to  the  use  of  the  computer  program  for  shock 
calculations. . 

Additional  confidence  is  gained  from  a  comparison  of  a  reflected 
pressure-time  trace  from  one  of  the  BRL  experiments  with  a  code 
predicted  reflected  pressure-time  history.  Figure  6  of  reference  (5) 
gives  an  enlargement  of  a  reflected  pressure-time  trace  from  a  1/8-lb 

l/o 

pentolite  test  at  a  scaled  distance  of  2.5  ft/lb  .  This  curve, 
which  is  free  of  extraneous  noise  and  oscillations,  is  shown  in 
Figure  4.11  as  the  solid  curve.  Shown  as  the  dashed  curve  is  the 
predicted  reflected  pressure-time  results  from  the  computer  code. 
Agreement  has  to  be  classified  as  excellent  in  light  of  all  the 
simplifying  assumptions  used  in  the  code  for  shock  calculations. 
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FIG.  4.4  CYLINDRICAL  CHARGE  EQUIVALENT  WEIGHT 
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CURVES  FROM  REF.  9 
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FIG.  4.6  VARIOUS  METHODS  FOR  PREDICTING  CASE  EFFECTS 


29 


FIG.  4.8  NORMAL  REFLECTION  FACTORS 
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SCALED  DISTANCE,  R/V/,/3  (FT/LB,/3) 


FIG.  4.10  REFLECTED  POSITIVE  IMPULSE-DISTANCE  CURVE  FOR 
PENTOLITE  IN  FREE-AIR  AT  SEA  LEVEL 
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CHAPTER  5 

CONFINED-EXPLOSION  GAS  PRESSURE  CALCULATIONS 


Phenomena  Description.  The  development  of  the  quasi-static  pressure 
that  exists  in  a  closed  structure  after  an  explosion  is  presented  in 
detail  in  references (17)  and  (18).  It  is  briefly  discussed  here. 

After  the  multitude  of  shock  reflections  from  an  explosion  in  a 
completely  closed  structure  have  dissipated,  there  exists  a  significant 
overpressure  in  the  structure.  A  tremendous  amount  of  heat  is  released 
from  the  chemical  decomposition  of  the  explosive  charge  and  from 
subsequent  reactions  with  oxygen  in  the  surrounding  air  in  the  structure. 
Mixing  of  the  extremely  hot  explosion  gas  products  with  the  initial 
gas  in  the  structure  results  in  an  elevated  equilibrium  temperature 
of  the  gas  mixture.  Since  the  volume  of  the  structure  remains  essen¬ 
tially  constant  during  the  explosion,  the  elevated  temperature  must 
be  accompanied  by  an  increase  in  the  equilibrium  pressure  of  the  gas 
mixture.  The  process  can  be  viewed  to  be  similar  to  a  combustion  test 
in  a  bomb  calorimeter.  The  pressure  will  slowly  decay  with  time  due 
to  heat  losses  to  the  structure  walls;  however,  in  comparison  with  the 
highly  transient  nature  of  the  shock  phenomena,  this  pressure  can  be 
truly  defined  as  quasi-static. 

Historically  in  the  literature,  this  quasi-static  pressure  has 
been  known  by  different  names  such  as  static  pressure,  steady  over¬ 
pressure,  internal  blast  pressure,  post-detonation  pressure,  and 
chamber  pressure.  It  is  assumed  that  the  reason  that  no  single  name 
has  evolved  is  because  there  has  been  misinterpretation  of  existing 
names  or  it  has  been  felt  that  existing  terms  do  not  adequately 
describe  the  phenomena.  Therefore  to  add  to  the  growing  list  of 
names  and  hopefully  to  clarify,  this  quasi-static  pressure  created 
by  mixing  the  hot  explosion  gas  products  with  the  initial  gas  in  the 
closed  structure  is  simply  defined  in  this  report  as  the  confined- 
explosion  gas  pressure. 
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Existing  Methcds  of  Calculation.  Currently  there  are  two  commonly  used 
methods  for  estimating  the  magnitude  of  the  confined-explosion  gas 
pressure;  that  proposed  by  Filler  in  references (17)  and  (18)  and  that 
proposed  by  Weibull  in  reference  (19).  Filler  proposed  that  the 
confined-explosion  gas  pressure  can  be  calculated  from  an  expression 
equivalent  to 

APg  «  (UhW)/V0 

where  AP  *=  confined-explosion  gas  pressure  (overpressure),  psi 

o 

h  =  heat  of  combustion  of  explosive,  cal/gm 

W  *  weight  of  explosive,  lb 

■a 

V  *  volume  of  closed  structure,  ftJ 

This  method  assumes  that  there  is  sufficient  oxygen  in  the  initial 
air  in  the  closed  structure  to  ensure  that  an  oxygen-deficient 
explosive  will  achieve  complete  combustion.  It  also  assumes  that 
the  specific  heat  of  the  gas  mixture  remains  constant.  This  approach 
was  verified  for  small  quantities  of  different  explosives  detonated 
in  a  large  air-filled  chamber  resulting  in  modest  confined-explosion 
gas  pressures  up  to  about  30  psi.  Realizing  the  deficiency  in  the 
use  of  the  heat  of  combustion  in  a  possibly  oxygen-poor  atmosphere, 
Filler  conducted  experiments  in  an  inert  atmosphere  and  found  results 
that  indicated  the  heat  of  detonation  yielded  accurate  agreement  for 
this  case,  as  expected.  Unfortunately  these  studies  did  not  determine 
analytical  relations  that  describe  the  phenomena  in  the  transition 
region  between  the  heat  of  combustion  and  heat  of  detonation.  Neither 
did  they  extend  to  the  high-pressure  region  where  the  effects  of 
variations  in  gas  specific  heats  could  be  observed  readily. 

Weibull  proposed  that  the  confined-explosion  gas  pressure  for  a 
TNT  charge  can  be  calculated  from  the  expression 

AP_  •  2410  (W/V_)0,72 

o  ^ 

This  method  was  an  empirical  fit  to  experimentally  measured  pressures 
from  TNT  explosions.  Unlike  Filler's  method,  there  are  no  means  of 
relating  this  equation  to  explosives  other  than  TNT.  However,  Weibull' s 
experimental  data  extends  into  the  high-pressure  range  (near  1000  psi) 
where  obviously  the  specific  heats  0:  the  gas  mixture  components  are 

changing  and  the  transition  between  heat  of  combustion  and  heat  of 
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detonation  can  be  observed.  Unfortunately,  this  study  was  limited 
to  an  empirical  approach  without  fully  exploring  the  underlying 
ohenomena. 


Need  for  Improved  Method.  Figure  5.1  gives  the  prediction  curves 
proposed  by  Filler’s  method  and  Weibull’s  method.  Weibull’s  extensive 
TNT  data  are  also  plotted  for  direct  comparison.  The  deficiencies 
of  these  two  methods  become  obvious  from  the  comparison.  Because 
complete  combustion  and  a  constant  specific  heat  of  the  gas  mixture 
were  assumed,  Filler’s  method  becomes  decreasingly  accurate  as  the 
pressure  level  increases.  Even  if  the  heat  of  detonation  is  used  with 
Filler’s  method  (the  lowest  curve  in  Figure  5.1),  the  deficiencies  in 
handling  the  transition  region  are  easily  recognized.  Weibull’s 
curve  approximates  the  TNT  data  better  than  Filler’s  method  over  the 
range  of  data,  but  it  is  all  too  clear  that  important  physical 
phenomena  are  being  glossed  over  in  the  empirical  treatment  of  the 
problem  that  makes  it  impossible  to  extend  this  method  to  any  explosive 
other  than  TNT. 

Since  the  confined-explosion  gas  pressure  is  believed  to  be  the 
most  important  loading  parameter  in  the  aircraft  internal  blast 
problem,  it  was  imperative  that  an  improved  method  for  calculating  the 
confined-explosion  gas  pressure  be  developed.  The  following  sections 
describe  the  technique  contained  in  the  computer  program  for  predicting 
the  confined-explosion  gas  pressure  and  comparing  code  results  with 
available  experimental  data. 


Description  of  Improved  Method.  The  improved  method  assumes  an 
explosion  in  a  closed  structure  of  volume,  V  ,  filled  with  air  at  some 
ambient  pressure,  P_,  and  temperature,  T  .  The  explosive  is  limited 
to  a  hydrocarbon  form  of  the  elements  C,  H,  0,  and  N  with  aluminum 
being  the  only  possible  metallic  additive.  Since  most  explosive 
compounds  are  oxygen-deficient,  it  is  assumed  that  the  reaction 
can  consume  all  of  the  oxygen  in  the  air  in  the  closed  structure,  if 
needed.  This  basically  is  assuming  optimum  mixing  and  reaction.  The 
code  calculates  the  number  of  moles  of  air  initially  in  the  closed 
structure  volume  from  the  perfect  gas  law.  One  mole  of  air  is  assumed 
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to  be  composed  of  0.21  mole  02  +  0.79  mole  N2.  From  the  C,  H,  0, 

N,  AL  composition  of  the  explosive  charge  given  as  weight  fractions 
in  Table  3.1,  the  code  calculates  the  number  of  moles  of  each  of 
these  elements. 

The  chemical  reaction  of  the  explosion  and  mixing  with  the  air 
in  the  closed  structure  creates  the  combustion  products  HgO,  AI^O^, 

CO,  C02,  02,  and  N2<  A  priority  in  the  reaction  is  assumed  as  follows; 
(1)  the  hydrogen  in  the  explosive  reacts  with  oxygen  such  that  all 
hydrogen  appears  as  H20,  (2)  the  aluminum  has  next  priority  on  the 
oxygen,  such  that  all  the  aluminum  appears  as  the  solid  ALgO^,  (3)  if 
there  is  an  overabundance  of  oxygen  in  the  explosive  and  structure 
air,  complete  combustion  occurs  such  that  all  carbon  appears  as  C02 
and  the  remaining  oxygen  not  needed  in  any  of  the  reactions  appears 
as  02,  (4)  if  there  is  insufficient  oxygen  in  the  system  after  the 
H20  and  AI^O^  reactions,  then  CO  and  C02  are  produced  in  quantities 
given  by  the  following  equations 

n(C )  +  m(0)  -►  a(C0)  +  b(C02> 

a  +  b  =  n  _  a  *  2n  -  m. 

or 

a  +  2b  =  m  b  *=  m  -  n 

where 

a  **  number  of  moles  of  CO  produced 
b  =  number  of  moles  of  C02  produced 
n  *»  number  of  moles  of  C 
m  =  number  of  remaining  moles  of  0 

and  no  02  exists  in  the  final  gas  mixture,  (5)  the  nitrogen  does  not 
participate  in  the  reaction  and  appears  as  N2  in  the  final  gas  mixture. 
From  the  above  calculations  the  number  of  moles  of  component  gases 
(H20,  CO,  C02,  02,  N2)  that  make  up  the  final  gas  mixture  in  the 
closed  structure  are  known. 

The  formation  of  H20,  AL20,,  CO,  and  jC>2  in  this  combustion-type 
process  releases  a  large  amount  of  heat  energy.  Respective  standard 
heats  of  formation  are  multiplied  by  the  moles  of  individual  gas 
components,  and  th°  sum  of  these  quantities  is  defined  for  use  jn  this 
report  as  the  heat  of  reaction.  The  heats  of  formation  for  the  gas 
products  are  negative  by  standard  thermodynamic  terminology,  i.e.,  if 
energy  is  released  to  the  surrounding  atmosphere,  the  heat  of  formation 
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is  negative.  Thus  the  heat  of  reaction  is  likewise  negative.  However 
for  convenience  in  this  report,  it  is  desirable  to  express  the  total 
amount  of  energy,  released  by  the  explosion  as  a  positive  quantity. 
The  heats  of  formation  of  the  gas  products  and  the  heat  of  reaction 
are  treated  as  positive  quantities  in  the  computer  program.  To  account 
for  the  heat  of  formation  of  the  explosive  compound  in  determining  the 
total  energy,  Q,  it  is  necessary  to  add  the  heat  of  formation  of  the 
explosive  compound  given  in  Table  3.1  to  the  heat  of  reaction.  (Signs 
of  values  in  Table  3.1  conform  to  standard  thermodynamic  terminology.) 

Ar  a  computational  model  only,  the  gas  components  of  the  final  gas 
mixture  in  the  closed  structure  are  assumed  to  exist  at  the  inicial 
ambient  pressure,  P  ,  and  temperature,  T  .  of  the  air  in  the  initial 

a  a** 

volume,  V  .  The  energy,  Q,  is  then  added  to  the  gas  mixture,  but  it 
is  added  in  100°F  steps  in  temperature. 

It  is  well  known  that  the  addition  of  heat  to  a  gas  in  a  constant 
volume  system  follows  the  perfect  gas  relation 

AQ  =  r.  Cy  AT  (5.1) 

where  AQ  «  heat  added 

n  =  moles  of  gas 

Cv  *  specific  heat  of  gas  at  constant  volume 
AT  =  change  in  temperature 

One  of  the  weaknesses  of  previous  methods  for  determining  the  confined 
explosion  gas  pressure  was  that  the  variation  in  Cy  with  temperature 
was  neglected.  Given  in  the  literature,  reference  (20),  are  equations 
relating  the  specific  heat  at  constant  pressure,  Cp,  with  temperature 
for  the  various  component  gases  in  the  final  gas  mixture.  With  the 
assumption  that  the  perfect  gas  relation 

R  =  C  -  C  (R.  *  universal  gas  constant) 

0  p  v  o 

can  be  used,  equation  (5.1)  becomes 

AQ  =  n  (C  -  RQ)AT  (5.2) 

and  direct  use  of  the  Cp  equations  in  reference  (20)  can  be  made 
For  convenience  in  calculation,  the  computer  finds  a  weighted  average 
Cp  to  be  used  in  equation  (5.2)  with  the  total  number  of  moles  of 
gas,  n,  in  the  final  mixture. 
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With  the  total  energy  released,  Q,  and  the  total  number  of  moles, 
n,  of  the  gas  mixture  known,  the  computer  uses  the  following  numerical 
procedure  to  determine  the  final  temperature  of  the  gas  mixture.  (The 
Initial  temperature  is  taken  at  T  *  T_  and  the  addition  of  Q  follows  a 

CL 

constant  volume  process.)  (1)  The  weighted  average  Cp  for  the  gas 
mixture  is  determined  for  the  temperature,  T.  (2)  For  a  temperature 
step  of  AT  =  100°F,  the  incremental  amount  of  heat,  AQ,  required  to 
change  the  temperature  by  100°F  is  calculated  from  equation  (5*2). 

(3)  The  temperature  of  the  gas  mixture  after  the  step  is  r  *  T  +  AT. 

(4)  The  incremental  energy,  AQ,  is  subtracted  from  the  total  released 
energy,  Q.  (5)  The  calculational  steps  (1)  through  (4)  continue  until 
all  of  the  total  released  energy,  Q,  is  used,  thus  the  final  tempera¬ 
ture,  T^,  is  calculated. 

With  the  final  temperature,  Tf,  determined,  the  perfect  gas  law 
gives  the  final  pressure  of  the  gas  mixture  in  the  closed  volume,  VQ, 
by  the  relation 

pf  »  n  Royvo  (5.3) 

Conventionally  this  pressure  is  expressed  as  an  overpressure,  so  that 
the  confined-explosion  gas  pressure,  P  ,  is  defined  as 

o 

APg  -  Pf  -  ra  WO 

It  should  be  restated  that  this  method  of  calculating  the  confined- 
explosion  gas  pressure  is  limited  here  to  C-H-N-O.  type  explosives 
with  aluminum  as  the  only  possible  metallic  additive.  Since  most 
common  explosives  used  as  fills  in  conventional  weapons  fall  into  this 
category,  this  limitation  is  not  considered  restrictive  to  the  general 
use  of  this  computer  program.  Also  this  improved  method  should. yield 
conservative  results  because  optimum  mixing  and  the  most  efficient 
chemical  reactions  are  assumed. 

Comparison  with  Experimental  Data.  Attention  is  now  directed  to  the 
adequacy  of  this  improved  method.  In  Figure  5.2  Weibull’s  TNT  data 
from  reference  (19)  are  plotted  as  indicated  by  circles.  The  computer 
code  predictions  are  given  as  the  solid  curve.  We  note  that  the 
agreement  with  the  data  is  excellent,  that  the  change  in  slope  of  the 
predicted  curve  follows  the  general  behavior  of  the  data,  and.  that 
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the  data  falls  either  oh  the  curve  or  slightly  below  it  which  demon¬ 
strates  the  conservat4 sm  of  the  new  method.  Prom  this  comparison 
alone,  it  is  concluded  that  this  technique  is  far  superior  to  the 
existing  methods  of  calculating  the  confined-explosion  gas  pressure 
for  TNT. 

Before  assuming  the  generality  of  this  improved  method,  it  is 
necessary  to  make  comparisons  with  experimental  data  from  different 
explosive  mixtures  and  different  initial  ambient  air  conditions.  The 
next  most  complete  set  of  data  is  found  in  reference  (21)  for  a 
RDX/WAX ,  89.5/IO.5  mixture  detonated  in  air  at  sea  level  conditions. 

A  plot  of  the  data  points  (circles)  from  this  study  and  the  code 
predictions  (solid  curve)  are  given  in  Figure  5.3.  Again  the  excellent 
agreement  and  the  conservatism  of  the  improved  method  predictions  are 
noted.  Reference  (21)  also  gives  data  for  this  same  RDX/WAX  mixture 
for  a  reduced  atmosphere  (P„  =  1  psia).  These  data  and  the  code 

o, 

predictions  are  given  in  Figure  5.4,  and  the  same  excellent  agreement 
and  conservatism  are  demonstrated. 

Other  assorted  data  for  different  explosives  were  found  in 
reference  (21),  and  some  aluminized  explosive  data  were  found  in 
reference  (17) — all  for  sea  level  conditions.  There  was  an  insufficient 
quantity  of  these  data  to  construct  curves,  thus  they  are  tabulated 
in  Table  5.1  along  with  calculated  code  predictions.  The  excellent 
agreement  is  again  noted,  especially  for  the  extremely  high-pressure 
PETN  data  and  the  aluminized  RDX  data.  An  interesting  observation 
can  be  made  with  the  aluminized  data.  As  the  percentage  of  aluminum 
increases,  the  overprediction  of  the  confined-explosion  gas  pressuure 
tends  to  increase.  But  even  for  the  unrealistic  mixture  containing 
50%  aluminum,  there  is  only  a  16%  deviation.  This  Increase  is 
believed  due  to  the  assumed  optimum  mixing  and  most  efficient  reaction 
in  the  code.  Evidently  the  aluminum  is  not  able  to  utilize  the 
oxygen  in  the  surrounding  air  to  the  maximum  extent  assumed  in  the 
code  calculation. 

It  is  concluded  from  these  comparisons  that  the  improved  method 
for  calculating  the  confined-explosion  gas  pressure  is  far  superior 
to  any  other  known  existing  technique.  Even  with  the  use  of  perhaps 
a  not-too-reallstic  combustion  type  model  and  the  liberal  use  cf 
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equilibrium  perfect  gas  relations  and  properties  for  high-pressure 
and  temperature  transient  conditions,  the  improved  method  appears 
to  perform  exceptionally  well.  From  these  comparisons  the  method 
appears  capable  of  handling  mono,  composite,  and  aluminized  explosives 
at  sea  level  ambient  conditions  and  reduced  atmospheric  ambient 
conditions.  Therefore  with  justifiable  confidence,  this  improved 
method  is  used  as  the  basis  for  confined-explosion  gas  pressure 
calculations  in  the  computer  program. 

Although  use  of  this  computer  program  has  been  consistently 
limited  to  C-H-N-0  explosives  with  aluminum  as  the  only  possible 
metallic  additive,  there  is  no  theoretical  reason  why  it  cannot  be 
adjusted  to  perform  well  with  other  metallic  additives  or  non-C-H-N-0 
explosives.  This  limitation  arises  only  because  there  exists  no 
experimental  data  on  confined-explosion  gas  pressure  for  these 
different  explosives  that  will  permit  the  establishment  of  a  set  of 
reaction  priorities  similar  to  those  for  C-H-N-0  explosives. 


CHARGE  WEIGHT/CvNFINED  VOLUME  (LB/FT3) 


FIG.  5.1  COMPARISON  OF  EXISTING  METHODS  FOR 

PREDICTING  CONFINED-EXPLOSION  GAS  PRESSURE 
WITH  EXPERIMENTAL  DATA  FOR  TNT  EXPLOSIONS 
IN  AIR  AT  SEA  LEVEL  CONDITIONS 


OVERPRESSURE  (PSI) 


OVERPRESSURE  (PSI) 
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TABLE  5.1 

MISCELLANEOUS  CONFINED-EXPLOSION 
GAS  PRESSURE  DATA 


CALCULATED 

EXPERIMENTAL 

TYPE  OF 

w/v. 

OVERPRESSURE 

OVERPRESSURE 

DEVIATION 

EXPLOSIVE 

(LB/FT3) 

(PSD 

(PSD 

<*> 

RDX/TNT 

60/40 

0.00221 

22.0 

19.9 

+10 

60/40 

0.00442 

41.0 

383 

+  7 

PETN 

0.182 

711 

725 

-2 

PETN 

0.304 

1089 

1110 

-2 

PETN 

0.405 

1405 

1400 

+  0 

RDX/AL/WAX 

98/  0/2 

0.00171 

15.3 

15.6 

-2 

76/22/2 

0.00171 

22.0 

213 

+  3 

63/35/2 

0.00171 

26.6 

243 

♦  9 

48/50/2 

0.00171 

303 

26.0 

♦16 

DATA  FROM  REFS.  17  AND  21 
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CHAPTER  6 


VENTING  CALCULATIONS 


Venting.  Inherent  In  the  preceding  section  was  the  assumption  of  a 
completely  closed  structure,  i.e.,  no  venting  occurs  before  the 
maximum  value  of  the  confined-explosion  gas  pressure  is  established. 
Therefore,  at  the  onset  of  venting  the  initial  conditions  of  the 
confined-explosion  gas  pressure  are  known  from  previous  calculations; 
Pf  (gas  pressure  in  absolute  units),  Tf  (temperature),  n  (total 
number  of  moles  of  gas),  Cp  (average  specific  heat  of  gas  at  Tf), 
and  VQ  (volume  of  gas).  A  combination  of  n,  V  ,  and  the  molecular 
weights  of  the  gas  mixture  components  yields  the  initial  density, p^, 
of  the  gas  mixture.  Gamma,  y ,  (the  ratio  of  specific  heats),  is 
found  with  the  known  value  of  Cp  from  the  perfect  gas  relation 

v  -  V°v  *  V(Ro  -  V  {6a) 


From  code  input  information,  the  constant  backpressure,  Pfa,  against 
which  venting  occurs  and  the  initial  vent  area,  Aq,  are  given. 

The  relations  governing  the  venting  process  have  been  derived 
in  reference  (22)  for  steady  isentropic  flow  through  a  perfect  nozzle. 
In  this  reference,  y  was  taken  to  be  1.4  which  permitted  the  relations 
to  be  expressed  in  closed  form.  However,  since  y  in  this  computer 
program  is  not  1.4  and  is  not  constant,  the  differential  form  of  these 
governing  equations  are  taken  from  reference  (22).  (There  is  a  typo¬ 
graphical  error  in  equation  (15)  of  reference  (22) — (y  -  1)  in  the 
denominator  should  be  (y  +  1).)  These  governing  equations  are: 


for  sonic  flow 
tPi 

3v-l  '  ‘ 
P1  2Y 


>1/Y/p0><7§T> 


A°  At 

if 

0 


(6.2) 
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for  > 


(*r 


for  subsonic  flow 

Y-l 

P1  “  y  AP1 


1=1  1=1 
Y  .  p  Y 
L  rb  i 


(6.3) 


for  Pb  <  P;L  < - 

»  _  \  ■ 


(f*1) 


Throughout  the  venting  process,  these  isentropic  relations  are 
assumed 

T1  =  To  (VPo)(Y“1)/Y 


P1  3  po  (VPo)1/Y 


(6.4) 


(6.5) 


The  terms  are  defined  as  follows 

AP,  *  (PQ  -  P^)  =  arbitrary  pressure  step  increment 

PQ  *  pressure  at  beginning  of  increment  step 

P^  *  pressure  at  end  of  increment  step 
Pb  =  ambient  backpressure  (constant  throughout  venting) 
g  **  acceleration  due  to  gravity 

Y  *  specific  heat  ratio  given  by  equation  (6.1)  (based 
on  Tq  and  assumed  constant  during  increment  step) 

PQ  =  density  at  beginning  of  increment  step 

P1  =  density  at  end  of  increment  step 

Tq  =  temperature  at  beginning  of  increment  step 
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*  temperature  at  end  of  increment  step 
VQ  =  volume  of  structure  (constant) 

Aq  «  vent  area  (constant) 

At  =  time  increment  (parameter  to  be  determined) 

Even  though  the  relations  are  derived  for  steady  flow,  it  is  assumed 
that  they  are  applicable  to  the  venting  problem  because  the  pressure 
step  used  in  the  numerical  solution  of  these  equations  is  sufficiently 
small  that  gas  mixture  properties  can  be  considered  constant  during  a 
single  incremental  step. 

It  is  assumed  in  the  venting  process  that  the  composition  of  the 
gas  mixture  remains  constant,  i.e.,  no  single  component  of  the  gas 
mixture  is  vented  preferentially.  The  pressure  increment  step  used 
in  the  computer  program  is  defined  as 

AP1  «=  (Pf  -  Pfa)/100  (6.6) 

i.-1.,  there  are  100  pressure  increment  steps.  Even  though  the  computer 
performs  100  steps  in  this  calculation,  only  every  tenth  step  is  printed 
as  output.  The  printout  can  be  increased  to  up  to  every  other  step  if 
desired.  During  the  venting  process  as  the  gas  density  is  decreasing, 
the  computer  keeps  a  running  account  of  the  quantities  of  each  component 
of  the  gas  mixture  remaining  in  the  compartment  structure.  The  need 
for  this  procedure  will  become  apparent  in  subsequent  sections. 

The  following  discussion  is  a  description  of  a  typical  venting 

calculation.  (1)  Starting  values  of  PQ,  pQ,  TQ,  y,  AQ ,  VQ,  and  P^  are 

known.  (2)  From  the  pressure  increment  step,  P^  is  calculated  and 
used  to  determine  if  flow  is  sonic  or  subsonic.  (3)  With  the  proper 
equation  (6.2)  or  (6.3)  chosen,  the  time  increment  At  is  determined, 
and  from  t^  «*  t  +  At  the  absolute  time  from  beginning  of  venting 
associated  with  P1  is  found.  (4)  From  equations  (6.4)  and  (6.5)  the 
temperature  and  density  at  the  end  of  the  increment  step  are  determined. 
(5)  From  the  density  the  number  of  moles  of  gas  mixture  components 
remaining  in  the  compartment  are  found.  (6)  From  the  temperature,  a 
new  average  C  is  calculated  from  which  a  new  y  is  determined. 

(7)  Values  at  the  end  of  this  increment  step,  P^,  T^,  P^,  and  new  y. 
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become  the  beginning  values  P  ,  T  ,  pQ,  and  y  for  the  next  step. 

(8)  The  above  procedures  form  a  loop  that  continues  until  the  99th 
increment  step  is  completed  which  is  the  step  immediately  before 
*  Pb .  The  program  is  stopped  here  because  equation  (6.3)  cannot 
be  solved  for  P^  =  Pb . 

Vent  Area  and  Volume  Changes.  The  orderly  procedure  given  above  can 
be  readily  interrupted  to  accommodate  vent  area  and  volume  changes 
in  accordance  with  input  failure  criteria  controlling  damage  propaga¬ 
tion.  By  constantly  monitoring  the  pressure-time  history  of  the 
venting  process,  the  computer  can  easily  adjust  to  changes  from  input 
of  the  type  presented  in  Figure  3*1.  For  a  given  pressure  or  time, 
an  adjustment  in  vent  area  can  be  made  simply  by  changing  the  value 
Aq  in  equations  (6.2)  and  (6.3).  However,  a  volume  change  requires 
not  only  the  change  of  Vq  in  equations  (6.2)  and  (6.3)  but  also  an 
adjustment  in  the  gas  mixture  pressure  because  the  volume  has  changed. 

Upon  wall  failure  in  the  initial  compartment,  it  is  assumed  that 
the  gas  mixture  in  the  initial  compartment  instantaneously  mixes  and 
comes  to  equilibrium  with  the  air  in  the  newly  available  compartment. 
The  conservation  of  energy  states  for  this  process  that  the  sum  of 
the  internal  energy  of  the  gas  mixture  immediately  prior  to  wall 
failure  and  the  internal  energy  of  the  air  contained  in  the  adjacent 
compartment  is  equal  to  the  internal  energy  of  the  new  gas  mixture 
after  the  mixing  process.  (No  further  chemical  reaction  is  assumed  to 
occur. )  In  equation  form,  this  concept  is  stated  as 


PlVo 


n  U 

(vpT)  TvT) 


PaVa  p2  <Vo  +  V 


(6.7) 


where  P^  =  pressure  of  gas  mixture  immediately  prior  to  wall  failure 
Y-^  =  gamma  of  gas  mixture  immediately  prior  to  wall  failure 
VQ  =  volume  of  gas  mixture  immediately  prior  to  wall  failure 
Pa  =  ambient  pressure  of  air  in  adjacent  compartment 
V  =  volume  of  air  in  adjacent  compartment 

cl 

Y_  e  gamma  of  air  (taken  to  be  1.4) 
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?2  -  pressure  of  new  gas  mixture 
Y2  s  gamma  of  new  gas  mixture 


Unfortunately  there  are  two  unknowns  in  equation  (6.7),  P2  and  y2* 

(y2  is  unknown  because  the  gas  composition  and  temperature  have 
changed. ) 

By  keeping  a  running  account  of  the  amounts  of  the  components 
of  the  gas  mixture  before  wall  failure  and  by  calculating  the  amount 
of  oxygen  and  nitrogen  in  the  air  in  the  new  compartment,  the  computer 
calculates  the  composition  of  the  new  gas  mixture  and  finds  the  total 
number  of  moles  of  the  new  gas  mixture,  n^.  Since  the  new  volume 
(VQ  +  V&)  is  known,  the  new  density,  p2,  is  calculated.  The  perfect 
gas  law 

p2  =  (n2F0T2)/(V0  ♦  Va)  (6.8) 

gives  a  second  relation  but  introduces  the  third  variable,  T2>  From 
the  programmed  Cp  equations  as  a  function  of  temperature  for  the  gas 
components,  the  computer  is  capable  of  generating  a  third  relation 
from  the  known  quantities  of  the  gas  components  in  the  new  mixture 

Y2  =  y2(t2)  (based  on  equation  (6.1P  (6.9) 

The  numerical  iteration  solution  of  equations  (6.7) — (6.9)  gives  the 
values  of  P2,  T2,  and  y2*  Since  p2,  n2,  and  the  new  gas  mixture 
components  are  known,  all  of  these  values  become  the  beginning 
parameters  for  the  next  increment  step  in  the  venting  calculation 
method.  Subsequent  wall  failures  controlled  by  input  damage  criteria 
are  treated  in  this  same  manner. 


Verification.  There  are  no  experimental  data  available  to  verify 
this  entire  venting  process  including  vent  area  and  volume  changes. 
There  are  only  limited  data  applicable  to  the  venting  process  without 
area  and  volume  changes.  These  are  given  in  reference  (22)  from  which 
the  venting  equations  were  taken.  Here  venting  of  the  confined- 
explosion  gas  pressure  from  a  test  facility  at  NOL  was  measured. 
Venting  gases  escaped  the  test  chamber  of  the  facility  through  a 
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"S"  shaped  labyrinth  passageway  out  either  a  partially  open  door 
(small  vent  area)  or  an  open  door  (medium  vent  area).  Agreement 
between  equation  predictions  and  experimental  data  was  very  good 
for  the  small  vent  area  v/here  flow  in  the  passageway  was  probably 
sufficiently  slow  not  to  induce  any  type  of  flow  losses.  Agreement 
for  the  medium  vent  area,  which  was  about  seven  times  greater  than 
the  small  vent  area,  was  only  fair  with  the  equations  underpredicting 
vent  times  by  20  to  302.  It  is  believed  that  the  seven- fold  increase 
in  vent  area  produced  relatively  high  flow  velocities  in  the  passage¬ 
way  from  v/hich  significant  losses  slowed  the  venting  process. 

In  reference  (22)  caution  is  expressed  in  using  these  venting 
relations  for  large  vent  areas.  However,  it  has  been  learned  that 

ft 

limited  unpublished  data  from  the  Naval  Ship  Research  and  Development 
Center  (NSRDC)  on  venting  explosion  gases  through  large  openings  agree 
very  well  with  predictions  from  the  venting  equations.  Therefore, 
with  only  limited  confirmation  of  the  venting  procedure,  this  method 
is  employed  in  the  computer  code  for  predicting  the  pressure-time 
history  of  the  confined  gas  mixture. 

Limitations .  The  venting  section  of  the  computer  code  has  not  been 
verified  experimentally  to  any  significant  degree.  Thus  experimental 
evidence  in  this  area  is  needed  to  assign  a  confidence  level  to  this 
section  equivalent  to  that  of  the  shock  and  confined-explosion  gas 
pressure  sections.  Four  assumptions  are  made  in  this  section  that 
need  additional  study.  First,  heat  losses  to  the  surrounding  walls  are 
neglected  as  a  significant  mechanism  to  reduce  the  gas  mixture  pressure 
Second,  a  constant  backpressure  against  v/hich  venting  must  occur  is 
assumed.  Third,  gas  mixing  and  the  establishment  of  pressure  equilib¬ 
rium  occur  instantaneously  with  compartment  wall  failure.  Fourth, 
no  chemical  reactions  occur  with  the  air  in  the  adjacent  compartments 
after  wall  failures.  In  terms  of  the  small  compartments  in  aircraft 
wings  and  significant  venting  to  the  atmosphere,  none  of  these 
assumptions  are  believed  to  be  restrictive  for  aircraft  applications. 
However,  for  large  explosions  and  large  structures  such  as  ship 

Information  received  from  J.  W.  Sykes  (NSRDC). 
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compartments  or  building  rooms,  they  may  indeed  be  restrictive  and 
may  require  additional  study  and  modification.  The  code  is  constructs 
in  a  manner  such  that  modifications  in  these  areas  can  be  easily 
accommodated. 
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CHAPTER  7 
USER’S  GUIDE 

Computer  Requirements.  The  computer  program  is  written  in  FORTRAN 
for  a  CDC  6400  computer,  and  it  should  work  without  change  on  other 
CDC  machines.  The  program  is  straightforward  and  can  be  adapted 
easily  to  other  computers.  The  major  change  that  may  have  to  be  made 
is  the  spreading  out  onto  individual  cards  of  the  statements  that 
are  now  placed  on  a  single  card  and  separated  by  the  $  sign.  Storage 
inquires  less  than  32,000  core  memory  words,  the  compilation  time 
is  about  15  seconds,  and  the  run  time  for  a  single  case  is  about 
1  second  on  the  CDC  6400. 

Program  Structure.  A  complete  flow  diagram  of  the  computer  program 
is  given  in  Appendix  A.  Detail  descriptions  of  the  input  cards  and 
format  are  given  in  Appendix  B.  A  complete  list  of  the  program 
variables  with  their  definitions  are  given  in  Appendix  C.  The  code 
consists  of  the  main  program  BLAST  and  six  subroutines.  The  functions 
of  these  sections  are  as  follows: 

BLAST:  reads  input  data;  does  venting  calculation;  does 
final  portion  of  the  shock  wave  calculation 

MIX:  supplies  new  conditions  (pressure,  volume,  temperature, 

gamma)  after  the  gases  of  two  compartments  are  mixed. 

HEDATA:  contains  tables  of  properties  of  explosive  components 
and  mixtures. 

GAMMA:  supplies  average  specific  heat  ratio  and  internal 

energy  for  a  gas  of  given  composition  and  temperature. 

GASES:  supplies  initial  conditions  in  the  compartment  immediately 
after  the  explosion  occurs  and  the  confined-explosion 
gas  pressure  is  developed. 
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TNT:  supplies  pressure,  distance,  arrival  time,  and  other  data 
for  a  spherical  1-lb  TNT  free-field  explosion  at  sea 
level. 

ARDC:  gives  standard-atmosphere  pressure  and  temperature  for  a 
desired  atltitude. 

A  complete  listing  of  the  entire  program  is  given  in  Appendix  D. 


Printed  Warnings  in  Output.  During  the  running  of  a  problem  with 
this  computer  program,  printed  diagnostic  or  warning  statements  may 
appear  in  the  output  to  alert  the  user.  These  are: 

(1)  "WFACT  NOT  KNOWN,  1.0  IS  USED.” 

This  means  that  no  energy  equivalent  weight  has  been  supplied 
for  the  shock  wave  calculation  with  the  desired  explosive.  The 
program  assumes  a  value  of  1.0,  thus  results  are  equal  to  those 
of  TNT. 

(2)  "CHARGE  SHAPE  CORRECTION  IS  CRUDE.  PSI  EXCEEDS  RANGE  OF 
EXPERIMENTAL  DATA." 

The  cylindrical  charge  equivalent  weight  depends  on  the  peak 
shock  pressure  level.  Above  100  psi,  no  experimental  data  were 
available  for  correlation  and  theoretical  techniques  were  used. 

The  warning  statement  is  printed  to  inuicate  that  shock  data 
for  the  particular  case  under  study  are  approximations. 

(3)  "CASE  WEIGHT  CORRECTION  IS  CRUDE.  PSI  EXCEEDS  RANGE  OF 
EX'!  FRIMENTAL  DATA." 

'.n<  method  for  calculating  the  casing  equivalent  weight  was 
bdC'-J  on  pressure  data  for  100  psi  and  below.  The  warning  state- 
;  is  printed  to  indicate  that  shock  data  for  the  particular 
-ase  under  study  are  approximations  because  the  peak  overpressure 
exceeds  100  psi. 

(4)  "CAUTION— CONTACT  SURFACE  HAS  ARRIVED.  DATA  ARE  CRUDE 
BEYOND  T(KSEC)  AFTER  SHOCK  ARRIVAL 

This  warning  statement  appears  during  the  shock  calculations 
if  the  contact  surface  reaches  the  desired  distance  being 
investigated.  It  indicates  that,  after  the  indicated  tim  ,  the 
shock  data  are  approximations. 
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Changes  to  the  Program.  There  are  several  items  in  the  computer 
program  that  probably  will  be  frequently  changed  by  the  user 
depending  upon  the  problem  under  consideration.  They  involve  the 
addition  of  new  explosives  to  the  data  table  in  the  program  and 
changes  in  the  amount  of  printout  for  the  shock  and  venting  calculations 
If  a  new  explosive  not  in  the  data  bank  is  frequently  used,  the  user 
may  wish  to  add  it  permanently  to  the  subroutine  HEDATA.  There  is 
room  for  11  new  explosives  in  this  table.  Beginning  with  index 
number  30  (cax’d  HEDA0675),  four  data  cards  using  the  format  for 
the  existing  explosives  can  be  inserted  to  input  the  new  explosive. 

The  amount  of  printout  desired  for  shock  and  venting  calculations 
may  vary  with  the  area  of  interest  for  a  particular  problem.  This 
is  easily  changed  by  varying  KMAX1  for  shock  calculations  or  KMAX2  for 
venting  calculations  on  card  BLAS0560.  The  shock  wave  calculation  is 
done  in  KMAX1  steps  in  time,  equally  spaced  within  the  positive  duration 
of  the  overpressure.  The  built-in  value  of  KMAX1  is  10.  If  more 
printout  is  desired,  change  KMAX1  to  20  or  40.  If  more  than  40  lines 
are  desired,  the  dimensions  of  PSI(40),  Tl(40),  T2(40),  and  PSIREF(40) 
must  be  increased.  Values  of  KMAX1  below  10  are  not  recommended 
because  the  numerical  integration  to  obtain  impulse  is  controlled  by 
the  number  of  these  steps. 

The  venting  calculations  are  performed  in  100  fixed  Integration 
steps.  However,  only  every  KMAX2-th  step  is  printed  as  output.  The 
built-in  value  of  KM AX 2  is  10,  giving  10  lines  of  venting  printout. 

If  more  venting  data  is  desired,  KMAX2  can  be  changed  to  5,  4,  or  2. 


Example  Problems.  To  demonstrate  the  use  of  the  computer  code  with 
its  different  options  and  features,  nine  sample  problems  have  been 
run.  They  are  varieties  of  the  following  base  problem. 

Consider  an  8  ft^  compartment.  A  projectile  has  penetrated 

p 

the  compartment  forming  an  opening  of  0.00545  ft  area.  The 
projectile  contains  0.0294  lb  of  explosive  of  composition  74% 

RDX ,  21%  AL,  and  5$  WAX.  It  has  a  length  to  diameter  ratio 
of  2.7  and  a  case  weight  to  charge  weight  ratio  of  4.24. 

Figure  7.1  shows  the  input  cards  for  the  nine  problems  that  have 
been  solved.  (Appendix  B  gives  the  description  and  format  of  input 
data  cards.)  Examples  1-3  show  the  three  ways  to  specify  the  explosive. 
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Examples  5-7  show  the  three  options  specifying  damage  criteria  for 
wall  failure  in  the  venting  calculations.  Example  8  shows  the 
method  of  a  shock  wave  calculation  with  several  distances  specified. 
Examples  4  and  9  demonstrate  problems  at  ap  altitude  other  than 
sea  level.  Specific  descriptions  of  each  example  are  given  in  the 
following  paragraphs . 

Example  1  involves  only  one  input  card.  The  explosive  number 
is  17 >  indicating  the  number  of  the  desired  explosive  in  the  list  of 
subroutine  KEDATA.  At  the  end  of  the  card,  N0PT=1  indicates  that 
only  a  venting  calculation  is  desired;  NV=0  indicates  that  the  chamber 
volume  and  vent  area  remain  at  their  initial  values  throughout  the 
problem;  and  NR=0  indicates  that  no  distances  are  specified  since 
this  is  only  a  venting  calculation.  The  results  of  this  problem  are 
shown  in  Figure  7.2. 

Example  2  involves  two  data  cards.  The  first  card  is  the  same 
as  in  Example  1  except  that  the  explosive  number  is  0.  This  causes 
the  second  card  to  be  read  in.  This  card  gives  the  energy  equivalent 
weight  *  1.30,  the  heat  of  formation  «  29.36  cal/gm,  and  the  weight 
fractions  of  C,  H,  0,  N,  and  AL.  The  results  are  the  same  as  for 
Example  1  which  are  given  in  Figure  7.2. 

Example  3  again  involves  two  data  cards.  The  first  card  has  an 
explosive  number  of  -1.  This  causes  the  second  card  to  be  read  in. 

This  card  gives  the  energy  equivalent  weight  =  1.30  and  the  weight 
fractions  of  the  desired  components  from  the  explosive  list  in 
HEDATA:  74$  number  27  (RDX),  21$  number  25  (AL),  and  5%  number  26  (V'AX). 
Again,  the  results  are  the  same  as  those  of  Example  1  which  are  given 
in  Figure  7.2. 

Example  4  is  the  same  as  Example  1  except  that  the  compartment  is 
at  altitude  rather  than  at  sea  level.  The  ambient  pressure  is  6.76 
psia  and  the  temperature  is  -24.6°C.  The  results  of  this  calculation 
are  shown  in  Figure  7.3. 

Example  5  returns  to  sea  level  but  the  compartment  volume  is 
allowed  to  change.  NV-1  means  that  one  card  of  volume  and  area 
change  data  is  to  be  read.  This  card  contains  the  following  data;  if 
the  confined-explosion  gas  pressure  in  the  tank  exceeds  30  psia,  the 

o  p 

volume  increases  by  4  ft-3  and  the  vent  area  increases  0.00545  ft  . 
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The  ambient  pressure  and  temperature  of  the  air  in  this  additional 
volume  and  vent  area  are  added  only  if  the  confined-explosion  gas 
pressure  exceeds  the  45-psia  level.  The  results  are  shown  in 
Figure  7. 4. 

Example  6  has  two  cards  of  volume  and  area  change  data.  The 
last  number  on  these  cards  is  2  which  indicates  that  the  changes 
of  volume  and  area  are  to  be  made  at  the  indicated  times  of  0.15  and 
0.60  sec.  No  tests  are  made  on  the  pressure,  so  that  the  changes  are 
made  at  the  indicated  times  regardless  of  the  pressure.  The  results 
are  shown  in  Figure  7.5* 

Example  7  has  three  cards  of  volume  and  area  data.  The  last 

number  on  these  cards  is  3  which  indicates  that  if  the  pressure 

exceeds  the  indicated  value  when  the  indicated  time  is  reached,  then 

the  volume  and  area  change  is  made.  For  example,  if  the  pressure  in 

o  p 

the  tank  exceeds  ^5  psia  at  0.15  sec,  4  ftJ  of  volume  and  0.005^5  ft 
of  vent  area  are  added.  The  results  are  shown  in  Figure  7.6. 

Example  8  is  a  shock  wave  calculation  only,  indicated  by  N0PT=2. 
NV«=0  since  no  venting  parameters  are  involved  in  a  shock  wave  calcu¬ 
lation,  and  NR=3  since  three  distances  are  desired.  The  second  card 
contains  these  three  distances:  0.667,  1.000,  and  1.333  feet  from 
the  center  of  the  charge.  The  results  for  the  single  distance  of  0.66 7 
are  shown  in  Figure  7-7. 

Example  9  is  the  same  as  Example  8  except  that  the  compartment  is 
at  altitude.  The  results  are  shown  in  Figure  7-8. 

Explanation  of  Typical  Output.  A  typical  example  of  the  printout  for 
shock  calculations  is  given  in  Figure  7.7  which  are  the  results  from 
Example  8.  The  index  number  and  properties  of  the  explosive  used  in 
the  calculation  appear  at  the  beginning  of  the  output.  The  two  warning 
statements  concerning  the  cylindrical  charge  equivalent  weight  and 
casing  equivalent  weight  are  noted.  Under  "SHOCK  WAVE  CALCULATION", 
the  left-hand  column  repeats  all  the  input  parameters  governing  the 
shock  problem.  In  the  right-hand  column,  certain  constants  derived 
by  the  computer  for  the  calculation  are  given: 

ADJUSTED  WT ( LB  TNT) — equivalent  TNT  sphere  from  equation  (^.10) 

HE  ENERGY  FACTOR — energy  equivalent  weight,  fg,  from  Table  3.1 
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CHARGE  WEIGHT  FACTOR  cylindrical  charge  equivalent  weight, 

f_,  from  Figure  4.5 

CASE  WEIGHT  FACTOR  casing  equivalent  weight,  from  equation 

(4.8)  or  (4.9) 

PRESSURE  SCALE  FACTOR  (Pg/Pa)  for  equations  (4.2)— (4.5) 

DISTANCE  SCALE  FACTOR  (Wg/W&)1/3  (Pa/Ps)1/3  for  equations 

(4.2)  — (4.5) 

TIME  SCALE  FACTOR  (Ws/Wa)1/3  (Pa/Pg)1/3  (Ta/Ts)1/2  for 

equations  (4.2)  — (4.5) 

NORMAL  REFL  FACTOR  normal  reflection  factor,  f^,  from 

equations  (4.11)— (4.13) 

The  tabulated  pressure-time  shock  data  is  noted  for  the  desired 
distance  of  0.667  ft.  Both  the  incident  and  normally  reflected 
overpressures  are  given  as  functions  of  time  where  time  is  measured 
from  the  instant  of  detonation  and  shock  arrival.  For  example,  the 
shock  arrives  at  the  distance  0.667  ft  in  0.07118  msec;  the  peak 
Incident  overpressure  is  317*6  psi  and  the  reflected  overpressure  is 
2144  psi;  and  the  positive  phase  of  the  shock  is  completed  0.192  msec 
after  detonation  or  0.1209  msec  after  shock  arrival.  Next  the  impulses 
for  the  Incident  and  reflected  waves  are  given.  Lastly,  the  warning 
statement  concerning  the  contact  surface  appears  which  states  that 
0.02787  msec  after  tne  snock  arrives,  the  pressure-time  data  are 
approximations . 

Figure  7.6  gives  printout  results  for  Example  7  on  the  confined- 
explosion  gas  pressure  venting  and  subsequent  changes  due  to  structural 
failures.  At  the  beginning  of  the  output  are  the  index  number  and 
properties  of  the  explosive  used  in  the  calculation.  Under  "VENTING 
CALCULATION"  a  repeat  of  input  parameters  is  given;  under  "BEGIN 
VENTING  CALCULATION"  the  input  failure  criteria  table  is  repeated. 

Under  "PROPERTIES  OF  GASES"  the  output  describes  the  condition  of 
the  confined-explosion  gas  in  the  initial  compartment  volume  before 
any  venting  has  occurred.  The  first  statement  indicates  that 
oxidation  was  complete,  i.e.,  sufficient  oxygen  was  available  to 
make  H2O,  AL2O2,  and  only  COg*  Had  there  been  insufficient  oxygen 
for  complete  oxidation,  the  output  would  have  indicated  the  name 
and  quantity  of  the  last  product  formed.  For  example,  if  all  H+HgO 
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and  AL  -*■  ALgO^  but  there  was  insufficient  oxygen  to  completely  react 
with  all  the  C  to  form  CO,  the  computer  would  print  "PERCENT  LAST 
PRODUCT  (CO)  =  (fraction  of  carbon  used)".  Next,  the  computer  prints 
the  maximum  temperature  of  the  confined-explosion  gas,  the  energy 
released  in  the  chemical  reaction  that  creates  the  confined-explosion 
gas  pressure,  the  specific  heat  ratio,  y,  and  the  maximum  value  of  the 
confined-explosion  gas  pressure  expressed  as  an  overpressure. 

Under  "BEGIN  VENTING  OF  GASES"  the  gas  pressure-time  data  are 
tabulated  along  with  the  amount  of  gas  in  the  confining  volume  (GASES), 
the  temperature  of  the  gas  (TEMP),  the  specific  heat  ratio  (GAMMA), 
and  an  index  (NEQN).  If  this  index  is  1  then  the  flow  velocity  is 
sonic;  if  2,  flow  velocity  is  subsonic.  The  beginning  time  is  zero 
for  this  calculation  which  is  set  arbitrarily  after  the  dissipation 
of  the  shock  wave,  and  the  overpressure  is  maximum  at  45.9  psi. 
Adjustments  made  with  compartment  failures  and  continued  venting  are 
noted.  For  example,  at  t=0.15  sec  the  gas  overpressure  is  36.5  psi 
which  is  above  45  psia  and  the  wall  fails.  A  new  pressure  of  36.75  psia 
or  22  psi  overpressure  is  calculated  for  the  new  volume  of  12  ft  ,  and 
venting  continues  through  the  new  area  of  0.0109  ft  until  t=0.6  sec 
when  another  failure  occurs.  The  code  readjusts  the  pressure  to 
accommodate  the  new  volume  and  venting  continues  until  the  overpressure 
is  essentially  zero  at  t=0.9  sec. 
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EXAMPLE 

1 

.0294 

17  2.7 

4.24 

8.  345-3 

14.7 

20. 

EXAMPLE 

2 

.0294 

0  2.7 

4.24 

8.  545-3 

14.7 

20. 

1.30 

29.36 

.163 

•  027 

•  280 

.320 

EXAMPLE 

3 

•0294 

-1  2.7 

4.24 

8.  545-3 

14.7 

20. 

1.30 

27  .74  25  .21 

26  .05 

EXAMPLE 

4 

.0294 

17  2.7 

4.24 

8.  545-3 

6.76-24.6 

EXAMPLE 

5 

•0294 

17  2.7 

4.24 

8.  545-3 

14.7 

20. 

30. 

0. 

4. 

•549-2 

14.7 

20. 

EXAMPLE 

6 

•0294 

17  2.7 

4.24 

8.  545-3 

14.7 

20. 

0. 

.19 

4. 

.545-2 

14.7 

20. 

0. 

.60 

4. 

.545-2 

14.7 

20. 

EXAMPLE 

7 

.0294 

17  2.7 

4.24 

8.  545-3 

14.7 

20. 

43. 

.15 

4. 

.545-2 

14.7 

20. 

20. 

4. 

•545-2 

14.7 

20. 

19. 

.8 

4. 

0. 

14.7 

20. 

EXAMPLE 

8 

•  0294 

17  2.7 

4.24 

8.  545-3 

14.7 

20. 

•447 

1.0 

1.333 

EXAMPLE 

9 

.0294 

17  2.7 

4.24 

8.  545-3 

6.76-24.6 

.647 

1.0 

1.333 

1  0  0 

1  0  0 

10  0 

1  0  0 

1  1  0 

1  2  0 

1  3  0 

2  0) 

2  0) 


FIG.  7.1  INFUT  CARDS  FOR  NINE  EXAMPLE  PROBLEMS 
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INTERNAL  BLAST  04MAGE  MECHANISMS  PROGRAM*  MAR  1972 
RPX/AL/wAX*  74/21/5 

EXPLOSIVE  PROPERTIES 

NUMBER  EQWT  EFORM  EXPLOSIVE  COMPOSITION  BY  HEIGHT 
KCAL/G  C  A  N  0  al 
IT  1,300  .029360  ,163  .027  .200  ,320  .210 

VENTING  CALCULATION 

CHARGE  WEIGHT (LR)  >  .2940E-01 

INlt  VOLUME (CU  FT)  ■  8,000 

INIT  VENT  AREA (SO  FT)  ■  .5450E-Q2 

AMBIENT  PRESSURE (PSlA)a  14.70 

AMBIENT  TEMP (C)  ■  20.00 

CHAMBER  PRESSURE (PSlA)a  14.70 

CHAMBER  TEMP (C)  ■  20.00 

NOPTa  1  NVa  0 

BEGIN  VENTING  CALCULATION 

PROPERTIES  OF  GASES— 

OXIDATION  COMPLETE 
TEMPERATURE*  0E6REES  F  a  1653*2 
ENERGY  RELEASE (KCAL/O)  a  3.4573 
SPECIFIC  HEAT  RATIO  a  1.3141 
«AS  OVERPRESSURE (PSI)  a  45.945 


BEGIN  VENTING  OF  GASES 


OVERPR(PSI) 

TIME(SEC) 

OASES (LB) 

TEMP(R) 

GAMMA  NEON 

45.94 

0. 

.6167 

2113. 

1.3141 

41.35 

.6936E-01 

,5809 

2074, 

1.3155 

1 

36.76 

.1453 

.5443 

2032. 

1.3167 

1 

32.16 

.2293 

.5070 

1986, 

1.3181 

1 

27.57 

,3229 

.4689 

1937. 

1.3197 

1 

22.97 

.4285 

.4297 

1884. 

1.3214 

1 

18.38 

.5496 

.3895 

1825. 

1.3235 

1 

13.78 

.6911 

,3479 

1759. 

1.3259 

1 

12. AT 

.7360 

.3358 

1739. 

1.3267 

1 

7,880 

.9175 

.2921 

1861. 

1.329A 

2 

3.286 

1,171 

,2462 

1570. 

1.333ft 

2 

.6943E-01 

1.690 

.2125 

1494. 

1.3380 

2 

FIG.  7.2  OUTPUT  RESULTS  FOR  EXAMPLES  1,  2,  AND  3 
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INTERNAL  8LAST  OAMAOE  MECHANISMS  MOOR AN t  MAR  107* 
ROX/AL/uAX.  74/21/8 

(XPLOSIVC  PROPERTIES 

NUNRER  COOT  CROON  explosive  composition  it  weioht 


KCAL/O 

C 

M  N 

IT  1.300  .0*0340 

•  143 

.0*7  «*R0 

VENT1N0  CALCULATION 

CNApflC  WCIOHT (Lfl) 

• 

•I040C-01 

INIT  VOLUME (CU  PT) 

9 

0.000 

INIT  VENT  ARM  ISO  PT) 

9 

•SAS0E«0* 

ARRICNT  PRESSURE(PSIA)* 

4.740 

ARRICNT  TEMP(C) 

9  ' 

■*4.40 

CRAMMER  PRESSURE IPSI A) ■ 

4.T40 

CNAMRER  TEMP (C) 

NOPf*  1  NVa  0 

•  < 

■*4.40 

•COIN  VCNTIN6  CALCULATION 


PROPCRTICS  op  oases— 

OXlOATION  COMPLETE 
TEMPERATURE.  OCOPCCS  P  a  *451.7 
CNCROV  RELEASE (KCAL/O)  •  3.4873 

SPCCIPIC  HEAT  PATIO  •  1 4*003 

•AS  OVERPRESSURE (PSD  «  A*. 014 


•COIN  VCNTINQ 
OVERRRIPSI) 

OP  OASES 

timciseci 

OASES (18 > 

TEMP(R) 

•AMMA  neon 

41.01 

0. 

•  34*4 

3114. 

1.2003 

30.4* 

.4T30C-01 

•  3104 

3083. 

1.2005 

34.33 

•  14*0 

•1080 

*008. 

1.2018 

30.04 

.**40 

.*714 

*01*. 

1.20*4 

PScTS 

.3*18 

•*400 

*031. 

1.X030 

*1.44 

.43*3 

•C*l* 

*741. 

1.2053 

jT.1T 

.5438 

.1040 

*440. 

1.2071 

1*.0T 

.7*41 

•  1473 

*8*3. 

1.2093 

0.5*3 

.0*04 

•  1304 

*303. 

1.302* 

8.447 

1.111 

•  1174 

2*40. 

1.3049 

1.354 

1.5*4 

•I400C-01 

*08*. 

1.3104 

.40A4C-Q1 

1.040 

.74S1E-01 

1070, 

1.3137 

2 

FIG.  7.3  OUTPUT  RESULTS  FOR  EXAMPLE  4 
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INTERNAL  BLAST  DAMAGE  MECHANISMS  PROGRAM*  MAR  1972 
R0X/AL/wAX«  74/21/5 

EXPLOSIVE  PROPERTIES 

NUMBER  EQWT  EFORm  EXPLOSIVE  COMPOSITION  BY  WEIGHT 
KCAL/G  C  H  N  0  al 
IT  1.300  .029360  .163  .027  .200  .320  *210 

VENTING  CALCULATION 

charge  weight <lb>  ■  .2940E-01 

INIT  VOLUmE(CU  FT)  ■  B.000 

INIT  VENT  ARFAJSO  FT)  ■  .54S0E-02 

AMBIENT  PRESSUOE(PSlA)«  14,70 

AMBIENT  TEMP (C)  ■  20.00 

CHAMBER  PRESSURE (RSIA)p  14.70 

CHAMBER  TEMP(C)  «  20.00 

NOPt«  1  NV»  1 

BEGIN  VENTING  CALCULATION 

TABLE  OF  VOLUME  AND  vent  AREA  CHANGES 

P(PSIA)  T (SEC)  V(CU  PT)  A (SO  FT)  PAMR(PSIA)  TAMB(C)  NoPTV 

30.00  0.  4,000  .5450E-02  14.70  20.00  1 


PROPERTIES  Op  OASES— 

OXIDATION  COMPLETE 
TEMPERATURE.  DEGREES  F  ■  1653.2 

ENERGY  RELEASE (KCAL/G)  •  3.6573 

SPECIFIC  HEAT  RATIO  •  1.3141 

GAS  OVERPRESSURE (PSD  •  45.945 

FAILURE  LEVEL  IN  TABLE  EXCEEDED. 
VOLUME  INCREASE (CU  FT)«  4,000 
NEW  TOT  VOL  (CU  FT)  m  12,00 

NEW  TOT  AREA  (SO  FT)  *  .1040E-01 

NEW  PRESSURE (P5IA)  ■  43,06 

NEW  GAMMA  ■  1,338 


begin  venting  of  gases 


OVERPR(PSI) 

TIME(SEC) 

GASES (LB) 

TEMP (R) 

GAMMA  NEON 

28.36 

0. 

.9162 

1500. 

1,3377 

25.52 

.5174E-01 

.8708 

1475. 

1.3409 

1 

P2.69 

.1077 

.8246 

1447, 

1.3424 

l 

19.85 

.16-6 

,7776 

1419. 

1.3439 

1 

IT. 01 

,2353 

.7296 

1388. 

1.3456 

1 

14,10 

,3090 

•  6806 

1355. 

1.3475 

1 

12.67 

,3515 

.6541 

1336. 

1,3488 

1 

9**3*. 

.4394 

.6032 

1299. 

1.3509 

2 

7.00« 

.5419 

.5508 

1250. 

1.3535 

2 

4,164 

,6696 

.4967 

1212. 

1,3565 

2 

1.329 

.0619 

.4406 

1161. 

1,3601 

2 

.1943 

1.037 

.4175 

1139. 

1.3622 

2 

FIG.  7.4  OUTPUT  RESULTS  FOR  EXAMPLE  5 
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INTERNAL  OLAST  DAMAGE  MECHANISMS  noMMi  NU  1*7* 

POX/AL/wAX,  74/21/5 

EXPLOSIVE  Memrltl 

NUMRER  rowr  EFORm  EXPLOSIVE  COMPOSITION  *V  WEIOHT 
KCAL/O  C  M  N  0  AL 
17  1.30*  ,M«MO  .U3  ,037  .1*0  .3*0  .*1* 

VENTINO  CALCULATION 

CNAAOE  WEIOMT (LM)  •'  .244SE«01 

INlt  VOLUMEieU  FT)  ■  0.000 

1NIT  VENT  AAEAISQ  FT)  ■  .545QE-0? 

OMOJENT  PRESSURE (PS1 A) ■  14.T0 

AWOIENT  TCMP(C)  ■  *0.00 

CNAMAEO  PRESSURE (PSl A) •  14.70 

CNAMOEP  TCMF<C)  ■  *0.00 

NOPT*  I  NV.  * 

■COIN  VENTINO  CALCULATION 

TOOLE  OF  VOLUME  AND  VENT  AREA  CMANOES 

P(PSIA)  T (SEC)  V(CU  FT)  A(S0  FT)  OAMOlOStA)  TOMO(C)  NOOTV 


0. 

0. 


.1500 

.0000 


PROPERTIES  OF  OASES— 
OXIDATION  COMPLETE 
TEMPEPATUPE.  OEOPEES  F  ■ 
ENEPOY  RELEASE IKCAL/O)  ■ 
SPECIFIC  HEAT  PATIO  ■ 
*05  OVERPRESSURE  (PSD  • 

•COIN  VENTINO  OF  OASES 
OVEPPPIPSI)  TIMEISECI 
45.04  0. 

01,35  ,A434E«01 

30.70  .1453 

30.4*  ,1500 

TINE  MAS  PCACHCO  TV(  !)■ 


FAILURE  LEVEL  IN  TASLE  EXCEEDED. 


4.000 

4.000 


1053.2 

3.0573 

1.3141 

45.445 


•ASCSILO) 

.0107 

.500* 

.5443 

.54*2 

.1500 


.54SSE-02 

,54S0E>0* 


14.70 

14.70 


*0.00 

to. SO 


TEMPI*) 

2113. 

2074. 

203*. 

202*. 


•ANNA  NEON 
1.3141 
1.3155  1 

1.3107  1 

1.310*  1 


VOLUME  INCREASE (CU  FT). 
NEW  TOT  VOL  (CU  FT)  • 
NEW  TOT  AREA  ISO  FT)  • 
NEW  PRESSURE (FSIA)  • 
NEW  OAHMA 


Ft.os 

14.05 

17.04 

15.44 
13.23 
12.04 

10.44 
P.2S3 
0,4*0 


.1500 

.1*04 

.*504 

.3000 

.3074 

.3031 

.451* 

,5*73 

.4000 


4.000 
12.00 
.I040E-01 
34.75 
1.343 
.0417 
,0034 
,7054 
.7204 
.0004 
.0700 
.0341 
.5445 
.SSST 


1*00. 

1370. 

1355, 

1330. 

1304, 

1240. 

1270. 

12*0. 

1*12. 


1.3430 

1,3404 

I.S47R 

1,344* 

1.3500 

1.3513 

1,35311 

1.3550 

1,3570 


TINE  mas  PEACHED  TV(  *)■ 

,4000 

failure  LEVEL 

IN  TAOLC  EXCEEDED. 

VOLUME  INCREASE (CU  FT)w 

4.000 

NEW  TOT  VOL  (CU  FT)  ■ 

10.00 

NEW  TOT  AREA 

(SO  FT)  ■ 

.1435E.S1 

NEW  PRESSURE (PSIA)  ■ 

17.04 

NEW  OAMmA 

• 

1.300 

F.3SR 

.0000 

.050* 

040.7 

1.304* 

2.12* 

,0170 

,0444 

045,4 

1,3000 

1  ,004 

.4351 

,0410 

0*2,1 

1.300) 

J.PS1 

.45*3 

.03*4 

030,7 

1,3005 

l.*15 

.4750 

.0*37 

035,3 

1,3 AO A 

1.1T* 

,0470 

.0151 

031,4 

1.3041 

.4*31 

,7**0 

,0043 

0*0,4 

I, 3044 

.707* 

,7511 

,7474 

0*4,4 

1,3047 

.471 A 

.7050 

,7000 

021,3 

1.340A 

.2354 

.0*40 

.7744 

017,7 

1.3403 

.23SRE.01 

.400* 

.7714 

0l*,4 

1,3404 

FIG.  7.5 

OUTPUT  RESULTS  FOR 

EXAMPLE  6 
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INTERNAl  BLAST  OAMAOP  M^CHtNI^MS  PROGRAM,  MAR  1972 
R0*/*L/wAX»  74/PI/5 

explosive  properties 

NU“»E»  FOwT  EPOHm  rxPLOSIVr  COMPOSITION  BY  WEIGHT 
KCAL/0  C  M  N  o  AL 
IT  1,3##  .029300  .103  ,02?  ,2R0  ,320  .210 

VERTING  CALCULATION 

Charge  weight (lb)  ■  ,2raoe»oi 

INI T  VOLUME  (fl)  FT)  •  A, 000 

IN1T  VENT  ARF A (So  ft)  »  .5450E-02 

AHOIENT  PRESSURE (PSIA).  14,70 

AMRTENT  TEMP)C)  •  10,00 

CHAMFER  PRESSURE (PSlA)a  14, to 

chamfer  TEMPJC)  ■  110,00 

NOPTm  I  NV«  3 

OEOIN  VPNTING  CALCULATION 

TARLE  OP  VOLUME  4N0  VENT  AREA  CHANGES 

P(PStA)  T (SEC)  VfCl)  FT)  A(SG  FT)  P4M«|PSIA)  TAHB(C)  NoPTV 

AS. 00  .1300  4.000  .S450E-02  14.70  20.00 

20.00  .0000  a. 000  .S450E-02  14. TO  20.00 

10.00  .0000  4,000  0,  14. TO  20.00 


PROPERTIES  OP  OASES— 

OXMATION  COMPLETE 

temperature,  degrees  r  ■  1053.2 

ENEROr  RELEASE (KCAL/0)  p  3,0573 
SREciric  heat  ratio  •  1.5101 
GAS  OVERPRESSURE  (PSD  •  43.045 


•eotn  venting  or  oases 


0VE42RIPSI) 

TI«CISEC) 

OASES (LB) 

TEHPIR1 

OAMMA  NEON 

43.04 

0. 

■  4107 

2113. 

1.314) 

41.33 

.OOlOEaOl 

.5000 

2074. 

1.3.53 

1 

30.70 

.105] 

.3443 

2032. 

1.3107 

1 

30.40 

.1500 

.4422 

2020. 

1.3100 

1 

TIME  mas  REACHED  TVf  1 ) a 

•  1SU0 

FAILURE  LEVEL  In  TABLE  EXCEEDED, 

VOLUME  INCREASE (CU  FT) a 

4,000 

NEW  TOT  VOL 

(CU  PT)  a 

12.00 

NEP  TOT  AREA 

(SO  rT|  a 

.lOOOCaOl 

NEW  PRESSURE f PSt A)  a 

30.75 

NEW  oammA 

a 

1.343 

22.05 

.1500 

.0417 

1400. 

1.3430 

1 

10.04 

.1004 

.0030 

1370. 

1.3404 

1 

1?. OA 

.2304 

.7034 

1353. 

1.3470 

1 

15,44 

.3000 

.7244 

1330. 

1,3492 

1 

13.23 

.3674 

.4004 

1304. 

1.3300 

1 

12.00 

.3031 

.4700 

1200. 

1.3313 

1 

10.40 

.4512 

.0301 

1270. 

1.3530 

2 

4.203 

.3273 

,5045 

1240. 

1.3550 

2 

0.420 

,4000 

.4307 

1212. 

1,3570 

2 

TIME  HAS  REACHED  TV!  2) a 

.0000 

PAILURE  LEVEL  IN  TABLE  ElCEEOEO. 

VOLUME  INCREASE ICU  FT).  4,000 
NE«  TOT  VOL  (CU  FT)  ■  10,00 

NO  TOT  AREA  ISO  FT)  ■  .  1435E*01 

NEW  PRFSSUREIPSIA)  a  17,0# 

NEa  gamma  a  1,300 

2.254  .0000  .0302  040.1  1.3004 

2.122  .0170  .0400  045,4  (.300ft 

1.080  .0331  ,0410  042.1  1,3043 

l.OS)  .0543  .0324  830.7  1,3005 

l.»14  .0750  .4237  035,3  1,3444 

1,170  .0074  ,R|51  031,0  1,300) 

.0431  ,7224  ,41)03  020,4  1.3004 

.7074  .7511  ,7074  024,0  1,1007 

.4710  .7050  ,700*  #21.3  t,30»fl 

.2344  ,4304  .7700  017,7  t.70|)7 

.23S4E-01  .0044  .7713  014.4  1.3000 


FIG.  7.6  OUTPUT  RESULTS  FOR  EXAMPLE  7 
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INTERNA).  BLAST  DAMAGE  MECHANISMS  PROGRAM*  MAR  1972 
EXPLOSIVE  PROPERTIES 

NUMBER  EQWT  EFOHm  EXPLOSIVE  COMPOSITION  BY  WEIGHT 
KCAL/G  C  H  N  0  AL 
17  1.300  .029360  .163  *027  .2fl0  .320  .210 

CHARGE  SHAPE  CORRECTION  IS  CRUOE.  PSI  ExCEEOS 
RANGE  OP  EXPERIMENTAL  DATA. 

CASE  WEIGHT  CORRECTION  IS  CRUOE.  PSI  EXCEEOS 
RANGE  OF  EXPERIMENTAL  OATA. 

SHOCK  WAVE  CALCULATION 


input  parameters 

CHARGE  WEIGHT  ADJUSTMENTS 

CHARGE  WEIGHT (LB)  * 

.2940E-01 

ADJUSTED  WT(LB  TNT) 

•6318E 

EXPLOSIVE  NUMBER  * 

17 

HE  ENERGY  FACTOR 

1.300 

L/D  RATIO 

• 

2.700 

CHARGE  SHAPE  FACTOR 

2.89* 

CASE/CHaRGE  hT  ratio  * 

4.2*0 

CASE  WEIGHT  FACTOR 

.5711 

CHAMBER  PRESSURE (PSIA)* 

14.70 

PRESSURE  SCALE  FACTOR 

.9997 

CHAMBER  TEMP(C)  « 

20.00 

OISTANCE  SCALE  FACTOR 

2.511 

ALTITUDE  (KFT)  * 

0. 

TIME  SCALE  FACTOR 

2.489 

NORMAL  REFL  FACTOR 

6.752 

OESIREO  DISTANCE (FT)  * 

.6670 

(CM)  « 

20.33 

time  after 

TIME  AFTER 

INCIDENT 

NORM  REFL 

EXPLOSION 

SHOCK  ARR 

OVERPRESS 

OVERPRESS 

(MSEC) 

(MSEC) 

(PSD 

(PSI) 

7.1180E-02 

0, 

317.6 

21*4. 

9.5353E-02 

2.4I73E-02 

100.2 

676.2 

.107* 

3.6259E-02 

63.11 

*26.1 

.1195 

*.83*S£-02 

♦0.98 

276.7 

.1316 

6.0*3l£-02 

26.81 

181.0 

.1*37 

7.251BE-02 

17.27 

116.6 

.1558 

8.*60*E-Q2 

10.6* 

71.8* 

.1079 

9.6690E-02 

5.91* 

39.93 

.1000 

.1088 

2.49* 

16.6* 

.1920 

.1209 

0. 

0. 

IMPULSE  (PSI. MSEC)  — 

INCIDENT  *  7,675 

REFLECTEC*  51.82 

CAUTION— CONTACT  SURFACE  HAS  ARRIVED. 

DATA  ARE  CRUOE  BEYOND  T(MSEC)  AFTER  SHOCK  ARRIVAL* 


FIG.  7.7  OUTPUT  RESULTS  FOR  FXAMPIE  8 
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INTERNAL  BLAST  DAMAGE  PEOrtrtl.lS.iS  PROGRAM.  MAR  1972 
EXPLOSIVE  PROPERTIES 

NUMBER  ECwT  EFOHM  EXPLOSIVE  COMPOSITION  bt  WEIGHT 
KCAL/G  t  H  N  0  AL 
17  1.300  ,02Y36u  .163  .02?  ,2«0  .320  .210 

CHARGE  SHAPE  CORRECTION  IS  CHUUE.  RSI  tXCEtOS 
RANGE  OP  EXPERIMENTAL  OaTm. 

CASE  WEIGHT  CORRECTION  IS  CKUDt.  PSI  EXCEEDS 
RANGE  OF  EXPERIMENTAL  OaTa. 

SHOCK  WAVE  CALCULATION 

input  parameters  charge  weight  adjustments 


CHARGE  HEIGHTU-B)  * 

«294uE»0l 

ADJUSTED  wT (Lb  TNT) 

.7406E 

EXPLOSIVE  NUMBER  * 

17 

he  energy  factor 

1  .'UO 

L/D  RATIO 

s 

2.700 

charge  shape  factor 

3.394 

CASE/CHaRGE  wt  RATIO  * 

4.240 

CASE  WEIGHT  FACTOR 

.5711 

CHAMBER  PRESSURE (PSIA)* 

6.760 

pressure  scale  factop 

2.174 

CHAHBER  TEMP (C)  • 

•24,60 

DISTANCE  SCALE  FACTOR 

1.638 

ALTITUDE  (KFT)  * 

0. 

TIME  SCALE'FACTOR 

1.707 

NORMAL* REEL  FACTOR 

7.652 

DESIRED  DISTANCE <FT>  « 

.6670 

(CM)  * 

20.33 

TIME  AFTER 

TIME  AFTER 

INCIDENT 

nohh  REFL 

EXPLOSION 

SHOCK  ARM 

OVtRPPt.SS 

OVERPRtSS 

(MSEC) 

(MSEC) 

(PSD 

(PSD 

S.8490E-0? 

U. 

276.5 

21-: 

8,50636-uZ 

2.6573E-02 

o7.1V 

66  4,? 

9,6349fc-02 

3.96S9E-02 

54.94 

*31.4 

•  1116 

5.3146E-02 

35.06 

2b0.1 

.1249 

6.6432E-02 

23.34 

163.2 

.1382 

7.971VE-U2 

15.04 

116.1 

•  1515 

9.3005E-02 

V.?64 

72.74 

.1648 

.1063 

6.14b 

40.42 

.1781 

•  1 196 

2.171 

17.05 

.IV 14 

.1329 

0. 

0. 

IMPULSE  (PST.MSEO  — 

INCIDENT  «  7.34b 

FEPlECTEOs  57,67 

CAUTION  -CONTACT  SURFACE  MAS  A-RlVtD. 

UA7,»  ARE  CHUL'E  bEYUNU  T(mSEU  aFTER  SHOCK  ARRIVAL*  1.5305E-32 


FIG.  7.8  OUTPUT  RESULTS  FOR  EXAMPLE  9 


=^£5^  S?!5 


NOLTR  72-231 


REFERENCES 


1.  Holland,  N.  0.  (Editor),  "Explosives — Effects  and  Properties  (U)", 
Naval  Ordnance  Laboratory,  NOLTR  65-218,  Feb  1967,  Confidential 

2.  Minzner,  R.  A.  Champion,  K.  S.  W. ,  and  Pond,  H.  L.,  "The  ARDC 
Model  Atmosphere,  1959  (0)",  Air  Force  Surveys  in  Geophysics  No.  115 
(AFCRC-TR-59-267) ,  Air  Force  Cambridge  Research  Center,  August  1959 

3.  Lehto,  D.  L.  and  Lutzky,  M.,  "One-Dimensional  Hydrodynamic  Code 
for  Nuclear-Exolosion  Calculations  (U)",  Naval  Ordnance  Laboratory, 

NOLTR  62-168,  Mar  1965 

4.  Granstron,  S.  A.,  "Loading  Characteristics  of  Air  Blast  from 
Detonating  Charges  (U)",  Transactions  of  the  Royal  Institute  of 
Technology,  Stockholm,  Sweden,  1956 

5.  Jack,  W.  il.,  Jr.,  "Measurements  of  Normally  Reflected  Shock  Waves 
from  Explosive  Charges  (U)",  Ballistic  Research  Laboratories, 

Memorandum  Report  Mo.  1^99,  Jul  1963 

6.  Ericsson,  U.  and  Edin,  K.,  "On  Complete  Blast  Scaling  (U)", 

Physics  of  Fluids,  Vol.  3,  pg«  893,  I960 

7.  Wisotski,  J.  and  Snyer,  W.  H.,  "Characteristics  of  Blast  Waves 
Obtained  from  Cylindrical  High  Explosive  Charges  (U)",  Denver  Research 
Institute,  University  of  Denver,  DRI  2286,  1965 

8.  Brinkley,  S.  R.,  Jr.  and  Kirkwood,  J.  G.,  "Blast  Wave  in  Air 
Produced  by  Line  Charges  (U)",  National  Defense  Research  Committee, 

NRDC  A-343,  Oct  1945 

9.  Fisher,  E.  M.,  "The  Effect  of  the  Steel  Case  on  the  Air  Blast 

from  High  Explosives  (U)",  Naval  Ordnance  Laboratory,  NAVORD  2753,  1953 

10.  Ministry  of  Supply,  '  ti  n  of  Excess  Blast  Pressure  and  of 
Positive  Blast  Impulse  -i.  i  y  for  Cased  Charges  of  RDX/TNT  60A0 
and  of  Mlnol  2  (U)M,  Ar  *  re-  *  .{esearch  and  Development  Establishment, 
1955,  Confidential 

11.  Parks j  D.  K.,  "Blast  Comparison  of  Plastic  and  Steel  Casings  with 
Two  Explosive  Loads  in  the  2.0  Inch  Gimlet  Warhead  (U)",  Denver 
Research  Institute,  University  of  Denver,  1958,  Confidential 

12.  Smith,  R.  and  Wangerin,  C.,  "Some  Preliminary  Results  on  the 
Effect  of  a  Steel  Case  on  the  Blast  from  Pentolite  Cylinders  (U)", 
Ballistic  Research  Laboratories  Memorandum  Report  NO.  863,  195^,  (Conf) 

13.  Heinemann,  R.  V/.,  "The  Effect  of  Casing  Materials  and  Explosive 
Composition  on  Blast  (U)",  Picatinny  Arsenal  Technical  Report  DR-TR  1-60, 
1961,  Confidential 


70 


NOLTR  72-231 


14.  Brode,  H.  L.,  "Reflection  Factors  for  Normally  Reflected  Shocks 
in  Air  (U)",  U.S.  Air  Force  Project  Rand  Research  Memorandum  RM-2211, 
Jul  1958 

15.  "Structures  to  Resist  the  Effects  of  Accidental  Explosions  (U)", 
Army  Technical  Manual  TM  5-1300,  Navy  Publication  NAVFAC  P-397, 

or  Air  Force  Manual  AFM  88-22,  Jun  1969 

16.  Glasstone,  S.  (Editor),  The  Effects  of  Nuclear  Weapons,  Revised 
Edition,  Department  of  Defense,  Atomic  Energy  Commission,  Apr  1962 

17.  Filler,  V/.  S,,  "Explosions  in  Enclosed  Spaces  II:  Measurements  and 
Theory  on  •Static*  Pressure  from  High  Explosives  Detonated  in  Air 

and  Nitrogen  Atmospheres  (U)",  Naval  Ordnance  Laboratory,  NAVORD  3890, 
Jan  1956. 

18.  Filler,  W.  S.,  "Post-Detonation  Pressure  and  Thermal  Studies  of 
Solid  High  Explosives  in  a  Closed  Chamber  (U ) " ,  6th  Symposium 
(International)  on  Combustion,  Reinhold,  N.Y.,  pg.  648,  1957 

19.  Weibull,  H.  R.  W.,  "Pressures  Recorded  in  Partially  Closed 
Chambers  at  Explosion  of  TNT  Charges  (U)",  Annals  of  the  New  York 
Academy  of  Sciences,  Vol.  152,  Art.  1,  pg.  357,  1968 

20.  Marks,  L.  S.,  Standard  Handbook  for  Mechanical  Engineers. 

7th  Edition,  pg.  4-12,  1967 

21.  James,  D.  J.  and  Rowe,  R.  D.,  "Measurement  of  Steady  Overpressure 
Loading  on  the  Wall  of  a  Spherical  Cavity  Resulting  from  the  Detonation 
of  a  Single  HE  Charge  ao  the  Centre  (U)",  Atomic  Weapons  Research 
Esatblishment ,  AV/RE  #E3/64,  1964 

22.  Proctor,  J.  F.,  "Structural  Analysis  of  NOL  Explosion  Testing 
Facilities  (U)",  Naval  Ordnance  Laboratory,  NOLTR  69-84,  Apr  1969 


NOLTR  72-231 


APPENDIX  A 

FLOW  CHART  FOR  COMPUTER  CODE 

The  following  pages  of  this  appendix  give  the  complete  flow 
chart  for  the  computer  program.  It  is  broken  into  three  logical 
sections,  input,  shock  wave  calculations,  and  venting  calculations. 


SUBROUTINE  HEDATA 
gives  equivalent  weight 
and  weights  of  C,H,N,0,AL 
in  desired  explosive 


SUBROUTINE  ARDC 
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SHOCK-WAVE  CALCULATION 


A-3 
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VENTING  CALCULATION  (CONT'D) 
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VENTING  CALCULATION  (CONT'D) 
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APPENDIX  B 

INPUT  DATA  CARDS 

This  appendix  provides  descriptions  and  explanations  of  all 
of  the  information  required  on  the  input  data  cards  for  this  program. 
Formats  of  these  cards  are  given,  and  the  sample  problem  input  cards 
in  Figure  7*1  can  be  used  as  guides. 
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DESCRIPTION  OF  INPUT  DATA  CARDS 


FIRST  DATA  CARD— FORMAT < c5. 2.  IS*  9E5.2,  315) 

NOTE  THAT  THE  THREE  QUANTITIES  INDICATED  BY  — MAY  CAUSE  ADDITIONAL 
CARDS  TO  BE  READ  IN. 

WLB  .WEIGHT  OF  EXPLOSIVE  CHARGE  (POUNDS). 

•NUMBER* I DENT  IF  I CAT  I ON  NUMBER  OF  DESIRED  EXPLOSIVE  IN  LIST. 

IF  THE  DESIRED  EXPLOSIVE  IS  NOT  IN  THE  LIST,  EITHER 

(1)  USE  THE  NEAREST  AVAILABLE  OHE  IN  THE  LIST,  OR 

(2)  ADD  THE  NEW  EXPLOSIVE  TO  THE  LIST,  OR 

(3)  READ  IN  THE  OESIREO  PROPERTIES  AFTER  THIS  CARD. 

RLOD  *LENGTH/D!AMETER  RATIO. 

CASE  *CASE  WEIGHT /EXPLOSIVE  WEIGHT  RATIO. 

VINIT  .INITIAL  VOLUME  OF  CHAMBER  ( .UBIC  FEET). 

AINIT  “INITIAL  VENT  AREA  (SQUARE  FEET). 

PAMS  .AMBIENT  PRESSURE  INTO  WHICH  VENTING  OCCURS  (PSIA). 

TAMB  “AMBIENT  TEMPERATURE  INTO  WHICH  VENTING  OCCURS  (CENTIGRADE). 
ALTKFT* ALTITUDE  (KILOFEET). 

IF  PAMB  AND  TAMB  ARE  BOTH  GIVEN  AS  0..  THE  CORRECT  VALUES  WILL  BE 
FOUND  BY  THE  PROGRAM  FROM  THE  ARDC  ATMOSPHERE. 

ALTKFT  IS  IGNORED  IF  PAMB  AND  TAMB  ARE  NOT  0. 

PCHAM  .INITIAL  AMBIENT  PRESSURE  IN  CHAMBER  (PSIA). 

TCHAM  .INITIAL  AMBIENT  TEMPERATURE  IN  CHAMBER  (Cl. 

IF  PCHAM  AND  TCHAM  ARE  0.,  THEY  ARE  ASSUMED  TO  EQUAL  PAMB  AND  TAMB. 
NOPT  «1  DO  VENTING  CALC.  *2  DO  SHOCK  P-T  CALC. 

•NV  “NUMBER  OF  CARDS  OF  VENTING  CHANGE  DATA  TO  BE  READ  IN. 

•NR  “NUMBER  OF  RADII  AT  WHICH  SHOCK  P-T  DATA  ARE  WANTED. 


SECOND  DATA  CARD.  OMIT  THIS  CARD  IF  NUMBER  IS  POSITIVE. 

THIS  CARD  HAS  TWO  POSSIBLE  FORMS  DEPENDING  ON  WHETHER  NUMBER  IS  0  OR  -1. 

IF  NUMBER* 0,  READ  IN  THE  FOLLOWING  EXPLOSIVE  DATA  FOR  AN  EXPLOSIVE 
NOT  APPEARING  IN  THE  LIST  IN  SUBROUTINE  HEDATA.  FORMAT ( 7E7.2 I 
WFACT  .BLAST  EQUIVALENCE  RELATIVE  TO  TNT  (USUALLY  ABOUT  1.0). 

EFORM  .ENERGY  OF  FORMATION  OF  THE  EXPLOSIVE  (CAL/GRAM). 

WFC  “WEIGHT  FRACTION  CARBON. 

wfh  “weight  fraction  hydrogen. 

WFN  “WEIGHT  FRACTION  NITROGEN. 

WFO  “WEIGHT  FRACTION  OXYGEN. 

WFA  “WEIGHT  FRACTION  ALUMINUM. 

(NOTE  THAT  THESE  ARE  WEIGHT  FRACTION,  NOT  WEIGHT  PERCENT.) 

IF  NUMBER* -1,  READ  IN  THE  FOLLOWING  DATA  FOR  PREPARING  A  MIXTURE 
OF  THE  COMPONENTS  IN  THE  LIST  IN  SUBROUTINE  HEDATA.  FORMAT ( E7.2 ,9 ( 13, FA. 3)  I 
WFACT  .BLAST  EQUIVALENCE  RELATIVE  TO  TNT  (USUALLY  ABOUT  1.0). 

NUMHE  (1 (“EXPLOSIVE  NUMBER  IN  THE  TABLES. 

HEFRAC( l ) “WEIGHT  FRACTION  OF  THIS  EXPLOSIVE. 

NUMHE  ( 2 1 *SAME  FOR  SECOND  COMPONENT. 

HEFRAC(2l“SAME  FOR  SECOND  COMPONENT. 

...  CONTINUE  FOR  AS  MANY  AS  9  COMPONENTS. 
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THIRO  DATA  CAROISI.  OMIT  IF  NV*0.  FORMAT ( 6E7. 2 , 1 7 ) 

THERE  IS  ONE  CARO  PER  N.  THERE  ARE  NV  OF  THESE  CAROS. 

THIS  IS  AN  ARRAY  OF  VENT  AREA  AND  VOLUME  CHANGES. 

PVINJ  “PRESSURE  AT  WHICH  A-V  CHANGE  IS  TO  OCCUR  tPSIA).  IF  THE  INITIAL 
CHAMBER  PRESSURE  EXCEEDS  PV(N),  A  NEW  CHAMBER  IS  ADDED. 

TV( N )  “TIME  AT  WHICH  A-V  CHANGE  IS  TO  OCCUR  (SEC). 

WIN)  “NEW  VOLUME  TO  BE  ADDED  (CUBIC  FEETJ. 

AVI N )  “NEW  VENT  AREA  TO  BE  ADDED  (SQUARE  FEET). 

PAV(N)=AM8I£NT  PRESSURE  IN  NEW  VOLUME  (PSIA). 

TAV(N)=AM8!ENT  TEMPERATURE  IN  NEW  VOLUME  CC). 

NOPTVIN) “CONTROLS  USE  OF  VENT  AREA  AND  VOLUME  CHANGE  TABLES. 

>1  BREAK  INTO  NEW  VOLUME  IF  INITIAL  PRESSURE  EXCEEDS  PV(N). 

«2  BREAK  INTO  NEW  VOLUME  IF  TIME  TV(N)  IS  REACHED. 

“3  BREAK  INTO  NEW  VOLUME  IF  PRESSURE  EXCEEDS  PVIN)  WHEN  TIME 
TVIN)  IS  REACHED. 


FOURTH  DATA  CARDIS).  OMIT  IF  NR“0.  FORMAT  I10E7.2) 

THERE  ARE  NR/10  OF  THESE  CAROS  WITH  10  R  VALUES  PER  CARD. 

TOTAL  OF  NR  ELEMENTS  IN  ARRAY. 

R ( I )  “ARRAY  OF  DESIRED  RADII  AT  WHICH  SHOCK  P-T  DATA  IS  WANTED  CFT>. 
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APPENDIX  C 

DEFINITIONS  0*  PROGRAM  VARIABLES 

A  complete  alphabetical  listing  of  all  program  variables  used 
in  this  code  is  given  in  this  appendix.  Also  a  definition  accompanies 
each  listed  variable. 
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DEFINITIONS  OF  PROGRAM  VARIABLES 


A  *CURRFNT  VALUE  OF  VENT  APEAISC  FT). 
AINIT  =  l  N I T I AL  Vc>NT  AREA  (SC  FT). 

AtRAO^in  f.'OLPS  OF  AIR  IN  NE.vL  Y  ADDE 
AIRVCL=Lo  I-CLEC  OF  .'IR  |f,  ALL  ACTIVF 
ALTKFT 
A  V  t  *1 ) 


ALTtt 

ALT2 

CASt 

DIF1 

OIF2 

DPI 

DO 

DT 

OTE.XP 

DTIVE1 

DUM 

DC 

D1 

EF 

EFOR.M 

EGUIVX 

EQUIV2 

ECmT 

FA 

FC 

FH 

FN 

FO 

FLEFT 

FRAC 

0 

GASL8 
GAS.'OL 
G  G 

GRAV 

GO 

Gl 

G2 

HEFRAC 

HF 

I 

IR 

J 

JJ 

K 


CHff-’sr.R. 

CHAMPERS. 

AL T I TUDt t A l LOr EE)  )  .  .\UT  USED  IF  PAMj,  TAMtJ  ARE  GIVEN. 
ARRAY  OF  MEW  VENT  AREAS  I  SO  FI). 


P. 

P. 


OLJPUTE  .UAL  ALT  1  T'j~t(  fit  I  )  . 

ALTITUDE (F'ETiRC  I  A^LVC  MEAN  SEA  LEVEL. 

CASE  WEIGHT /L-lARij;:  .vEIoriT  RATIO. 

PREVIOUS  OIFFTRFNCE  BETWEEN  ?  WAYS  TO  CALC 
CURRENT  DIFFORE'JC-  BETWEEN  2  WAYS  TO  CALC 
INTEGRATION  INTERVAL  IN  PRCSSURE ( PSF  )  . 

ENCRGY  ADDED  KCAL )  IN  INTEGRATION  STEP. 

INTEGRATION  INTERVAL  IN  FINDING  TEMPERATURE. 
TEMPERATURES)  INTERVAL  IN  T.P.G  ITERATION. 

TIME  INTERVAL  IN  VENTING  INTEGRATION  (SEC). 

DUMMY  VARIABLE  IN  SUBROUTINE  CALL. 

DENSITYILB/CJ  FT!  AT  START  OF  VENTING  STEP. 
DENSITY(LS/Cj  FT)  AT  END  OF  VENTING  STEP. 

ARRAY  OF  HE  ENERGY  OF  FORMATION  DATA  (CAL/G). 

ENERGY  OF  FORMATION  Or  DESIRED  EXPLOSIVE  0CCAL/G) . 
INTERMEjIATc  GUANTITY  IN  CHARGE  SHAPE  CORRECTION. 
*RRAY  Cr  CYL-SPH  LUUIVALtfl'.t  FACTORS. 

EQUIVALENT  WEIGHT  REFeRRuu  TO  TNT  (FOR  SHOCK  CALCS). 
WEIGHT  FRACTION  AL. 


OF 

OF 

OF 

OF 

OF 


ARRAY 
ARRAY 
ARRAY 
ARRAY 
ARRAY 
LOSS  FRACTION 
INTERPOLATION 


WEIGHT 

WEIGHT 

WEIGHT 

WEIGHT 


FRACTION 
FRACTION  I 
FRACTION  I 
FRACTION  : 
FOR  VuNTIrtO 
FACTOR. 


•.ASS  CHANGES. 


CURRENT  VALUE  OF  SPECIFIC  HEAT  RATIO. 

POUNDS  OF  GAS  REMAINING  IN  ACTIVE  CrIAviiRS. 

LB  MOLES  OF  GAS  REMAINING  IN  ACTIVE  CHAvgERS. 

( fi-1  .  )  !'• 

ACCELERATION  OF  GRAVITY  *32.2  FT/SEC/SEC. 

SPECIFIC  HEAT  RATIO  AT  START  Ur  INTEGRATION  STEP. 
=  SPEC IF'I C  HE/T  RATIO  AFTER  MIXING  GASEj. 

SPCC I F I HEAT  RATIO. 

ARRAY  OF  WEIGHT  FRACTIONS  (USED  WITH  NLKHE). 

ADC REV J  AT  I CN  FOR  HFFRACII). 

GENERAL  DO-LOOP  INDEX. 

INDEX  FOR  l-C-LCCP  ON  R(IR).  RANGE  IS  1  TO  NR. 
DO-LOOP  INDEX 

DESIRED  TNT  DATA  LIE  oETwt'EN  RIJJJ  AND  RIJJ-1). 
DO-LOOP  INDEX 


57 V  W»* ‘W' 
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UMAX  1  -DESIRED  NUMoER  CF  P-T  POINTS  oEFORE  HE  GASES  ARRIVE. 

TYPICAL  VA!  UKS  ARF  10  TO  40. 

KMAX2  -VENTING  PRKIOUT  CJ.sI,wL.  PRli.T  AtCUT  100/K.MAX2  LINES  OF 
DATA.  TYPICAL  VALJZS  ARE  2  TO  10. 

KOUNT  -LUUNTtR  FOR  PRINT  i.vG  POKING  Vt.NI  INvj. 

L  -DO-LOOP  INOFX 

MCO  -ED  VCLES  OP  C"  I“!  THE  CHAVBFR. 

MC02  -LB  MOLES  OF  CC2  Pi  THE  CHAFFER. 

MH2  -LR  -'OLES  OF  H2  IN  THE  CHA  -‘ER. 

VH20  «LB  MOLES  CF  H2C  IN  THE  CHAMBER. 

MN2  -LB  MOLES  CF  N2  P.  TML  C11AM1IER. 

M02  «LB  MOLES  OF  02  IN  THE  CHAMBER. 

Mi  -GRAM  MOLES  OF  C02  FORV-tD. 

V2  «GRAV  VOLES  OF  H2  PQP-'EC. 

M3  «NOT  USED. 

"4  -NOT  USED. 

MS  -GRAM  VOLES  OF  AL203  FORMED. 

M6  -NOT  USED. 

M7  -NOT  USED. 

M8  -GRAM  VOLES  OF  H20  FORMED. 

M9  -GRAM  VOIES  CF  CO  FCRV£0. 

NAME  -ARRAY  FOR  NAVES  CF  EXPLOSIVES. 

NAMES  -ARRAY  FC9  NAVE  OF  DESIRED  EXPLOSIVE. 

N  -EXPLOSIVE  NUMBER  IN  TA3LE.  OR  CHA«SCR  BREAKING  INDEX. 

NN  -INDEX  FOR  VENTING  LOOP. 

NEON  -1  FOR  P.GT.PC.  >2  FOR  P.lT.PC  (Cr.OOSES  VENTING  EQUATION). 
NOPT  -1  DO  VENTING  CALC.  *2  UC  SHOCK  P-T  CALC. 

NOPTVIN ) -ARRAY  OF  WALL-BKEAKING  OPTIONS. 

-1  BREAK.  WALL  IF  PRESSURE  EXCEEDS  PVIN). 

-2  BREAK  WALL  WHEN  TIMC  TV(N)  IS  REACHED. 

-3  BREAK  WALL  IF  PRFSSURE  EXCEEDS  PVIN)  AT  TIME  TV(N) • 

NR  -NUMBER  CF  ELEMENTS  IN  USE  IN  S  ARRAY  II  TO  100).  ' 

NSAVE  -LAST  LINE  OF  VENT  DATA  USED  IN  INITIAL  BREAKS. 

NUMBER-NUMBER  OF  OtSIREO  EXPLOSIVE  IN  uATA  LIST. 

NUcMHE  -ARRAY  Or  tXPLOSIVc  NUMBERS  FOR  mRLTRAKY  MIXING  UP  OF  HE. 

NV  -NU.M-ER  OF  ELEMENTS  IN  VENTING  ARRAY. 

Ni  -GRAM  .iOI.ES  OF  C  IN  THE  EXPLOSIVE. 

N2  -GRAM  VOLES  OF  H2  IN  THE  EXPLOSIVE. 

N3  -GRAM  "OLES  OF  N2  IN  THE  EXPLOSIVE. 

N4  -GRAM  VOLES  OF  02  IN  THE  EXPLOSIVE. 

N5  -GRAM  MOLES  OF  AL  IN  THE  EXPLOSIVE. 

N6  -GRAM  VOLES  OF  N2  IN  IHE  CHAMBER  AIR. 

N7  -GRAM  MOLES  OF  02  IN  IHE  CHAMBER  AIR. 

0VEWP2* ARRAY  OF  CYL-GPH  cQUl VALENCt  CVERPRlSSURES  (PSI). 

OVPSI  -OVERPRESSURE ( PS l )  AT  START  CF  VENTING. 

OVPSI 1-CURRENT  OVERPRESSURE ( PS f )  • 

OVPO  - -OVERPRESSURE  IPS  I)  AT  STAPT  OF  VENTING  STEP. 

0VP1  -OVERPRESSURE I P 5 1 )  AT  END  CF  VENTING  STEP. 
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P<  !  i 
PA 

PAMB 
PAV(-N) 
PC( JJ) 
PCHAM 
PC  TNT 
PI 

PINIT 

PSCAL 

PSFAMp 

PSFCH 

psnn 

PSIFEr 
PTNT 
PTSCAL 
PV  ( N) 
PO 

POPSI 

PI 

P1PSI 

P2 

P2A 

P2B 

0 

OPERG 

01 

R 

RC  II 
RC  y. 

REF 

RESULT 

RLOD 

9R 

RSCAL 

RTNT 

SIGMA 

T 

TAMB 

TAU 

TAVI N  t 

TCHAM 

TCI  JJ) 

TCMSEC 

TCTNT 

TEMP 

TEMPO 

TEMPI 

TEMP  2 

TF 


■ARRAY  OF  INCICCNT  OVCRPR-.SSURG (PS I  )  FOR  TNT. 

•  INITIAL  i:vERPRES5.XE(PSn  IN  CH.*.:-:;.FR  WITHOUT  GAMMA  CORR. 

■  OUTSIDE  A.-.filENT  PRESSURE  I  ®S!  A )  . 

■AR1 AY  OF  AMBIENT  PRESSURt S< PS1 A )  IN  NEW  CHAMBERS. 

■ARRAY  OF  Lu-.TACT  SUNFACc.  ARRIVAL  TIMEISECl  DATA  FOR  TNT. 
■INITIAL  AIR  PNESSUKlIPAIAJ  IN  ORIGINAL  CHAVdER. 

■  CONTACT  SURFACE  OVlRP.RESSURE (PS  I )  AT  RTNT. 

■  OVSRPR/PtA*.  OVCNP  RATIO  FROM  FITTING  cOUAT  tON. 

■  INITIAL  CMA'VCR  PRESSURFlPS^A)  BEFORE  WALLS  aREAK. 

■SC' L I NG  FACTOR  FOR  REDUCING  PRESSURES  TO  SEA  LEVEL. 

■OUTSIOF  AMP If NT  PRFSSURE ( PSFA) . 

■CHAMBER  PRESSURE (PSFA) . 

■  AVAY  OF  OVERPRESSURE (PS  1 1  AT  DESIRED  RADIOS  R  AT  TIME  Tl(I). 
(  I  )*REFLlCT-^  OVLRPRlSSUREIPSI I  CDRRESP.  TO  PSI(I). 

■  PEAK  CVERP<£SSURE(PSI )  AT  RTNT. 

*PSf al*tscal 

■ARRAY  CF  PRESSURES(PSIA)  WALLS  CAN  WITHSTAND. 

«PRCS5URE ( PSFA )  AT  «=TART  OF  VENTING  STEP. 

■PC  IN  OVERPRESSURE (PS I  I. 

■PRESSURE ( PSF A )  Ai  END  OF  VENTING  STEP. 

■PKESSURE(PSIA)  AT  tl.J  OF  VENTING  STEP. 

■PRESSURE(PSFA)  AFTER  .ALL  BREAKS. 

■PRESSURE  FROM  FIRST  EON  (FOR  MIXING  T*0  CHAMBERS) » 

■  PRESSURE  FROM  SECOND  EON  (FOR  -'*1X1  NG  TWO  CHAMBERS). 
■ENERGY(KCAL)  RELEASED  3Y  EXPLOSION. 

■ENERGY  RELFiSFOItCCAL/GSA*-'). 

■CUMULATIVE  ENERGY  DURING  INTEGRATION  FOR  TEMPERATURE. 

■  GRAM  MOLES  02  LEFT  IN  CrtAKcER  (CALLED  RR  IN  BLAST),, 

■  ARRAY  OF  L-LoIREd  L  I srANCLS (  F  T  )  .  NR  ELEMENTS  IN  THIS  AKRAY. 
■DISTANCE  R ( I )  CONVERTED  TO  CM. 

■OVERPRESSURE  REFLECTION  FACTOR. 

■CJANTITICS  BEING  PRINTED  IN  SUBROUTINE  GASES. 
■LE'IGTH/OIA.vetER  RATIO  OF  CHARGE. 

■(IN  BLAST  ONLY )  GRAM  MOLES  CF  C2  LEFT  IN  THC  CHAMBER. 

■  SCALING  FACTOR  FUR  RE.JUCl.NG  RADII  TO  SEA  LEVEL. 

■R< IP )  RtU^CEU  TO  1  Lo  TNT  AT  SEA  LEVEL. 

■  SHAPE  PARA.’-ir'.TtR  IN  MTTIuG  tOUATJCN  FOR  P-T  DATA  FOR  TNT. 

•  T£vPtRAT'jRE(  R  )  IN  SUBROUTINE  GAFMA. 

■OUTSIDE  AMBIENT  TEMPER A TURl I C )  . 

■FRACTION  OF  POSITIVE  DURATION  RASED  ON  65.  CM  CURVE  SHAPE. 

■  ARRAY  OF*  AMBIENT  TEMPS(C)  IN  Nf»  CHAMBERS. 

■  INITIAL  TEMP(C)  IN  ORIGINAL  CHA  -liiEtR  BEFORE  EXPLOSION. 

«ASR AY  OF  TIME(SEC)  GET *Es.N  SHOO.  AND  CS  ARRIVAL  FOR  TNT. 
■TCTNT  IN  MSEC. 

■CONTACT  SURFACE  ARRIVmL  TIME<SEC)  AT  RTNT. 

■  TE-'PERA TUNE  IN  ARCC  GvoKuuTINc. 

■  GAS  T  c”.P  ( R )  AT  START  OF  VENTING  STEP. 

■  GAS  TC.''P ( R )  AT  END  CF  VENTING  STEP. 

■  GAS  TfJf'fR)  AFTER  NEW  CHAMPFR  IS  ADDED. 

■GAS  TEMPERATURE(FAhRFNHCIT). 
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r  1  me  1  *Ti'-,£cr>ECi  at  of  venting  stop. 

TOTMOL*TOTAL  PO-j.'.D  MOLES  OF  GAELS  IN  CHAMBER. 

TPt JJ)«ARRAY  OF  POSITIVE  OVERPR  JURAT ICNI SEC )  FOR  TNT. 

TPTNT  «POSt  T  IVC  Phase  O'JnAT  I JN ( StC  )  AT  RTNT • 

TR  ‘GAE  TFVPIU)  IN  CHAMBER. 

THCH  «A:-wILNT  Ib.'-.PI.R)  !.  L-RIGI..mL  CriA.-oLR  UEKJRE  EXPLOSION. 

TS ( JJ ) *ARRAY  OFSHUCx  FRONT  ARRIVAL  TP’CCSECJ  i/ATA  FOR  TNT. 

7SCAL  ‘SCALING  FACTOR  FOR  REDUCING  TIMES  TO  SEA  LEVEL. 

TSTNT  «SHCC<  F;<"iT  ARPIVAL  TIMF.ISEC)  AT  RTNT. 

TV(N)  rAPPAY  OF  TIMfStsrr)  FOR  i. PEAKING  WALLS. 

Till)  ‘ARRAY  OF  TIME (SEC)  AFTER  DETONATION. 

T2II)  ‘ARRAY  OF  TIME! SEC  I  AFTER  SHOCK  ARRIVAL. 

U  ‘SPECIFIC  MEAT  OF  GAS  MIXTURE. 

01  ‘SPECIFIC  MEAT  OF  CC2. 

02  ‘SPECIFIC  MEAT  OF  M2. 

UA  ‘SPECIFIC  HEAT  OF  02. 

UA  ‘SPECIFIC  HFAT  OF  N2. 

UA  ‘SPECIFIC  HEAT  CF  H2C. 

UR  ‘SPECIFIC  HEAT  OF  CO. 

V  ‘AC1IVE  VOLUME (CU  FT). 

VINIT  ‘INITIAL  CHAMBER  VOLUMEIC'J  FT). 

VV(N)  ‘ARRAY  OF  NE.V  CHAMBER  VOLUMES < GU  FT). 

VO  ‘CHAMBER  VOLUMcICU  FT)  AT  START  OF  VENTING  STEP. 

VI  ‘CHAMBER  VOLUME ( CU  FT)  AT  END  OF  VENTING  STEP. 

V2  ‘CHAMBER  VOLUME  I CU  FT)  AFTER  tiCW  CHAMBER  IS  ADDED. 

WA  ‘POUNDS  CF  * L  IN  THE  EXPLOSIVE. 

WC  ‘POuNDS  OF  C  IN  THE  EXPLOSIVE. 

WFA  ‘KEIC-hT  FRACTION  OF  AL  IN  THE  EXPLOSIVE. 

wfact  ‘Charge  energy  relative  to  eoual  ..eight  of  tnt. 

wFC  ‘WEIGHT  FRACTION  CF  C  IN  THE  EXPLOSIVE. 

WFH  ‘/.EIGHT  FRACTION  Or  H  IN  THE  EXPLOSIVE. 

WFN  ‘WEIGHT  FRACTION  OF  N  IN  THE  EXPLOSIVE. 

WFO  ‘■.•.EIGHT  FRACTION  OF  0  IN  THE  EXPLOSIVE. 

WFT  ‘APPROXIMATE  CHARGE  LtNGTnlFT). 

WH  ‘POUNDS  CF  H  IN  THE  EXPLOSIVE. 

WLB  ‘WEIGHT < Lb )  OF  EXPLOSIVE  CHARGE. 

WLBA  ‘ADJUSTED  CHARGE  /.EIohT(Lu). 

WN  ‘POUNDS  OF  N  IN  THE  EXPLOSIVE. 

WO  ‘POUNDS  CF  0  IN  THE  EXPLOSIVE. 

’/.•PERL  ‘APPROXIMATE  CHARGE  WEIGHT  PER  UNIT  LENGTHCLU/FT ) . 

XCASE  ‘CORRECTION  FACTOR  FOP  CASE  EFFECT. 

XJMP1  ‘.SIDE-ON  POSITIVE  IMPULSE  BEFORE  HE  GAS  ARRIVAL  IPSI.MSEC). 
XIMP1R*REFLECTEJ  OVLRPR  POSITIVE  IMPULSE  aEFORE  GAS  ARR  (PSI.MSEC). 
XSHAPE‘CCRRCCT ION  FACTOR  FOR  CHARGE  CYL INDRICI TY. 

XSHAPL«CORRtCTION  FACTOR  FwR  INFINITE  LINL  CHARGE. 

X2  ‘GRAM  MOLES  OF  GAS  IN  IN  CHAM5ER. 

X3  ‘GRAM  MOLES  OF  GAS  IN  CHAMBER  WITHOUT  THE  AIR. 
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APPENDIX  D 

FORTRAN  LISTING  OF  PROGRAM 

This  appendix  gives  the  complete  FORTRAN  listing  of  the  computer 

« 

program.  All  seven  sections  are  labeled  with  appropriate  card 

i  ' 

numbers . 

/ 

BLAS0010  -  BLAS3160 
MIX  0010  -  MIX  0270 
HEDA0010  -  HEDA1U00 
GAKM0010  -  GAMM0160 
GAS  0010  -  GAS  0890 
TNT  0010  -  TNT  1330 
ARDC0010  -  ARDC0H10 
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FORTRAN  LISTING  OF  PROGRAM 


PROGRA''  BLAST)  I'lPOT  ♦OUTPUT) 

COMMON /DAT  41 /nL'>..Nq.'-’  •  R  .  R  U  U  ,  C  A  S  E  .  V I .  >  t  T  .  I N  l  T  .  P  A  M :  .TAMB.ALTKFT . 

1  PCHAV.TCSV.f  v?>T.  .V.iYi.  '•FACT.uCLKM 
COMi’iUiN/ TM  | *./R  TuT  */.i  A XI 

CU.-'.-.vf./TNt..LT/F  r..T.rul..T.tP).,T,lLT..T.PCT;.T.PliI(4G>.TlC40I.T4:<4C) 
l.JJ.Xt'Pl 

*.  t.(t  •<  C v  •  *  C  u2  ti  i2q  1H2 
COMMON /  H  E-V  453/  v.C  . .-.•4..iN«»w ,  ..A 

common / 'x T f 1? f r * v. e*-,* r.N,.- ro«  fa 

COMMON/ GAS/M  •  =  •2  ».\3»\4»:«5...  6.  NT. 

1  mi.m2.mj. a, r?R«o«x2 

common/ vent/  pvi  10  *rv«  ;*))  »vvc  i'll  .avii^i  .pavi  1 0 >  * t v ( 1 0 *  ,.»uPTvt  io 

lt*i 

COWSON/HE/w  .Hl<9)  »HEF\Al.(9) 

DIMENSION  PSlP’-r  (40I.UI  ICC) 

C 

OtMFNSIO':  OVFpP2(24|,COUtV2(24) 

PEAL  ‘-’2  «?  x«vA 

REAL  .•••0? »"N2»”C0 .  •••CO 2 ,,'H?''.XH2 
r 

C  TABLES  FOR  CYLIMRICAL  Vj  SPHERICAL  CHARGE  EQUIVALENCE. 

DATA  ECU  I V2 '  5. 12  ,4.  aO  »4»CC  »■>»  po»3 .16 .2.  Jy  »C.**o»2. 17  *  1.93  » 1  •  7 1  * 

1  1.53.1.35.1.20*1 .05. .80.. 45,. 24, .158.. 085.. 037, 

2  .042. .04'',. "iS0. .OS')/ 

DATA  0VERP2/«P»;»  »3<’.0»  .400.  .50C..6CG*  »?CC»  *800.  *900.  *  1000.  ,1 100.  » 

1  lPOO,,^",,!  a*".  ^OCi  ♦?0('0.  » 

2  60no..55*f'.,7r>')'»..l.E6/ 

r 

300  FORMAT (C5.2.15.9E3.2.3I5) 

305  FOR.-.AT(l7. 13. F4.j|) 

310  FORVAT( 10E7.1) 

320  FORMAT (*0FHCC<  WAVE  CALCULATION*/ 

1*CINPUT  PARAMETERS  »,  16X,  *CrlARGE  ivEIGHT  ADJUSTMENTS*/ 


1*CINPUT  PARi-iETCRS  » »16X,  *CrlARGE  i-.EIGHT  ADJUSTMENTS*/ 

2#  CHARGE  V.f  I'--HT(Li-)  **C12 .4.4X. *A1JUGTED  .-.TIL  t*  TNT)  **G12.4/ 
3*  EXPLOSIVE  NUMBER  «*I6.10X»  #HE  LMERGY  FACTOR  «*Glc.4/ 
4*  L/i)  RATIO  *«G12.4»4X.*ChARGE  SHAPE  FACTOR  «*uU.4/ 
5#  CASE/CHARGE  /.T  RATIO  *«ul2 .4.4A.*C/Uu  ..EIuHl  FACTOR  **GU.<f/ 
6*  CHAMBER  PRESSURE  t  PS  I A  )  **ul2. 4. 4X. ‘PRESSURE  SCALE  FACT'  .<*#Gic.4/ 
7*  CHA»ER  TE-'PIC)  »*Gl2.4.4X.*t>f3TANCc  3CALC  FACTOR** 312,4/ 


S*  ALTITUDE  (K.FT)  »*G12.4,4X,*T  |;.£  SCALfc  FACTOR  «*G12.*»/ 

0  36X,  4X,#N0Pi'AL  HEFL  FACTuK  *#G12.*») 

U?r>  FORMAT  (»1  INTERNAL  uwAST  DAMAGE  Kf^wAf, ISMS  PROGRAM,  MAR  1972#) 

430  FOR*1AT<#09C3IRED  **  l  ST  ANCF;  ( FT  |  -#1PG12.4/ 


1# 

450  FORMAT ( 

1*0  TIME  AFTER  TI.-.C.  AFTER  I.NCIucNT  NORM  RtF’L  */ 

2*  EXPLOSION  3HuC<  ARK  LVEnPRESG  OVERPKESS  »/ 

3*  (MSEC)  (.•'SEC)  (PSD  (PSD  *) 

AOG  RCIO^AT  ( 1X,1°4G12,4) 

500  FORMAT)  #«fvs'.LSE  IPSI  .MSEC )—  */ 

1*  INCIDENT  =*IP012,4/«  REFLECTCC*#1PG12.4| 

52')  FORMAT (  •OCA.'TJuN — COhTACT  SURFACE  HAS  ARRIVED.*/ 

1*  DATA  APE  CRUDE  uCYOND  T(r.jiC)  AfIlK  SHOCr*  ARKI  VAL**1PG12.** ) 
533  FORMAT  ( *  Cr'-RG-  3H„Pc  iOkKlCT  ION  IS  uRUDl.  P3I  EXCtEOS  •/ 

1  «  RANGE  OF  EXPERIMENTAL  DATA.# ) 

540  FORMAT ( *  CASE  /.EIGHT  CORRECTION  IS  CMJuE,  PS)  EXCEEDS  */ 


**1PG 12  » 4 ) 


NCT.i-i  rvtFL 
OVERPKESS 
(PSD 


uLASOOlO 
cLA 30020 
OLA30030 
uLAiOO  «0 
qL-oOOpO 
ivLAsOObC 
oLA5G07Q 
DLASOObO 
HLA30090 
uLASOlOO 
SLAS0110 
)  BLA30120 
oLAoOlSO 
DLAS0135 
BLASOIhO 
oLASOl 53 
6LASU160 
SLA30170 
9LAS0160 
5LAS01 °C 
BLASO^OO 
BLAS0250 
DLA50252 
GLAS02  34 
BLASC260 
3LiSf)262 
OL-S0264 
8LA3029C 
6LA50300 
BLAS0uC3 
BLAS03 10 
oLASO 320 
C.LAS0330 
alA503*0 
bLASCuaO 
OLA30360 
oLA 30370 
aLASOidO 
ULA30390 
dLASO'.OO 
ruLA304l0 
6LiS0*»20 
BLAS0430 
bLAS04**0 
BLA30450 
oLAS0*»o0 
t>LA30**70 
t)LA30480 
PHS0490 
BLA 30300 
BLAS0510 
OLASC520 
uLA 30330 
tiLASOaub 
BLA30537 
BLA30P40 
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1  *  RAN?"  OF  rX°fRt**rM'-L  f'ATA.at 

KMAX 1 « 1 0  «*■.•■*  A  X2«lC 

READ  INPUT  •'AT-. 

570  READ  300*  >‘L"  «Nv--»;Lrt .RLGO.C Aj- , V I <« I  T  *>. INI  T *PA.a<o*TAIii.*ALTNFT  » 

1  PCHAp  •TCHAA'(.‘4uP  T  .KVtuR 

IF  ( t.Ltl.LE.'*.  )  jK'P  SPRINT  -2C  SIF C KONuER.tU.O)  GO  TO  600 
JF  CF.UW,,ER.'  E.-l  I  00  TO  5  '* 

READ  305.  .FACT, iN-J, vhp (  t  )  »H--FP.‘C  t  ! )  .  I*t  »«> 

5«6  CALL  HEOAT.'.  SCO  TO  605 

600  READ  3I0.V.FACT.CFO«W...FC,v»  H..-.I  N..-.F  C..- FA  AtFOR.i«EFuRM/1000« 

C  FOUNDS  OF  EACH  ILEF i'«T 

60*  *'C«'..FC*«Li:  ».»r»«.«FH**L!»  i*i\«,.r..*«Lcj  .»mU»  *FO**Lii  i»A»lr.F'A**Lo 

IF ( aFACT.GT .0. )  00  TC  615 
■••■FACT.!.  ‘PRINT  612 

612  FORVATI*  •.-.•FACT  f.CT  l.'  IS  USED.*) 

615  PRINT  620  »  NUNB  ER  .  *F  AC  T  .  C  FOR’-'.*  »FC  .  .<F  h  » :lf N  .  *F 0  •  i.  F  A 
620  F03-t4Tc«3£xPLCS!Ve  PROPERTIES*/ 

1*  NUMBER  E.'^T  EFOKrl  tXPL-GIVE  COMPOSITION  oY  ..EIGHT*/ 

2*  <CAL/G  C  H  M  0  AL  */ 

31H  »l4,F7.3.Ffl.6,3F6.3) 

C  FI  NO  PA'^  AN!)  T  IF  NCT  GtVF.N. 

IFtPAY?.EC.:>.l  CALL  Afi!'CCALT<FT .PAf'.fi »TAVP. ) 

IF (PCHAK.LE.O,  I  PC»-A.v«PAMP  SIFCTCHAV.LE.O.)  TCHAf-UTAMB 
C  DO  VENTING  CALC  IF  ?.CPT«1  AM)  1>C  SHOCK.  ..AVt  CALC  IF  NGPT-2. 

IFCNOPT.LT. 2)  GO  TO  1540 
C 

C  itGIN  3H0CN-..AVE  PKCPc.RT I cS  SCCTICN. 

730  READ  310*(R(IR)*IK*1»NR) 

00  1463  IR.l.AR  MFC  13. £0.1)  GO  TO  760 

PRINT  420  JPHINT  620».\U!oER,».FACT.CF0RM,VFC.a'FH,v.FN.WF0»WFA 
760  RC«*RC IR)*?0.4P 
C  ADJUST  CHARGE  A'lGriT. 

WLBA«.VLO*vFACT  tPA.".6«PCHA"  JTAi-'b»TCHA.V 
RSCAL- ( 1. /■•  LCiA*PAMI/14. 6961  ?*)*•.  333  J3333 
RTNT»RC-'*RSCAL  *CAlL  TNTCCI 
C  YAXE  CHARGE  SHAPE  CORRECT ICN. 

XSHwPE-»l.  S1KRL0J.LT. 2.  I  C-0  TC  1140 

IF( PS IClI.GT. 100.1  PRINT  535 
C  CHARGE  HEIGHT  PcR  UMT  LENGTH  JF  CYLINDRICAL  CHARGE. 

WPERL«C3.!414*100./(4.»RLOJ**2I}«*. 33033333  **LUA**. 66666666 
C  MAKE  CHARGE  SHAPE  CORRECTION. 

IFlPStm.GT.20.)  GO  TO  900  iXSHAPt«l .4 5  SGO  TO  1140 

900  XSMAPE*0»6! 3*P5 1 1 1 ) **.26  7 
C  FIND  INF  IN  I Tt-CYL I HvuK  CHARGE  S"APt  CoHRtCT ION. 

IFCPSim.LT. 200.)  GO  TO  11*0 

DO  940  I«2.24  SIFCP51CD.lt. 0VERP2CI))  GO  TO  950 
940  CONTINUE  «f»24 

950  FRAC* IPSI C I J-0VEPP2C I- l ) 1/C0VERP2C I )-0VERP2C 1-1 ) ) 

ECUIVX.E0UIV2C 1-1 ) .FRAC*  <  LwOl V2 1 1 I-L0UIV2C l-l ) ) 

XSHAPL«E0UIVX*<.PEKL**1.5/WLr.A 

tFCXSHAPL.LT. XonAPL I  XSHAPl«X5mAPL 

C  TEST  IF  otSIRCS  JlsTANct  IS  CLCGt  L..O0GH  TO  CHAKGt  FOR  GOOD  RESULTS. 
C  APPROX.  CHARGE  LENGTH. 


I’L  AS0342 
r>LASC360 
ULAS0P62 
6LAS0P65 
CLAS0370 
6LAS0P72 
bLAS0374 
3LA50576 
6LAS0580 
CLAS0596 
bLASObOO 
0LASO0C2 
LLASC605 
ULAS0608 
3LA50610 
GLAS0612 
GLAS0615 
bLAS0o20 
SLAS0625 
OLAS0630 
BLAS0640 
BLAS0650 
BLAS0660 
BLAS0670 
BLAS0690 
GLAS0700 
ULAS0710 
BLASU720 
bLAS0730 
BLAS0740 
OLAS075C 
BLAS0760 
DLAS0770 
6LAS0700 
tcAscaoo 

oLASOdlO 

ULAS0630 

bLAS0B40 

bLAS0B50 

GLAS0660 

0LAS0870 

bLAS0680 

□LAS0690 

ULAS0900 

LLAS0910 

OLAS092C 

BLAS0930 

BLAS0940 

6LAS0950 

BLAS0960 

ULAS0970 

bLASORbO 

ULASIOIO 

BLAS1020 


r\  r\  n  r 
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l"?**  »FT«WPF9L*‘U* 

f  VA<F  CASE  ••FI-r.MT  CORRECTION.  CASE«CASf /charge  weight  ratio. 
1140  XCAS£«0.477''.53/(1.7CA6t) 

IF(CASt»LT.0.53)  XCASG*i.-CA:ie**2/(  i.+CASEI 
WL?A«/(L-3*»XSf-AP£*XCASF 
IF!PSH1).GT.100.)  PRINT  340 
C  FIND  SCALING  faptgrs. 

PSCAL*I4.AOM7*/PAvp 


RSCAL*!  ’•/'.•L~A»PAKH/14.60M7R)»*.33?34*33 
TSCAL«SSCAL*SORT(  <273.167TA.MA|/288.16) 

C  FIND  DESIRED  RADIUS  FOR  1  LB  TNT  AT  SEA  LEVEL. 

RTNT»RC  »«SCAL 
C  FIND  DATA  IN  TA-  LES. 

CALL  TNT (  1  I 

C  CALCULATE  NORMAL  REFLECTION  FACTOR. 

IFIPSim.6T.2nG,)  GO  TO  1280 

REF«!7.«14.S'>417374.*P5I!1  )  )/<7.*l4.6G61787PSt(l)  )*2.»G0  TO  1290 
1280  RFF«-3.i8.?,37*ALOGlO(PSI(in  l IF ( REF.GT.13. )  REF»13. 

12«o  xi;apip»xI':o!.ref 

C  SCALE  RESULTS  TO  DtSIREO  ChARGc  ocIGhT  AND  mMoIENT  CONDITIONS* 

C  CHANGE  TIMES  FRC"  SECONDS  TO  MSEC. 

00  1340  <«1  .XMAXt 

PSI!<I*PSI  (M/PSCAL  SPSIH£FU)«PSlCO#REF 
Tl(K»«Tl(K)/TSCAL*10no.  ST2 ! < ) «T2 U ) /TSCAL*1000* 

1340  CONTIMi'F 


PTSCAL«PSCAi.*TSCAL 

XI MP1*X IMP 1*1000. /PTSCAL  SX!MP1R*XI.MP1R*1000. /PTSCAL 

C  PRINT  THE  RESULTS. 

PRINT  320»WLa..vLbA.  f*o.v.3ER»*FACT *  RLOO.XSHAPE.  CASt.XCASt. 
1  PAMO.PSCAL.  TA.'-io .KSCAL »  AlTaFT  .TSCAL.  REF 
PRINT  43n,P(  ir >  ,RC'M  SPRINT  450 


PRINT  4«>A«,  (Tl(K).T2(4)(P4I(N).PStREF(<).K«l  .KM  AX  I  » 
PRINT  500, XI MP1.XIMP1R  3TC’>SEC«TCTNT»1000./ TSCAL 
IFcRTNT.LT. 65. I  PRINT  520.TCMSEC 
1480  CONTINUE 
C-C  TO  570 

’  END  SrtOCK-rfAVE  PROPER UtS  SECTION. 


BEGIN  VENTING  SECT t ON. 

1540  PRINT  1550. ‘-.LP. VIM  T. AIN!  T.PA.VB,TAM9,PCMAM,TCHAM.NCPT.NV 

1550  FORMAT! *OVEMt.NG  CALCULATION*/ 
l*OCHAR«iE  Wc!GmT(LB)  **012.4/ 

3*  IMT  VOLUME ( CU  FT)  *»G12.4/ 

4*  INIT  VENT  AREA! SO  FT)  »*G12.4/ 

5*  AMJIENT  PRESSURE! PSIA)*«G12.4/ 

6*  AvntrM  T*'<p ( C I  ■*0)2,4/ 

7*  CHA‘*rro  r*r<sijoF(  f»«;  1  a  )  ■*C-1 2 .4/ 

F*  CHA-NER  T F MP ! C )  *«GI2.4/ 

9*  NCPT** 1 3  »*  KV«*t?) 

1650  FORMAT ! 6E7.C. l 7 ) 

1660  FDR!  AT!«OTA-LE  CF  VOLUME  AND  VENT  AREA  CHANGCS*/ 

1*  P(PSIA)  T(StC)  V ! CU  FT)  A(SU  FT)  PAMBCPSIA)  T 

2"3 ! C )  NOPTV*/10(1H  ,6012.4,1?/)) 

1 6CC  FORMAT ( *0?C31N  VENTING  CALCULAT ION* J 


E-LAS1030 
CLASl 130 
fcuASl 140 
SLASH  50 
C-LAS1160 
ELAS1162 
FLASH  70 
PLASllBO 
6LAS1190 
SLAS1200 
ELASU10 
oLA31c20 
ULAS1230 
CLAS1240 
bLASl/30 
BLAS1260 
5LAS1270 
CLAS1280 
6LAS1290 
bLASlAOO 
6LAS13 10 
FLAS1320 
BLAS1330 
6LAS1340 
PLAS1350 
BLAS1355 
3LAS1360 
BLAS1 380 
9LAS1390 
BLAS1400 
5LAS1410 
BLAS1420 
0LAS143O 
SLAS1440 
5LAS1480 
3LAS1h90 
BLASlsOO 
5LAS1510 
*LAS1520 
6LAS1530 
•3LAS1540 
3LAS1550 
HLAS1560 
6LAS1580 
SLAS1P90 
BLAS1600 
BLAS1610 
3L*S1620 
BLAS1630 
OLAS1640 
3LAS1650 
SLAS1660 
AoL AS1670 
BLAS1680 
5LAS1690 


^  ^  U  U  u  O  U  u  u  u 


*  >35**1  /Ni > ►,* w.a yv>  *:*"  a  **  v--"*'-'r*vc-*^ >^/yt v» <-*£> sWvThm^^*^ 3V*  !/uk*-'i“ "w ’A^oYs*,:  1  ;y V!^-^ p>tJ>  -  <>*??%  V\ *.* -‘r- >  a  <iSws*7*f*=* 
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c 

»«A!N!T 
PRINT  1690 

PSFCH*PCHAM*144.  »TRCH*(  TCnAM«-273 . 16) « 1 .8  SPSFAMe*PAMa*144. 

C  POUND  MOLES  CF  AIR. 

A I RPOL*PSFCH*V IN IT/(1545.*TRCH) 

IFINV.LE.OJ  r. o  TO  1790 

READ  1650.<PV(N).TV(S) .VV(N J.AVIM .PAV(N) .TAV(N) .NGPTVCN) .,'I*1.NV) 
PRINT1660«(PV(N) .TV IN) «VV(N) «AV(N) .PAV(N) .TAV(N) .NOPTV(N)  .N*1,NV) 
1790  CALL  GASES (V!\ IT .AIRMCL.OVPC.GC.Tfi) 

GA5LB*32.*M02+  2  3.  2*  28.*>.CC*  44.*MC02*  18.*MH20*  2.*MH2 

GASM0L«M02+,t:.2+.:C0*;c  32VH20MH2 
C  PI Nl T* INI T IAL  C:iAM_,;.R  PRCSSURr.  (PSFA)  AFTER  EXPLOSION. 

PINIT«0VPC'144.  +PSFCH 

PO*P INI T  SV9*V  INI  T  STE.MPO*TR  SNSAVE*i 
IFINV.LE.*')  oo  TO  2110 
I F { NOP  TV ( 1 ) «  NE  » 1 )  GO  TO  2110 


BREAK  WALLS  IF  PO  EXCEEDS  TABULATED  VALUES. 

DO  2070  N* 1 »NV  SNSAVE«N 
I F ( NOP T V ( N ) «  NE . 1 )  GO  TO  2110 
IF(P0.LT.PV(NJ*144. )  GO  TO  2070 
AsA+AVtNl 

CALL  MIX(PC.VO.TE‘«PO.GO.GASMOL,  P2 .V2  .TEMP2.G2 ,AI RADD) 
P0«P2  SV0*V2  STE.MP0«TEMP2  SG0*G2 
PO.  VO.  TEMPO.  GC  ARE  NO.v  AFTER  NEW  VOLUME  IS  ADDED. 

ADD  NEV/  AIR  TO  MCLES  OF  N2  AND  02. 

VN2-MN2*.  7<>0>*A  I  RADD  SM02  *MC2*  .  2095*  A  I  RADD 

GASM0L*M02*MN2*MC0*MC02*Mh20*Mh2 

POPSt*PO/144. 


PRINT  20!rtwiN),V2,A,P('PSI.G2 
I  FORMAT C*OFAILURF  LEVEL  IN  TAELE  EXCEcDFD.*/ 
1*  VOLUME  INCREASE (CJ  FT»**G12.4/ 

2*  NEW  TOT  VJL  <CU  FT)  **012.4/ 

3*  NEW  TOT  AREA  (SO  FT)  «*G12.4/ 

4*  NEtv  PRESSURE  <  PS  I A )  «*G12.4/ 


5*  NEW  GAMMA  **G12.4) 

070  CONTINUE 

INITIAL  BREAKING  INTO  NEW  CHAMPERS  IS  NOW  COMPLETED. 


NO  VENTING  IF  t REA *0. 

110  IF(A.EO.O.)  00  TO  S70 

PO* INITIAL  PEAK  PRESSURE  (PSFA). 

V0*INITIAL  VCLUKfc  (Co  FT). 

TO* I  NIT IAL  TEMP ( R  J . 

PRINT  216C 
2160  FORMAT  (*r>ncr,t, v  VENTING  CF  GASES*/ 

1*  OVERPR(PSI)  T  I‘-’E  ( SEC )  GASES(LB)  TEMP(R)  GAMMA  NEON*) 
OVPS I  *  (  P9-P5EAMU  )  /  1 44.  ST  I-VE 1  *0.  SG*GO 

GASLR*32.*M02+  28.».VN2*  28.*MCC*  44.*.v.CC2+  I8.*MH2Q*  2.*MH2 
PRINT  3030.  OVPS  I  .TIi'-'El.G  ASLli.  TEMPO.  G 
PA*PSFAMn  SGRAV*32. 2 
C  DENSITY  (Li/CU  FT). 

D0«GASL9/V0 

C  CRITICAL  PRESSURE  (PSFA). 


eLAS1700 

BLAS1710 

BLAS1720 

BLAS1730 

BLAS1740 

BLAG1750 

BLAS1760 

oLAS:.770 

BLAS1780 

BLAS1790 

BLAS1800 

8LAS1610 

BLAS1820 

BLAS1830 

BLAS1840 

BLAS1850 

BLAS1860 

6LAS1870 

BLAS1680 

BLAS1S90 

BLAS1900 

SLAS1910 

BLAS1920 

6LAS1930 

BLAS1940 

BLAS1950 

BLAS1960 

BLAS1970 

8LAS1980 

5LAS1990 

BLAS2000 

3LAS2010 

BLAS2020 

BLAS2C30 

BLAS2040 

BLAS2050 

BLAS2060 

6LAS2070 

BLAS2080 

BLAS2090 

BLAS2100 

BoAS21 10 

BLAS2120 

BLAS2130 

BLAS2140 

5LAS2 1 50 

BLAS2160 

BLAS2170 

BLAS2180 

BLAS219C 

BLAS2200 

BLAS2210 

BLAS2220 

BLAS2230 

BLAS2240 


PCaPA/t  t?./rr*i.n**«6/(f.“i.n  i  " 

C  PRESSURE  INCREMENT 
DP1»IP0-PA1*.?1 

N£0N*1  $TIVE1*0.  SKOUNT* 1 

Vl*VO 
C 

N«NSAVF 

C  PP6IN  VCNT I Nfi  LOOP. 

DO  3110  «N«1.1*00  SP1*P0-DP1  SG*GO 
IrlPl.GT.PA)  GO  TO  2360 
KOUNT *99  $P1«P0  SCO  TO  3010 

2360  GO  TO  <2380.24301 .NEON 
(  VENTING  FOR  Pl.GT.PC 
2380  IFIP1.LE.PC.AND.P0.GT.PC)  P1*PC 

otimei*(PO-p:i/pi«*i(a.*g-i.i/c2.»gi)*vi/a/ 

1  S0RT(GRAV*G**3»P0»*( l./Gi/DO* <2./<G*l.1 t«*t <G*1. 1/16-1. 11 1 
GO  TO  2470 

C  VENTING  FOR  Pl.LT.PC 
2430  IFIP1.LT.PA.AN0.P0.GT.PAI  P1*PA 
GG*(G-1.I/G 

DT  JME1* (PC-P1 J / (Pl**GG* (P1**GG“PA**GG J •*.5 1* VI /A/ 

1  SCRT(GRAV*G**3*2./( j-l.I*(PC»PA**2/0O»*G>**(l./G»I 
2470  TIME1*TIMF1.2TTME1 

IF  I NV.EQ.O I  00  TO  2620 
IFtfi.GT.NVl  GO  TO  2920 
IF(  NOPTVl.N  I .lT .2  1  GC  TO  2920 
C  CHECK  TIME  AGAINST  VENTING  TAeLE. 

IF(TIMEi.LT.TVIN) I  GO  TO  2920 
IF ( NOPTVIN 1  .£0.2 1  GO  TO  2560 
IFIP1.LT.PV(M*144.J  GO  TO  2920 
C  ADJUST  TI'-’El  TO  EQUAL  TV. 

2560  DT  t MEi *TV(  Nl-*(  T tMEl-DTIVL'l ) 

TIMF1*TV(NI 

GO  TO  I  2590. 26201  . NEON 

2590  DPI  »P1**( < 3.*G-1 • 1 / <2.*G1 1*A/V1*DTIME1* 

1  SORT I GRAV*G**3*I!,0**I i./GI /OO  J*( 2./(G+l.))**l(G+l» >/<  G-l. I ) J 
GO  TO  2650 
2620  GG*(G-1.)/G 

CPl*Pl**GC»1?l**GG-PA«*GG)*».5*A/Vl«0riMCl» 

1  SORT (gRav*g**?* ( 2. / c G~1 . )  i * i Po  »P4**->/oo**r,)**ci,/r,>  i 
2650  Pl*PO-DPl  30VP1*P1-PA  10VPSI 1*0VP1/144. 

C  REDUCE  MASSES  CUE  TO  VENTING. 

01*00*(Pl/P0l**Il./G) 

TEMPl'T£vpO»(Pl/PO) ** ( <  G-l. )/G) 

C  FRACTION  LEFT  AFTER  THIS  VENTING  STtP. 

FLEFT*01/D0 

V02«MG?*FLFFT  $M:i2»MN2*FLEFT  3VCO*MCO*FLEFT 
MC02*MC02*FLEFT  SMm20*VH20*FLEFT  r-'H2»KH2*FLEFT 
GASLB*D1»V1 

GASM0L*.v02*"f.2^.vC0+v,c02*.MH20+MH2 

PRINT  3030.0VPS11 .  THE!  .GAM.6.TEKP1  .G. NEON 

PRINT  2770.M.TIME1 

2770  FORi'.AT ( *  TIKE  HAG  REACHED  TV(»12*J**G12.4) 

C  BrGIN  VOLU'T-AREA  CHANGE  SECTION. 
a.a+AVINJ 


9LAS2250 
9LAS2260 
3LAS2270 
BLAS2280 
BLAS2290 
6LAS2300 
PLAS2310 
PLAS2320 
6LAS2330 
BLAS2340 
3LAS2350 
BLAS2360 
6LAS2370 
3LAS2380 
6LAS2390 
9LAS2400 
BLAS2410 
BLAS2420 
BLAS243C 
3LAS2440 
SLAS2450 
BLAS2460 
BLAS2470 
0LAS240O 
3LAS249C 
BLAS2500 
BLAS2510 
BLAS2520 
BLAS2530 
BLAS2540 
BLAS2550 
3LAS2560 
•U.AS2570 
BLAS25S0 
5LAS2590 
3LAS2600 
6LAS2610 
6LAS2620 
9LAS263C 
PLAS2640 
BLAS2660 
3LAS2660 
BLAS2670 
BLAS2680 
BLAS7690 
BLAS2700 
F.LAS2  71C 
3LAS2720 
BLAS2730 
5LAS2740 
BLAS2750 
SLAS2760 
BLAS2770 
BuAS2  78C 
9LAS2790 


I 


?  *?*•'  '^'i>«,4;^'^-*"  4>**|27 
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t\ 


CALL  MIX(P1.V1.TEMP1.G.GASM0L.  P2.V2.TEMP2.G2.AIRADD) 

PI  »P2  1V1»V2  <TEMP1«TEis*P2  *6*61«G2 

VN2*.vN2*.7<505*A  TRADD 
M02«MC2*.2095*AIRADD 

GASLB*32.*M02*  28.*KN2+  28.«MCO.  44.*MC02«-  1B.*MH20*  2.*MH2 
V1«V?  J°1 °St  *P! /1 44. 

Oi-GASLO/VI  «''P1«{01-PA)*."1 
PR  TNT  2010»VV(N),Vl.A.PlPSItGl 
N«N+1  SKOUNT »<M*X2  SGO  TO  2970 

C  END  VOLUME-AREA  CHANGfc  SECTION. 

C 

C  REDUCE  MASSES  DUE  TO  VENTING. 

2920  DlO0*(Pl/P0)*«U.AU 

TEMPI “TEMPO* < PI /Pfl) **( (G-l. )/G) 

FLEFT«D1/D0  SGASLB“01*Vl 

M02«M02*FLEFT  SMN2«.*'N2*  FLEET  SKCO»MCC*FLEFT 

MC02«MC02»FLEFT  JMH20«MH2C*FLEFT  SKH2«MH2*FLEFT 
2970  IF (PI. EC. PC  1  KCUNT«KMAX2 
IF ( Pl.EO.PA J  K0UNT*KMAX2 
C  PRINT  EVERY  KMAX2-TH  LINE. 

IF (KOUNT .NE.KMAX2 I  GO  TO  3050  SK0UNT«0 
3010  0VPSI1*<P1-PAJ/IA4. 

PRINT  3030 .OVPS II . T I ME 1 *G ASLB . TEMP 1 »G .NEON 
3030  FORMAT ( 1H  .4G12.4.F7.4, 14) 

IF ( KOUNT. E0.901  GO  TO  3130 
3050  P0«P1  JKCUNT“K0UNT*1  SD0«D1  *TEMPO«TEMPl 
IF (Pl.LE.PA )  GO  TO  3130 
CALL  GA.05A (TEMPI. GO »OUMJ  SG»G0 
C  RECALCULATE  PC  WITH  NEW  G. 

IF(Pl.GT.PC)  PC«PA/((2./(G*l.) )**(G/(G-1.) )) 

IF (Pl.LE.PC)  NF1N«2 

;iir  continue 

C  END  VENTING  LOOP. 

3130  IF ( NR.GT.O )  GO  TO  730 
GO  TO  570 

C  END  VENTING  SECTION. 

END 


BLAS2800 

8LAS2610 

BLAS2820 

BLAS2630 

BLAS2640 

BLAS2850 

BLAS2860 

BLAS2870 

BLAS2680 

BLAS2690 

5LAS2900 

BLAS2910 

BLAS2920 

BLAS2930 

BLAS2940 

BLAS2950 

8LAS2960 

BLAS2970 

BLAS2580 

BLAS2590 

BLAS3000 

BLAS3010 

BLAS2020 

BLAS3030 

BLAS3040 

BLAS3050 

BLAS3060 

BLAS3070 

BLAS3080 

BLAS3090 

BLAS3100 

BLAS3U0 

8LAS3120 

BLAS3130 

BLAS3140 

BLAS3130 

BLAS3160 


I 
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SUBROUTINE  FIX<P3.VO. TEMPO. GO.GASKCL.  P2.V2.TEMP2.G2.AIRADDI  MIX 

MIX  THE  GASES  I N  ADJACENT  CHAMBERS.  MIX 

COMMON/VENT/  PV< 1CI .TV< 1CJ  tVVIiO ) »AV< 10 • .PAV< 101 »TAV< 10 >  »NOPTV< 10IMIX 
1»N  MIX 

V2«V0n'VtNJ  MIX 

MOLES  OF  AIR  IN  THE  ME*  VOLUME.  MIX 

AIRADD*  <  PAV  <N|  *144.)*VV<N)/<1 545. *<  TAV<N)42?3.16I *1 »B)  MIX 

ITERATE  TO  FIND  NEW  T.P.G.  MIX 

TOTMOL*GASf-'OLfAIRADO  MIX 

DIF1«1.£10  SDT£MP*TEKP0/1C0.  MIX 

PO,  VO.  TEMPO.  GO  are  BEFORE  NEW  VOLUME  IS  ADDED.  MIX 

P2«  V2.  TEMP2.  G2  ARE  AFTER  NEa  VOLUME  IS  ADDED.  MIX 

FIRST  GUESS  FOR  TEMP2.  MIX 

TEMP2«TEMP0  MIX 

00  240  J-1.100  MIX 

P2A«<T0TM0L*1545./V2I*TEMP2  MIX 

CALL  GAMMAITEMP2.G2.0UMI  MIX 

P2B«<G2-1.I/V2  *(<1.4-1.»*P0*V0  •MG0-1.I*PAVIN)»VVIN))/  MIX 

1  <<1.4-1. I*<G0-1. II  MIX 

0IF2«ALS<P2A-P2B)  MIX 

IF<DIF2.GT.0IFll  r0  TO  250  MIX 

CONTINUE  SEARCH  FOR  CORRECT  TEMP2.  MIX 

TEMP2«TEVP2-DTEVP  SDIF1«0IF2  MIX 

240  CONTINUE  MIX 

250  P2«<P?A+P?9l/2.  MIX 

RETURN  MIX 

END  MIX 
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SUBROUTINE  HEOATA 
C  TABLES  OF  EXPLOSIVES  DATA. 

COMMON/DATAl/WLB, NUMBER ,RLOD,CASE,V INI T,AIN I T .PAMB.TAMB .ALTKFT , 
1  PCHAM »T CHAM, NOP T  «NV*NR*  WFACT  ,£FORM 
COMMCN/WTFRAC/WFC  »WFH*WFN»WFO  *WFA 
COMMON/HE/NUMHE 1 9  >  »HEFRAC {  9  J 

DIMENSION  EQWT ( 40  1  »EF { 40 ) »FC( 40 1 »FH 1 40 1 »FN ( 40 1 *FO ( 40 1 »FA(40 1 
DIMENSION  NAME (40*5) *NAMES< 51 
C  EOWT ■EQUIVALENT  WEIGHT  REFERRED  TO  TNT, 


EF«ENERGY  OF  FORMATION 


DATA! (NAME!  1 

.1), 

I* 

110H 

» 10H 

, 

DATA  EOWT (  1) 

,EF( 

1 

1  /l.OO, 

-78. 

40 

DATA! (NAME!  2 

•  I ) . 

I* 

110H 

,  10H 

, 

DATA  EQWTt  21 

,EF  ( 

2 

1  /I . 13, 

-307 

.1 

DATA! (NAME!  3 

,1), 

I* 

110H 

,  10H 

, 

DATA  EOWT (  3) 

,EF( 

3 

1  /0.85, 

-386 

.3 

DATA! 'NAME!  4 

,n. 

I« 

HOHi 

,  10H 

* 

DATA  EOWT (  4) 

,EF( 

4 

1  /l. 17, 

-242 

•  8 

DATA! (NAME!  5 

.1). 

I* 

110H52/48 ) 

,  10H 

, 

DATA  EOWT (  5) 

,EF( 

5 

1  /0.90. 

-238 

.5 

DATA ((NAME I  6 

,n«t* 

110H 

,10H 

, 

DATA  EQWT (  6) 

,EF( 

6 

1  /l. 14, 

22. 

79 

DATA( <  NAME  t  7 

,n. 

I® 

110H39. 6/1.0) 

,  10H 

, 

DATA  EQWT (  7) 

,cFl 

7 

1  /1.10. 

4. 

33 

DATA ((NAME!  8 

.11, 

I* 

110H 

*  10H 

, 

DATA  EOWT (  8) 

»EF  ( 

8 

1  /1.19, 

57. 

00 

DATA! (NAME!  9 

,1), 

1  = 

110H 

,  10H 

, 

DATA  EQWT (  9) 

,EF  ( 

9 

1  -  /1.09, 

24. 

93 

DATA! (NAME! 10 

,1), 

1  = 

110H 

,  10H 

, 

DATA  EOWT(IO) 

, EF ( 10 

1  /1.23. 

-276 

.4 

DATA( (NAMEtll 

,1), 

1  = 

110H 

,  10H 

, 

DATA  EOWT ( 11 ) 

,EF ( 11 

1  /la  18, 

-239 

.5 

HEDA001C 
HEDA0020 
HEDA0030 
HEDA0040 
HEDA0060 
HEDA0065 
HEDAOOTO 
HEDA0075 
KEOA0080 

CAL/GJ.  HEDAC090 

*5 1 *10H 1 3 7H  TNT  »10H  »\0H  »HED».009S 

3  »F6»  3 1 )  HEDA0096 

*FC (  1>,FH(  1»*FN|  1 ) ,F0(  1),FA(  li  HEDA0100 

» 370 *  .022  *  .185.  .423.  .000/  HEDA0110 

.5 (»10H( 37H  TNETB.10H  .10H  .HEDA0115 

3.F6.3) ]  HEDA0116 

»FC I  2),FH<  2) »FN(  2),F0(  2)»FA(  2>  HEDA0120 

.186.  .017.  .217.  ,580.  .000/  HEDA0130 

.5 ) «10H( 37H  LXPLO.lOHSi VE  D  »10H  .HECA0135 

3  «F6. 3 ) t  HEDA0136 

.FC (  3I.FHI  3 ) »FN(  3),F0I  3I.FAI  3)  HEDA0140 

.293.  .025.  .227,  .455,  .000/  HEDA0150 

,5)«10H<37H  PENT0.10HLITE  {PETN.10H/TNT.50/50.HEDA0155 

3.F6.3M  HEDA0156 

.FC(  4I.FHI  4)  ,FN(  4),F0t  4I.I/.I  4)  HEDA0160 

.280,  .024,  .182,  .514,  ,000/  HEQAOl?0 

,51 *10H{ 37H  PlCRA.iCHTOL  IEXPL0.10HSIVE  D/TNT .HEDA0175 


3.F6.3J » 

»FC (  5>.FH(  5 )  ,FN(  5) ,F0(  5),FA(  5) 
.329,  .024,  .207,  .440,  .000/ 


MEDA0176 

HEDA0180 

HEDA0190 


*5 1 *10H( 37H  CYCLO.IOHTOL  (RDX/T.10HNT, 70/301  .NEDA0195 


3 ,F6.3 ) ) 

,FC(  6I.FHI  6I.FNI  6I.F0I  6>,FA(  61 
.225,  .026,  .320,  .429.  .000/ 


HEDA0196 

HEDA0200 

HEDA0210 


,5 ) *10H( 37H  COMP  ,10HB  (RPX/TNT.10H/WAX,59.4/»HEDA0215 
3 ,F6, 3 ) )  HEDA0216 

,FC (  7) ,FHI  7),FNI  7) ,1  X  7),FA(  7J  HEDA0220 

.252,  .026,  .298,  .424,  .000/  HEDA0230 

,5)*10H(37H  RDX/W.10HAX,  98/2  ,10H  ,HEDA0235 

3 , F6.3 ) )  HEDA0236 

»F C<  8I.FHI  8 ) vFN (  8).F0(  B),FA(  81  HEDA0240 

.176,  .030,  .371,  .423,  .000/  KEDA0250 

,5 ) « 10H< 37H  COMP  ,10HA-3  (RDX/W.10HAX.91/9I  .HEDA0255 

3.F6.3JI  HEPA0256 

,FC (  91.FHI  9I.FNI  9>,F0(  9I.FAI  9)  HEDA0260 

.225,  .038,  .344,  .393,  .000/  HEDA0270 

, 5 ) * 10H( 37H  TNETB.10H/AL.  90/10, 10H  .HEDA0275 

3.F6.3I)  HEDA0276 

.FCUOI.FHUOJ.FNIJOI.FOUO)  »FA(10.  HEDA0280 

.168,  .014,  .196,  .522,  .100/  HEDA0290 

, 5 ) * 10HI 37H  TNETB.10H/AL,  78/22, iOM  .HEDA0295 

3.F6.3) I  HEDA0296 

,FC(ll),FH(ll),FN(ll),FO(ll),FA(ll)  HEDA0300 

.146,  .012,  .170,  .452,  .220/  HEDA0310 
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lDjTAUNAMEn2:n,,:l.^o„,,7H  TNETB.lOH/AL,  72/28.10H 

DATA  EOWT(12».EF(12),FCU2!!FHU2),FNa2),FOa2>,FA(12> 

oATAMNiS^.rff  t1:  :13?*  *011*  •l57‘  *4ie*  -«o/ 

w  ‘  :KS!!"  TNETB«10m/al»  65/35. 10H 

OATA  EQWT<13>,EFU3>,FCU3J,FHU3>,FNU3>,FOtl3>»FAll3> 

^DATa^inamf#  :12J*  *010*  ,U2’  *377*  .330/ 

i0ATAnNAHE(U.n.I-U5,.lOH«37H  TRJ  TO.IOHNAL  < TNT/A.10HL.80/20) 

OATA  e0WT«14>,EF(14»,FCCl4UFH«lA),FNtUj,F0U4|,FA<I4| 

DATAnN?MF,1«i‘'??,?2.,  22<,f*  *018*  •I48‘  ‘338,  *200/ 

^AtA.IMHE, 15.1).,. 1.5,., 0H.37H  «0X/A,1<,HL/»AA.  M/.MH10/2 

OATA  EOWT(15).Fr|(j5l,FCn5J.FH(13»,FNa5),FOa5l.FAa5» 

J  /1.30.  50.38,  *160*  .027.  .333,  »3B0«  .iaa/ 

OA7A<(NAMEn6,r),f-j,5,,loH(37H  ROX/A.lOHl/WAX,  78/.10H20/2 
-  fl0H  »13»F6.3)| 

OATA  EOWT(16),EF(16),FC(16),FHa6),FNJ16).FOn6).FAt16» 

OATAMNAMFM7  *°24»  *295*  *337*  ‘*00/ 

i0ATA»«NAME(l7,n,I.l,3,.i0H(37H  RDX/A.10HL/WAX.  74/.10H21/5 

OA7A  EOWT(17>,EF(17>,FCa7l,FH(17).FNa7J,FO(17).FA(17) 

^DATAnNAMF. i(  f?*?6:  :163*  •°27*  •280»  «320i  .2iO/ 

1  ini.  *^AME( I8»II»I«l»5)«10HO7H  RDX/A.lONL/WAX,  74/.10H22/4 
....  -  ,10H  *  13 *F6*3>  > 

DATA  EQWT( 16  >  <EF ( 18| ,FC { 18) ,FH(18) ,FN{ 18) »FO( 18) «FA( 18) 

nSJA‘  ‘NAME«  19.J  >  !  1.  j,5)«lOH<;7H6RDX/A!liHE/jAX:  6^%H33/5 

nS^nSf  '<"'*•*'>'«  ■■  («ox„ioh,n„al..3, 

iMI./!ow,:2,,.ep,|o;.:c,|owph,2,u«,„weo,|„,.:a,!o> 

mi%!«"E,2‘:!w';!;^:;os;;’7H  h-‘  <'1°h»°x/7»t/al.,0h/»ax,«,» 

iBArA/M,7.22,.|r,2!,.rc,«,.f„,«„™,„„„,„„F,,„1 
data  EOWTt23),EF(23),FC(23),FH(23J,FN(23),FO(23|,FA(23) 

DATAlJliSp#**"?  !  :2°°*  ‘°22’  *17l»  *257’  .380/ 

110H  ^  jloH^ ‘"ijUSlSlI3™  TNETB'l°H'Rt>X/AL.  3.10H9/26/35 

DATA  EOWT(24).EF(24),FC(24),FH(24),FN(24),FO(24),FA«24) 

X  /1.24,  -102.6.  .115*  ,013.  ,184,  ,338.  .350/ 


,NAME‘25*1 > » l*1.5)«10H<3?H  ALUHI .10HN0M 
110H  »10H  .13.F6.3}) 


»10H 


.HEDA0315 

HEDA0316 

HEDA0320 

HEDA0330 

.HEDA0335 

HEDA0336 

HEDA0340 

HEDA0350 

.HEDA0355 

HE0A0356 

HEDA0360 

HEDA0370 

.HE0A0375 

HEDA0376 

HEDA0380 

HEDA0390 

.HEDA0395 

HEDA0396 

HEDA0400 

HEDA0410 

.HEDA0415 

HEDA0416 

HEDA0420 

HEDA0430 

.HEDA0435 

HEDA0436 

HEDA0440 

HEDA0450 

.HEDA0455 

HEDA0456 

HEDA0460 

HEDA0470 

.HEDA0475 

HEDA0476 

HEDA0480 

HEDA0490 

.HEDA0495 

HEDA0496 

HEDA0500 

HEDA0510 

.HEDA0516 

HEDA0516 

HEDA0520 

HE0A0530 

.HE0A0535 

HEDA0536 

HEDA0540 

HEDA0550 

.HEDA0555 

HEDA0556 

HEDA0560 

HEDA0570 

.HEDA0575 

HEDA0576 


D-10 


teliliilitfTiiii'ltlh'ftiifil 


vAwniritttt 


DATA  EQWT (25) 
1  /  O.t 

DATA  1 1 NAME (26 
110H 

DATA  EQWT (26) 
1  /  O.t 

DATA! (NAMEI27 
110H 

DATA  EQWT (27) 
1  t  O.  * 

DATA( (NAME ( 2U 
1 10H 

DATA  EQWT ( 28 ) 
1  /  0.. 

DATA( (NAME(29 
110H 

DATA  EQWT (29) 
1  /  O.t 

DATA( (NAME(30 
110H 

DATA  EQWT ( 30) 
1  t  0.* 

DATA) (NAME! 31 
X10H 

DATA  EQWT(31) 
1  /  O.t 

DATA) (NAME) 32 
110H 

DATA  EQWT (32) 
1  /  0.. 

DATAt (NAME) 33 
110H 

DATA  EQWT) 33) 

1  /  O.t 

DATA( (NAME) 34 
110H 

DATA  EQWT (34) 
1  /  0.. 

DATA) (NAME) 35 
110H 

DATA  EQWT (35) 

1  /  -o.. 

DATA( (NAME) 36 
UOH  , 

DATA  EQWTI36) 
1  /  O.t 

DATAt (NAME(  37 
UOH 

DATA  EQWT (37) 
1  /  O.t 

DATAt (NAME(38 
11GH 

DATA  EQWT (38) 
1  /  O.t 


♦EFI25 

0. 

.!).!« 
tlOH  t 
tEF ( 26 
-392. 
tilt!* 
tlOH  t 
.EF (27 
66.16 
•  I ) 1 1* 
tlOH  t 
tEF( 28 
-<*07.1 
«l)tl* 
tlOH  t 
tEF( 29 
16.26 
<().!> 
tlOH  t 
.EF ( 30 
0. 

tilt  I* 
tlOH  t 
tEF ( 31 

0. 

.(>tl* 
tlOH  t 
tEF ( 32 
0. 

i.!)tl» 
tlOH  t 
tEF ( 33 

0. 

tl)tl« 
tlOH  t 
tEF (34 
0. 

't  I  )  »  I* 
tlOH  t 
.EF ( 35 
0. 

>tl)tl« 

tlOH  t 
.EF ( 36 
0. 

tl)tl« 
tlOH  t 
tEF ( 37 

0. 

ltl)tl« 
tlOH  t 
tEF ( 38 
0. 


NOLTR  72-231 


»FC(25).FH(25)»FN(25).FO(25)tFA(25) 
O.t  Ott  O.t  O.t  1./ 

.5)»10H(37H  WAX  tlOH  tlOH 

3tF6.3) ) 

tFC(26> .FH(26).FN(26).F0(26)tFA(26) 
.856.  .144.  O.t  O.t  0*/ 

.5 ) «10H( 37H  RDX  tlOH  tlOH 

3tF6.3) ) 

tFC(27)tFH(27).FN(27).F0(27)tFA(27) 
.162.  .027.  .379.  .432.  0./ 

t 5 ) «1QH( 37H  PETN  .10H  tlOH 

3.F6.3)) 

»FC(28)tFH(28).FN(28).FO(28)tFA(28) 
.190.  .026.  .177.  .607.  0*/ 

.5 ) »10H( 3 7H  TETRY.10HL  tlOH 

3.F6.3) ) 

»FC( 29 ) .FH(29 I .FN(29) »F0( 29) »FA(29) 
.293.  .017.  .244.  .446.  0./ 

. 5 ) ■ 10H  ,10H  tlOH 

3.F6.3) ) 

tFC(30) .FHI30) .FN(30).FO)30)tFA(30) 

O.t  O.t  O.t  O.t  0*/ 

. 5 ) *10H  tlOH  tlOH 

3.F6.3)) 

:FC(3U.FH(31),FN(31)tFO(31).FA(3U 
O.t  O.t  O.t  O.t  0./ 

.5 ) »10H  tlOH  tlOH 

3.F6.3H 

.FC(32).FH(32) .FN(32).FO(32).FAl32) 
O.t  O.t  O.t  O.t  0./ 

.  .5 ) *10H  tlOH  tlOH 

3.F6.3)) 

tFC(33).FH)33) t FN ( 33 ) .FO ( 33 5 . FA( 33 ) 
O.t  O.t  O.t  O.t  0./ 

,.5)«10H  tlOH  tlOH 

3.F6.3)  ) 

.FC ( 34) »FH ( 34 ) ,FN( 34 ) ,F0( 34) ,FA( 34) 
O.t  O.t  O.t  O.t  0./ 

.  .5 ) *10H  tlOH  tlOH 

3.F6.3)  ) 

.FC(35).FH(35).FN(35).FO(35)»FA(35) 
O.t  O.t  O.t  O.t  0./ 

,  .5 ) *10H  tlOH  tlOH 

3tF6.3) ) 

»FC(36)tFH(36)tFN(36)tFO(36)tFA(36) 
O.t  O.t  O.t  O.t  0./ 

.  .5 1 ®10H  tlOH  tlOH 

3.F6.3) ) 

tFC(37) tFH(37).FN(37)tFO(37)tFA(37> 

<  O.t  O.t  0»t  O.t  0./ 

.  1 5 ) * 10H  tlOH  tlOH 

3tF6.3)  > 

.FC(38)tFHi38) tFN(38) tFO(38).FA(38) 
O.t  O.t  O.t  O.t  0./ 


HE0A0580 
HEDA0590 
.HEDA0595 
HEDA0596 
HEDA0600 
HEDA0610 
tHEDA0615 
HEDA0616 
HE0A0620 
KEDA0630 
tHEDA0635 
HE0A0636 
HEDA0640 
HE0A0650 
tHEOA0655 
HE0A0656 
HEDA0660 
HEDA0670 
tHEDA0675 
HEDA0676 
HEDA0680 
HEDA0690 
tHE0A0695 
HEDA0696 
HEDA0700 
HEDA0710 
tHEDA0715 
HEDA0716 
HEDA0720 
HEDA0730 
tHEDA0735 
HEDA0736 
HEDA0740 
HEDA0750 
tHEDA0755 
HEDA0756 
HE0A0760 
HEDA0770 
tHEDA0775 
HE0A0776 
HEDA0780 
HEDA0790 
tH6DA0795 
HEDA0796 
HE0A0800 
HEDA081C 
tHEDAOSl 5 
HEDA0816 
HEDA082C 
HEDA0830 
tHE0A0835 
HE0A0836 
HEDA0840 
HE0A0850 


r\r% 
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DATA( (NAME) 39* I) • l«l *5 ) *10H  »10M  »10H 

110H  *  10H  .I3.F6.3M 

OATA  EOWT ( 39) .EFI 39) *FC(39) »FH<  391 »FN( 39) »FOt  39) »FA( 39) 

1  /  0*.  0*.  0*.  0*.  0*.  0*.  0./ 

DATA! (NAME (40. I).I*1.5I»10H  *10H  »10H 

X10H  .10H  .I3.F6.3)) 

OATA  EQWmo>.EF<40».FC<40>.FH<40)«FNt40).FO(*0)»FA{40> 

1  /  0*.  0*.  0*.  0..  0..  0..  0./ 

c 

IF (NUMBER. EO.-l )  CO  TO  1200 

N-NUMBER  SWFACT«EOWT(N)  SEFORM«EF(N)/1000. 

C  WEIGHT  FRACTIONS 

WFC*FC(N)  *WFH«FH(N)  SWFN«FN(N)  SVF0*F0(N>  SWi“A«FACN) 
DO  1030  L«1.5 

1030  NAMES ( D -NAME (N.L I  SPRINT  NAMES 

RETURN 

MIX  UP  AN  EXPLOSIVE  FROM  COMPONENTS  IN  LIST. 

1200  EFORM»WFC«WFH«WFN»WFO«WFA»0.  SPRINT  1370 
00  1290  1*1 .9  SN*NUMHE( I )  SHF-HEFRACI 1 1 
IF(N.EQ.O)  GO  TO  1300 
DO  1225  L«1.3 

1225  NAMES! LI *NAML (N.L)  SPRINT  NAMES. N.HF 

EFORM*EFORM+HF*EF (N)  SWFC«WFC+HF*FC(N) 

WFH*WFH4HF*FH<  N )  SWFN«WFN4HF*FN ( N I 

WF0«WF04HF*F0 ( N )  SWFA«WFA4HF*F A I N I 

1290  CONTINUE 

1300  EFORM.EFORM/IOOO.  SRETURN 

1370  FORMAT(*0MAK£  UP  SPECIAL  HE  MIXTURE—*/ 

1*  NAME*  »29X»*NUMBER  WT  FRAC*I 
END 


.HEDA0855 

HEDA0856 

HEOAO86O 

HEDA0870 

•HEOA0875 

KkOA0676 

HEOAO88O 

HEOA0890 

HEDA0990 

HEDA0995 

HEOAIOOO 

HEDA1010 

HEOA1020 

HEOA1025 

HEOA1030 

HEOA1050 

HEDAU90 

HE0A1195 

HE0A1200 

HE0A1210 

HEDA1220 

HE0A1222 

HEOA1225 

HE0A1230 

HEDA1240 

HE0A1250 

HE0A1290 

HEDA1300 

HE0A1370 

HE0A1375 

HEDA1400 
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SUBROUTINE  GAMMA! T.G.U)  GAMMOOlO 

COMMON /MOLGAS/M02 »MN2  .MC0.MC02 .MH20.MH2  GAMM0020 

REAL  M02  »MN2 .MC0.MC02  *MH20*MH2  GAMM0030 

C  M02  ETC  ARE  LB  MOLES  BUT  UNITS  CANCEL  HERE*  GAMM0040 

U2»5. 76+20. /T#*.54.000578*T  GAMM0060 

U4-11. 515-172. /T**.5+1530./T  GAMM0070 

U6«9.47-3470./T+l  1-60000./T++2  GAMM0080 

U8-19.86-597. /T##. 5+7500. /T  GAMM0090 

U9«9. 46-3290. /T41070000./T**2  GAMM0100 

Ul« 16. 2-6530. /T41410000./T**2  GAMM0110 

U«U2*MH2+U4*M02+U6*MN2+U8*MH20+U9*MC0+U1*MC02  GAMM0120 

U«U/ ( MH2+MO2+MN2+MH2O+MC0+MC02 )  GAMM01 30 

G«U/(U-1 .987 )  GAMM0140 

RETURN  GAMM0150 

END  GAMM0160 
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SUBROUTINE  GASE5I V.AIRMOL.P2.G.TR)  GAS 

C  INPUT  DATA  ARE  NPRINT.V.AIRMOL.WC.WH.WN.WO.WA.  GAS 

C  OUTPUT  DATA  ARE  P2.G.TR.Q.M02.MN2.MC0.MC02.MH20.MH2.  GAS 

C  V«CHAMBER  VOLUME  (CU  FT  I.  GAS 

C  AIRMOL»LB  MOLES  OF  AIR  IN  CHAMBER*  GAS 

C  P2«OVENPRES5URE  IN  CHAMBER  DUE  TO  HE  GASES  (PSD*  GAS 

C  G-SPECIFIC  HEAT  RATIO  OF  HE  GAS-AIR  MIXTURE  IN  CHAMBER*  GAS 

C  TR» TEMPERATURE  (liANKINE ) •  GAS 

C  STATIC  CHAMBER  PRESSURE  CALCULATION  ACCOUNTING  FOR  AVAILABLE  02*  GAS 
C  COMBUSTION  PRODUCT  SEQUENCE  IS  H20*  Al203.  CO.  C02.  GAS 

COMMON/OATAl/WLB. NUMBER. RLOO. CASE. VIN I T .AINI T .PAMB. TAMB. ALTKFT .  GAS 

1  PCHAM.TCHAM.NOPT.NV.NR,  WFACT.EFORM  GAS 

C0MM0N/M0LGAS/M02.MN2  *MC0*MC02 .MH20.MH2  GAS 

COMMON/HEMASS/WC.WH.WN.WO  «WA  GAS 

COMMON/GAS/N1.N2.N3.N4.N5.N6.N7.  GAS 

1  M1.M2.M3.M4.M5.M6.M7.M8.M9.  R.0.X2  GAS 

REAL  N1.N2.N3.N4.N5.N6.N7  GAS 

REAL  Ml .M2.M3 .M4.M5 .M6.M7 .M8.M9  GAS 

REAL  M02.MN2.MC0.MC02.MH20.MH2  GAS 

PRINT  210  GAS 

210  FORMAT  OOPROPERT I ES  OF  GASES— #1  GAS 

Q"R*M2«M5«M7*M8*M9*Ml»0*  GAS 

C  GRAM  MOLES  C.H2.N2.02.AL2  GAS 

Nl«WC*4S3.6/12.  GAS 

N2«WH*453.6/2.  GAS 

N3«WN*453.6/28.  GAS 

N4«W0*453.6/32.  GAS 

N5-WA*453. 6/53. 963  GAS 

N6-N3+.7905*AIRM0L*453.6  GAS 

N7«N4*.2095*AIRMOL*453.6  GAS 

C  R-GMOLES  02  LEFT.  GAS 

R"N7-N2/2.  SIFIR.GT.O.)  GO  TO  340  GAS 

M8*2.*IR+N2/2. I  SM2-N2-M8  $R«0«  SQ«57.80»M8  SGO  TO  450  GAS 

340  M8-N2  SQ>N2*57. 80  SR«R-1.5*N5  SIFIR.GT.O.)  60  TO  370  GAS 

M5*lR+i.5*N5)*l2./3. J  SO«Q-.M5*400.3  SR»0.  SGO  TO  430  GAS 

370  Q»Q+N5*400.3  SR.R-N1  SIFIR.GT.O.)  GO  TO  510  GAS 

R»R*Nl/2.  SIFIR.GT.O.)  GO  TO  420  GAS 

M9«2.»IR+Nl/2.)  $Q«Q+M9*26. 42  SR«0.  SGO  TO  470  GAS 

420  M1«2.*R  SM9«N1-2.*R  SQ«Q+M1*94.05+M9*26.42  SR>0.  SGO  TO  490GAS 

430  RESULT-M5/N5  SPRINT  440.RESULT  SGO  TO  540  GAS 

440  FORMAT!*  PERCENT  LAST  PROOUCT  IAL203)  •»  G12.5)  GAS 

450  RESULT*M2/N2  SPRINT  460. RESULT  SGO  TO  540  GAS 

460  FORMAT I*  PERCENT  LAST  PROOUCT  IH20)  •»  G12.9)  GAS 

470  RESULT-M9/N1  SPRINT  480. RESULT  SGO  TO  540  GAS 

480  FORMAT  I  *  PERCENT  LAST  PRODUCT  ICO)  »*  G12.5)  GAS 

490  RESULT-M1/N1  SPRINT  500 .RESULT  SGO  TO  540  GAS 

500  FORMAT  I*  PERCENT  LAST  PRODUCT  IC02)  ■*  G12.5)  GAS 

510  PRINT  520  GAS 

520  FORMAT I*  OXIDATION  COMPLETE*)  GAS 

Ml-Nl  S0«0+94.05*N1  GAS 

540  X2*R+M2+N6*M8+M9*M1  GAS 

X3«X2-A!RMOL*453.6  GAS 

Q«Q*.592*X3  GAS 

C  O-ENERGY  RELEASED  (KCAL).  GAS 
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Q*Q+WLB*453.6*EF0RM 
C  LB  MOLES  OF  GASES. 

M02*R/453.6  SMN2*N6/453.6  »MCO*M9/453,6  SMC02«Ml/453.6 

MH20-M8/453.6  SMH2-M2/453.6 

0T«100.  SQ1*0.  $U«7.  $T*1»8*( TCHAM+273.16 ) 

C  INTEGRATE  UNTIL  T  IS  FOUND  SO  ENERGY  EQUALS  Q. 

00  760  J*1.100  $T«T+DT 

CALL  GAMMA t  T.G.U)  S0Q«<U-1.987)*DT*X2*. 0005556 

Q1*Q1+DQ  SIF(Ql.GE.O)  GO  TO  780 
760  CONTINUE  SPRINT  770 
770  FORMAT ( *  T  HAS  REACHED  UPPER  LIMIT.*) 

C  CORRECT  T  SO  Q1  HITS  Q  EXACTLY. 

780  T*T-( Ql-Ql /DQ*DT  SG«U/{ U-1.987) 

C  ABSOLUTE  PRESSURE  (PSIA) 

P2-X2/453. 6*1. 987*778. /144.  *T/V  STR*T  SP2-P2-PCHAM 
TF«T-460.  SPRINT  810. TF 
810  FORMAT (*  TEMPERATURE.  DEGREES  F  «*G12.5) 

QPERG*Q/ I453.6*WLB)  SPRINT  814.QPERG 
814  FORMAT ( *  ENERGY  RELEASEIKCAL/G)  «*G12.5) 

PRINT  830. G 

830  FORMAT ( *  SPECIFIC  HEAT  RATIO  »*G12.5) 

PRINT  870. P2 

870  FORMAT (*  GAS  CVERPRESSURElPSI )  «*G12.5) 

RETURN 

END 


GAS  0580 
GAS  0600 
C-AS  0610 
GAS  0620 
GAS  0630 
GAS  0640 
GAS  0650 
GAS  0660 
GAS  0750 
GAS  0760 
GAS  0770 
GAS  0775 
GAS  0780 
GAS  0785 
GAS  0790 
GAS  0800 
GAS  0810 
GAS  0812 
GAS  0114 
GAS  0820 
GAS  0830 
GAS  0860 
GAS  0870 
GAS  0880 
GAS  0890 


;vs  <s , 
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SUBROUTINE  TNT(L) 

TNT 

0010 

c 

POSITIVE-PHASE  PROPERTIES  FOR  1  LB  TNT  IN  SEA- 

LEVEL 

AIR. 

TNT 

0020 

COMMON/TNT IN/R TNT. KMAX1 

TNT 

0030 

COMMON/TNTOUT /PTNT  »TSTNT 

.T'ATNT.TCTNT  , PC  TNT 

,PSI (40) *71(40 t ,T2(40) 

TNT 

0040 

1* 

JJ.XIMP1 

TNT 

0050 

DIMENSION  PI108I 

*R< 1081 , 

TS1108), 

TP ( 108 ) , TC ( 23 ) 

TNT 

0060 

c 

RADIAL  DISTANCE  FROM 

CHARGE 

CENTER  (CM). 

TNT 

0070 

OATA  R/  4.054, 

4.680, 

5.700, 

7.120, 

9.200, 

10.6, 

12.4, 

TNT 

0080 

I 

14.7,  17.6, 

19.0, 

21.7, 

25.0, 

28.1, 

30.0, 

32.0, 

TNT 

0090 

2 

34.5.  37.4, 

41.0, 

46.0, 

48.8  , 

52.4, 

57.3, 

65.  , 

TNT 

0100 

3 

66.4,  72.0, 

79.0, 

85.3, 

89.0, 

93.0, 

98.4, 

105.. 

TNT 

0110 

4 

114.,  125., 

133., 

142., 

154., 

171., 

180., 

196.. 

TNT 

0120 

5 

218.,  23°. , 

251., 

269., 

289., 

316., 

350., 

402., 

TNT 

0130 

6 

438.,  483., 

546., 

640., 

6C8.  , 

781., 

920., 

1070*. 

TNT 

0)40 

7 

1160.,  1260., 

1400., 

1580., 

1830., 

2200., 

2450*. 

2790*. 

TNT 

0150 

8 

3220.,  3900., 

4250., 

4960., 

6010>, 

7030., 

7710., 

8540*, 

TNT 

0160 

9 

9600.,l.l0E+*» 

1.29E+*, 

1.57E+4, 

1.77E+4.2. 

04E+4, 

2.40E+4, 

2.92E+4, 

TNT 

0170 

13 

.21E**.3.79E+4, 

4.64E+4, 

5.48E+4, 

6.04E+4.6. 

73E+4. 

7.60E+4. 

8.76E+4, 

TNT 

0180 

1103600. *127400., 

144200., 

166400., 

197100. ,242700., 

267800., 

317600*. 

TNT 

0190 

2391300., 464300., 

512500., 

572500., 

649100. ,750400., 

891100., 

1.10E+6, 

TNT 

0200 

31 

•248E6, 1.444E6* 

1.716E6, 

2.120E6, 

2.342E6/ 

TNT 

0210 

c 

INCIDENT  OVERPRESSURE  CPStl. 

TNT 

0220 

DATA  P/  7800., 

7000., 

6000., 

5000., 

4000., 

3500., 

3000* . 

TNT 

0230 

1 

2500.,  2000., 

1800*, 

1500., 

1200*, 

1000., 

900., 

800*. 

TNT 

0240 

2 

700. ,  600., 

500*, 

400., 

350., 

3001, 

250., 

190., 

TNT 

0250 

3 

180.,  150*, 

120., 

100., 

90., 

80., 

70., 

60*  , 

TNT 

0260 

4 

50.,  40., 

35.. 

30., 

25., 

20., 

18., 

15., 

TNT 

0270 

5 

12.,  10*. 

9.0, 

8.0, 

7.0, 

6.0, 

5.0, 

4.0, 

TNT 

0280 

6 

3.5,  3.0, 

2.5, 

2.0, 

1.8, 

1.5* 

112, 

1.0, 

TNT 

0290 

7 

.90,  .80, 

.70. 

.60, 

.50, 

.40, 

.35, 

.30. 

TNT 

0300 

8 

.25,  .20, 

.18, 

.15, 

.12, 

.10, 

.09, 

.08  , 

TNT 

0310 

9 

.07,  .06, 

.05, 

•  04, 

•  035, 

.03, 

.025, 

.02, 

TNT 

0320 

1 

.018,  .015, 

.012, 

.010, 

.009, 

.008, 

.007. 

.006  • 

TNT 

0330 

1 

.005,  .004, 

.0035, 

.003, 

•0025, 

.002. 

.0018. 

•0015, 

TNT 

0340 

2 

.0012.  .0010, 

.0009, 

.0008, 

.0007, 

.0006, 

•0005, 

•0004, 

TNT 

0350 

1 

.00035.  .0003, 

.00025, 

.0002, 

.00018/ 

TNT 

0360 

c 

SHOCK 

FRONT  ARRIVAL 

TIME  (SECT. 

TNT 

0370 

DATA  TS/  l.E-10. 

•78E-6, 

2.25E-6, 

4.52E-6.8. 

30E-6, 

10.9E-6, 

14.8E-6, 

TNT 

0380 

r< 

..  0E-'‘.-'7.4E-6, 

31.3E-6, 

39.3E-6, 

49.9E-6.60 

•7E-6, 

68.3E-6, 

76.2E-6, 

TNT 

0390 

<£06  •  1 100«E-6  9 

113. E-6, 

146. E-6, 

163. E-6, 186. E-6, 

220. E-6, 

281. E-6, 

TNT 

0400 

3293 #E-6t 343 #E-6 • 

410. E-6, 

479. E-6, 

5 17. E-6, 564. E-6, 

629. E-6, 

714. E-6. 

TNT 

0410 

4840#E-6  9 1  #01  E-*3  u 

1.14E-3, 

1.30E-3, 

1. 51E-3 , l • 

84E-3, 

2.02E-3, 

2.36E-3, 

TNT 

0420 

52 

•84E-3.3.32E-3, 

3.58E-3, 

4.04E-3, 

4.52E-3.5, 

19E-3, 

6.06E-3. 

7.40E-3, 

TNT 

0430 

68 

, 38E-3.9.59E-3, 

I1.3E-3, 

13.9E-3, 

15.2E-3.17 

•8E-3, 

21.7E-3, 

26.0E-3, 

TNT 

0440 

7 

'.0286,  .0315, 

.0355, 

.0407, 

.0479, 

.0587, 

.0659, 

.0758, 

TNT 

0450 

8 

.0884,  .1082, 

.1185, 

.1392, 

.1700, 

.1998, 

.2198. 

.2438, 

TNT 

0460 

9 

.2746,  .3157, 

.3714, 

.4536, 

.5123, 

.5916, 

.6973, 

•8500, 

TNT 

0470 

1 

.9352,  1.106, 

1.355, 

1  .602, 

1.766, 

1.969, 

2.225, 

2.566, 

TNT 

0480 

1 

3.036,  3.735, 

4.229, 

4.631, 

5.783, 

7.123, 

7.861, 

9.324, 

TNT 

0490 

2 

11.49,  13.64, 

15.05, 

16.83  - 

19.06, 

22.04, 

26.18, 

32.32, 

TNT 

0500 

3 

36.66,  42.42, 

50.42, 

62.29, 

68.81/ 

TNT 

0510 

c 

DURATION  OF  POSITIVE 

OVERPRESSURE  I SEC). 

TNT 

0520 

c 

FIRST 

23  VALUES  ARE 

DURATION 

IF  CONTACT  SURFACE  HAD 

NOT  ARRIVED. 

TNT 

0530 

DATA  TP/3.90E-5, 

4.30E-5, 

5.00E-5, 

5.90E-5.7. 

20E-5, 

8.00E-5, 

9.00E-5. 

TNT 

0540 

11 

.02E-4.1.19E-4, 

1.28E-4, 

1.41E-4, 

1.60E-4,!. 

76E-4, 

1.86E-4, 

1.98E-4, 

TNT 

0550 

."ye,”-*"?.  #>  iA***Y®J  M>* ' 
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22.10£-4,2.27E-4,2,45E-4,2.74E-4,2.89E-4,3.0aE-4.3.35E-4,3.80E-4, 
33.88E-4,4.20E-4,4.70E-4,  .5003.  .52E-3,  .55E-3.  .58E-3.  .62E-3. 
4  .68E-3,  .76E-3.  .83E-3.  .90E-3.  .99E-3.1 .12E-3tl,19E-3*1.29E-3, 
51.4ie-3,i.51E-3»1.57E-3»1.65E-3.1.78E-3.1.84E-3.1.96E-3.2.X3E-3. 
62.24E-3,2.36L-3.2.52E-3.2.72E-3,2.82E-3,3.00E-3.3.l7E-3,3.32E-3. 
73.39E-3.3.45E-3.3.54£-3,3.62E-3,3.72E-3»3.84E-3»3.90E-3.3.97E-3. 
84.05E-3.4.15£-3,4.19E-3,4.25E-3.4.34E-3.4.40E-3.4.45E-3,4.49E-3, 
94.52E-3*  0<«  0  •  •  0..  On  0*  •  0**  0.  * 

1  0 • * 0 • » 0 • • 0 • » 0 • *0 • * 0. ♦ 0 • *  0 • » 0 • *0 • * 0 • *0. .0* .0. .0. • 

X  0 • *o. *0 • »o • »o • *o • * 0. *0 •  •  c».o.»o. »o. *o»/ 

C  TIME  (SEC)  BETWEEN  SHOCK  FRONT  ANO  CONTACT  SURFACE  ARRIVAL* 

DATA  TC/  l.E-10,  .IE-6,  .24E-6,  .46E-6.  .B0E-6.1.10E-6,  1.5E-6, 

12.24E-6.  3.3E-6.3.95E-6.  5.2E-6.  7.5E-6 .10.9E-6 . 1 3.2E-6 • 15.9E-6 , 
220o2E-6 * 26*2E-6 »  34.E-6*  49.E-6.  58.E-6.  77. E-6,116. £-6*452. E-6/ 
C 

C  T1»TIME.  T2-TIME  SINCE  SHOCK  ARRIVAL. 

FUNl(Xl,X2)«EXP(AL0GtXl)+FRAC*AL0GlX2/Xl)) 

XIMP1-TCTNT-0. 

C  FIND  DESIRED  POINT  IN  TABLES. 

DO  810  JJ-2.108 
IFIRTNT.LT. R(JJ))  CO  TO  830 
810  CONTINUE  5JJ-108 

C  DESIRED  DATA  LIES  BETWEEN  R(JJ-l)  AND  R(JJ>. 

830  FRAC*ALOG( RTNT/R <  JJ-1 ) ) /ALOG< R ( JJ ) /R ( JJ-1 ) ) 
PTNT«FUN1(P(JJ-1),P(JJ) )  SPSI(1)*PTNT 
C  CALC  ONLY  PSI(l)  IF  L-0  . 

IF(L.EQ.O)  RETURN 
TSTNT-FUN1 (TS(JJ-l).TS(JJ) ) 

TPTNT-FUNl (  TP( J J-l ) .TP  ( JJ ) ) 

IF( JJ.LE.23)  TCTNT-FUN1 (TCI JJ-1) »TC( JJ) ) 

C  FIND  SHAPE  PARAMETER  SIGMA. 

SIGMA-228. /RTNT-. 95 
IF (JJ.LE.23)  SIGMA-228. /65. -.95 
T1 ( 1 ) -TSTNT  ST2(1)«0. 

PS  I  ( KMAX1 )  -0.  STKKMAX1 J «TS TNT+TPTNT  ST2  ( KMAX1 )  -TPTNT 

DO  1030  K-2.KMAXX  S IF ( K.EO.KMAX1 )  GO  TO  1020 

TAU-FLOAT ( K ) /FLOAT ( KMAX 1 1 
PI«(l.-TAU)*EXP(-TAU*(l.*SIGMA/( ,5+TAU) ) ) 

PS  I ( K ) *PTNT*PI  ST2 ( K ) -TAU* TPTNT  J T1 ( K ) -TSTNT+T2 ( K) 

1020  XIMP1-XIMP1+.5*(PSI(K)+PSI(K-1))*(T2(K)-T2 (K-l ) ) 

1030  CONTINUE  JRETURN 
END 


TNT  0560 
TNT  0570 
TNT  0580 
TNT  0590 
TNT  0600 
TNT  0610 
TNT  0620 
TNT  0630 
TNT  0840 
TNT  0650 
TNT  0660 
TNT  0670 
TNT  0680 
TNT  0690 
TNT  0740 
TNT  0750 
TNT  0760 
TNT  0770 
TNT  0780 
TNT  0790 
TNT  0800 
TNT  0810 
TNT  0820 
TNT  0830 
TNT  0840 
TNT  0850 
TNT  0860 
TNT  0870 
TNT  0880 
TNT  0890 
TNT  0940 
TNT  0950 
TNT  0955 
TNT  0960 
TNT  0970 
TNT  0980 
TNT  0990 
TNT  1000 
TNT  1010 
TNT  1020 
TNT  1030 
TNT  1330 


n  n  n  r* 
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SUBROUTINE  AROC(ALTKFT.PAMB.TAMB) 
ALTKFT -ALTITUDE  (KILOFEFT). 
PAHB»AMBIENT  PRESSURE  (PSD. 
TAMB-AMBIENT  TEMPERATURE  (Cl* 


AROCOOIO 

ARDC0020 

ARDC0030 

ARDC0040 

ARDC0050 


ALTZ- ALTKFT* .304BE3  SALTH-6356766.0*ALTZ/ <6356766. O+ALTZ)  ARDC0060 

IF (ALTH.GT, 11000. )  GO  TO  100  ARDC0070 

TEMP-266. 16-0. 0065*ALTH  ARDC0060 

PAMB-14. 696178/ <288. 160/ <288. 160-0. 0065*ALTH)I**5.25612218SGOT0400ARDC0090 
100  IF (ALTH.GT .25000. I  GO  TO  130  ARDC0100 

TEMP-216.66  ARDC0110 

PAMB-3. 28254528/ ( 10.**(0.068463253E-3*(ALTH-1 1000.0) I I  8G0T0  400ARDC0120 

130  IF (ALTH.GT.47000. I  GO  TO  170  AR0C0130 

TEMP-2 16. 66+0 • 003* (ALTH-25000.0)  AR0C0140 

PAMB-O. 36094654/ ((141. 660*3. 0E-3*ALTH|/216. 661**1 1.38826473  AR0C0150 

GO  TO  400  ARDC0160 

170  !F< ALTH.GT. 53000. I  GO  TO  200  AR0C0170 

TEMP-282.66  ARDC0180 

PAMB-0.0174686/(10.**(0.0524926823E-3*(ALTH-47000.0m  SGOTO  400ARDC0190 

200  IF( ALTH.GT .79000. )  GO  TO  230  AR0C0200 

TEMP-282. 66-0. 0045* (ALTH-53000.0)  ARDC0210 

PAMB«8.40408E-3/( (282. 66/TEMP)*»7. 5921761  SGOTO  400  AR0C0220 

230  IF (ALTH.GT. 90000.)  GO  TO  260  ARDC0230 

TEMP-165.66  ARDC0240 

PAMB-1 ,46198E-4*EXP  (-0. 0341647942*1 ALTH-79000. 01/165. 66ISGOTO  400ARDC0250 
260  IF(ALTH.GT. 105000.1  GO  TO  290 

TEMP- 165.66+0.0040* ( ALTH-90000.0 1 
PAMB-1. 5519E-5*(165.66/TEMP)**8. 541198 
290  IF ( ALTH.GT. loOOOO* I  GO  TO  320 
TEMP-225.66+0. 02*(ALTH-105000.0) 

PAMB- 1 . 04442E-6* ( 225.66/TEMP ) **1 . 708239 
320  IF(ALTH.GT. 170000. I  GO  TO  330 

TEMP- 1325. 66+0.01 *(ALTH-160000.0 I 
PAMB«5.14015E-8*( 1325. 66/TEMP)**3. 416479 
350  IF(ALTH.GT. 200000. I  GO  TO  380 

TEMP-1425. 66+0. 005*(ALTH-170000. 01 
PAMB-4.0654E-8*( 1425. 66/TEMP  1**6.832958 
380  TEMP-1575. 66+0. 0035*( ALTH-200000.0) 

PAMB-2.0595E-8*( 1575.66/TEMP) **9.761369 
400  TAMB-TEMP-273.16  SRETURN 
ENO 


ARDC0260 

ARDC0270 

SGOTO 

400 

ARDC0280 

ARDC0290 

ARDC0300 

SGOTO 

400 

AR0C0310 

ARDC0320 

ARDC0330 

SGOTO 

400 

ARDC0340 

ARDC0350 

AROC0360 

SGOTO 

400 

AR0C0370 

ARDC0380 

ARDC0390 

AR000400 

ARDC0410 

